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EXECUTIVE  SUMMARY 


This  final  report  summarizes  work  done  under  AFOSR  Grant  Number  F49620-92-J-0152 
during  the  period  February  1,  1994  to  March  15,  1995  on  the  three-year  research  program 
entitled  “Plasma  Microwave  Electronics:  Studies  of  High  Power  Plasma-Loaded  Backward 
Wave  Oscillators”. 

The  purpose  of  our  research  in  Plasma  Microwave  Electronics  was  to  study  the  funda¬ 
mental  concepts  of  high  power,  efficient  microwave  generation  by  using  relativistic  electron 
beams  in  plasma  loaded  microwave  devices.  The  thrust  was  to  carry  out  an  intensive  experi¬ 
mental  and  theoretical  study  of  improved  plasma-filled  Backward  Wave  Oscillators  (BWOs) 
operating  at  X-Band  (frequency  of  8.5  GHz)  in  order  to  identify  the  nature  of  the  interaction 
in  the  three  component  system  (beam/plasma/e.m.  wave)  and  optimize  it.  We  wished  to  a 
develop  a  systematic  understanding  of  the  operation  of  these  devices  and  their  scaling  with 
parameters  that  would  allow  for  the  future  design  of  larger  diameter  overmoded  devices  to 
obtain  even  higher  output  powers  (1-10  GW). 

During  this  reporting  period,  significant  progress  was  made  both  in  the  experimental  and 
theoretical  aspects  of  the  program. 

A  brief  summary  of  the  work  performed  during  this  reporting  period: 

1)  We  successfully  developed  and  implemented  a  technique  to  determine  the  reflection 
and  the  electromagnetic  field  profile  in  an  empty  cavity  which  is  open,  i.e.  couples  to 
an  antenna. 

2)  A  wide  band  microwave  matching  section  was  designed  and  fabricated.  It  will  be 
installed  between  the  cavity  and  the  output  antenna  in  order  to  improve  the  coupling 
betweeji  the  two.  This  novel  matching  section  is  now  being  tested. 

3)  A  novel,  long  Langmuir  probe  was  developed  in  order  to  accurately  measure  the  back¬ 
ground  plasma  density  and  temperature  inside  the  backward  wave  oscillator  tube.  A 
series  of  experiments  was  completed  successfully  in  order  to  accurately  map  the  density 
profile  in  both  the  axial  and  radial  direction. 

4)  We  successfully  developed  and  implemented  a  microwave  technique  to  measure  the 
dispersion  characteristics  of  smooth  and  corrugated  cavities  loaded  with  time  varying, 
radially  inhomogeneous  plasma.  Using  this  technique,  the  dispersion  characteristics  of 
TMoi  modes  in  a  plasma  loaded  cavity  were  measured  over  a  wide  range  of  background 
plasma  densities,  from  10®  to  10^^  cm“^. 

5)  A  code  was  developed  to  analyze  the  frequency  upshifts  of  TMon  modes  in  a  corrugated 
cavities  due  to  the  presence  of  a  radially  inhomogeneous  plasma  background. 

6)  We  began  a  series  of  hot  tests  of  an  improved  plasma  loaded  backward  wave  oscillator. 
The  improvements  include:  a)  preventing  microwave  breakdown  near  the  collector, 
which  can  cause  microwave  pulse  shortening,  using  an  improved  output  section;  b)  novel 
techniques  for  plasma  injection  and  filling  of  the  tube;  and  c)  electron  gun  protection 
against  adverse  plasma  effects. 
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7 )  We  developed  the  nonlinear  model  of  the  parametric  excitation  of  plasma  mode  in  the 
relativistic  BWO.  This  model  was  used  to  compare  with  experimental  results.  Particu¬ 
lar  attention  was  paid  to  confirmation  of  the  parametric  decay  hypothesis.  Predictions 
of  the  model  with  regard  to  the  output  spectrum  of  the  device,  the  excitation  of  plasma 
waves,  and  the  dependence  on  plasma  density  were  determined.  The  development  of 
the  model  proceeded  through  stages.  Particle  methods  were  introduced.  A  code  was 
developed  in  which  the  plasma  electron  motion  is  one-dimensional  owing  to  the  strong 
applied  magnetic  field.  The  radiation  is  still  modeled  by  an  envelope  equation  allowing 
the  code  to  run  much  faster  than  a  full  PIC  simulation.  Results  of  this  code  indicated 
that  parametric  decay  was  possible  in  principle,  but  not  likely  for  the  parameters  of 
the  experiments. 

8)  We  undertook  an  extensive  comparison  of  the  results  of  model  simulations  of  the  type 
discussed  above,  and  full  PIC  calculations  using  MAGIC.  This  work  is  ongoing.  We 
have  developed  numerical  tools  which  extract  the  EM  fields  from  MAGIC  and  integrate 
the  orbits  of  beam  electrons  in  these  fields.  This  effort  will  ultimately  determine  the 
range  of  validity  of  various  models  of  BWO  operation. 

Publications  and  Conference  Proceedings: 
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Electromagnetic  Properties  of  Open  and 
Closed  Overmoded  Slow- Wave  Resonators  for 
Interaction  with  Relativistic  Electron  Beams 

W.  Main,  Y.  Carmel,  Senior  Member,  IEEE,  K.  Ogura,  J.  Weaver,  G.  S.  Nusinovich,  Senior  Member,  IEEE, 

S.  Kobayashi,  J.  P.  Tate,  J.  Rodgers,  A.  Bromborsky,  Member,  IEEE,  S.  Watanabe,  M.  R.  Amin,  K.  Minami, 
W.  W.  Destler,  Fellow,  IEEE,  and  V.  L.  Granatstein,  Fellow,  IEEE 


Ai^s/racr— Specific  slow  wave  structures  are  needed  in  order 
to  produce  coherent  Cherenkov  radiation  in  overmoded  rela¬ 
tivistic  generators.  The  electromagnetic  characteristics  of  such 
slow  wave,  resonant,  finite  length  structures  commonly  used  in 
relativistic  backward  wave  oscillators  have  been  studied  both 
experimentally  and  theoretically.  In  experiments,  perturbation 
techniques  were  used  to  study  both  the  fundamental  and  higher 
order  symmetric  transverse  magnetic  (TM)  modes.  Finite  length 
effects  lead  to  end  reflections  and  quantization  of  the  wave 
number.  The  effects  of  end  reflections  in  open  slow  wave  struc¬ 
tures  were  found  from  the  spectral  broadening  of  the  discrete 
resonances  of  the  different  axial  modes.  The  measured  axial  and 
radial  field  distributions  are  in  excellent  agreement  with  the 
results  of  a  2-D  code  developed  for  the  calculation  of  the  fields 
in  these  structures. 

Index  Terms--  Periodic  structures,  slow  waves,  dispersion 
curve,  resonators. 

I.  Introduction 

VARIOUS  electrodynamic  structures  capable  of  support¬ 
ing  the  propagation  of  slow  electromagnetic  waves  are 
widely  used  for  microwave  generation  and  for  the  acceleration 
of  charged  particles.  These  slow-wave  structures  are  designed 
to  match  the  phase  velocity  of  the  propagating  electromagnetic 
waves  to  the  speed  of  electrons  in  the  same  structure  in  order 
to  facilitate  an  effective  beam/wave  interaction.  For  interaction 
with  weakly  relativistic  electron  beams  it  is  therefore  necessary 
to  slow  down  the  wave’s  phase  velocity  significantly  while  for 
interaction  with  relativistic  electron  beams  the  phase  velocity 
required  is  only  slightly  smaller  than  the  speed  of  light.  The 
electromagnetic  properties  of  structures  with  very  different 
phase  velocities,  obviously,  would  also  be  expected  to  be  very 
different.  Therefore,  in  spite  of  the  availability  of  detailed 
studies  of  structures  intended  for  operation  in  conventional, 
weakly  relativistic  microwave  tubes  [1],  [2]  there  is  still  a  need 

Manuscript  received  September  30,  1993;  revised  May  30,  1994.  Part  of 
this  work  was  presented  as  an  invited  talk  at  the  IEEE  ICOPS  1993.  This 
work  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research 
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The  authors  are  with  the  Institute  for  Plasma  Research,  University  of 
Maryland,  College  Park,  Maryland  20742-3511  USA,  except  as  noted: 

W.  Main  is  with  ACCURAY,  Inc.,  Santa  Clara,  CA  USA. 

K.  Ogura,  S.  Watanabe,  M.R.  Amin,  and  K.  Minami  are  with  Niigata 
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for  more  detailed  studies  of  slow  wave  structures  intended  for 
operation  with  relativistic  electron  beams. 

A  second  motivation  for  the  work  presented  in  this  paper  is 
the  recent  interest  that  has  developed  around  the  generation 
of  high  power  microwave  radiation  in  devices  employing 
overmoded  slow  wave  structures.  Despite  traditional  concerns 
about  mode  competition  in  overmoded  devices,  recent  experi¬ 
ments  [3],  [4]  have  shown  that  efficient,  high  power  operation 
can  be  achieved  in  such  devices  operating  in  a  single  mode. 
The  present  study,  therefore,  has  been  undertaken  in  part 
to  accurately  determine  the  electromagnetic  characteristics  of 
overmoded  slow  wave  structures  to  aid  in  the  linear  and 
nonlinear  analysis  of  advanced  microwave  sources  employing 
such  circuits. 

Finally,  such  issues  as  finite  length  of  the  periodic  structures 
and  finite  reflections  at  both  ends  should  be  studied  in  more 
detail.  For  example,  in  a  finite  length  structure  the  axial 
wave  numbers  of  the  electromagnetic  modes  are  quantized 
affecting  the  spectral  characteristics  of  device  operation.  Also, 
the  amount  of  reflection  at  both  ends  of  the  structure  affects 
the  quality  factor  of  each  of  the  quantized  axial  modes  in 
the  multiple  resonance  “slow- wave”  cavity  in  a  unique  way. 
This,  in  turn,  strongly  influences  the  interaction  of  the  electrons 
with  the  electromagnetic  waves  associated  with  each  mode. 
Especially  important  is  the  effect  of  reflections  on  the  non¬ 
stationary  operation  of  short  pulse  relativistic  backward  wave 
oscillators  (BWO’s)  [5],  [6]. 

Perturbation  techniques  are  available  for  measuring  the 
spatial  distribution  of  fields  in  resonant  cavities  [1],  [2].  These 
techniques  have  been  primarily  applied  to  the  fundamental 
modes  in  accelerator  cavities,  which  are  electromagnetically 
closed  (shorted)  at  both  ends  and  typically  have  quality  factors 
(Q)  over  4000.  In  contrast,  the  Q  factors  of  spatially  periodic 
slow-wave  structures  intended  for  operation  with  intense  rel¬ 
ativistic  electron  beams  are  much  lower  several  hundred) 
and  the  electromagnetic  properties  of  these  structures  have 
not  been  studied  in  detail  until  recently  [7]-[ll].  Due  to  the 
high  fields  achieved  in  these  devices,  microwave  power  is 
usually  extracted  via  open-ended  matching  sections  that  lower 
the  total  Q  and  can  cause  asymmetric  axial  field  distributions. 
To  accurately  model  this  type  of  device  it  is  important  to  know 
the  reflection  coefficient  at  the  end  of  the  structure  and  the  field 
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distribution.  Tb  date,  parameters  have  usually  been  obtained 
through  numerical  calculations,  which  are  more  complicated  as 
the  Q-factor  decreases.  In  this  study  we  present  experimental 
methods  for  determining  these  parameters  and  for  checking 
the  results  of  numerical  calculations. 

To  properly  address  these  issues  we  have  developed  a 
formalism  describing  the  fields  inside  a  finite  length,  spatially 
periodic  structure  and  methods  for  measuring  the  fields  and 
end  reflections.  For  this  study  we  used  a  periodic  structure 
in  the  shape  of  a  cylindrical  waveguide  with  a  sinusoidally 
varying  conducting  wall.  A  schematic  diagram  of  this  Appara¬ 
tus  appears  in  Fig.  1(a).  The  dispersion  diagram  for  the  first 
six  symmetric  transverse  magnetic  modes  in  an  infinitely  long 
slow  wave  structure  with  the  same  dimensions  is  shown  in 
Fig.  1(b).  The  periodic  nature  of  the  structure  leads  to  a  band 
pass  characteristic  of  each  mode.  The  passbands  are  given 
an  index,  starting  with  1  for  the  lowest  and  increasing  by 
1  for  each  subsequent  passband.  Thus,  the  first  is  referred 
to  as  TMoi,  the  second  as  TM02,  and  the  nth  as  TMon-  In 
this  paper,  we  have  studied  the  fundamental  (TMqi)  and  two 
higher  order  (TM02,  TM03)  transverse  magnetic  modes  in  a 
spatially  periodic,  sinusoidally  corrugated  structure  of  finite 
length.  Preliminary  data  for  TM04  is  also  available. 

This  paper  is  organized  as  follows.  In  Section  II  we  discuss 
the  formation  of  axial  modes  associated  with  each  transverse 
mode  in  a  finite  length  periodic  resonator  and  the  relation  of 
these  modes  to  the  dispersion  diagram  of  the  corresponding 
infinite  length  structure.  This  dispersion  diagram  is  useful  for 
modeling  some  aspects  of  BWO  operation,  such  as  the  relation 
of  beam  energy  to  the  approximate  operating  frequency.  To 
model  the  more  complex  nonlinear  behavior,  it  is  necessary 
to  numerically  simulate  the  operation  of  the  BWO.  For  that 
purpose  one  needs  to  know  the  amplitude  and  phase  of  the 
reflection  coefficient  from  the  structure  ends,  as  well  as  the 
beam-wave  coupling  coefficient.  The  coupling  coefficient  can 
only  be  calculated  once  the  field  profile  in  the  structure  is 
known.  In  Section  III  we  describe  how  the  electromagnetic 
fields  are  calculated  by  expanding  the  fields  in  a  spatially 
harmonic  series.  In  Section  IV  we  show  how  the  resonant  fre¬ 
quencies  as  well  as  the  radial  field  profile  associated  with  each 
axial  mode  of  the  closed  slow-wave  structure  was  measured 
for  the  three  TM  modes  considered.  In  Section  V  we  describe 
how  we  experimentally  determined  the  reflection  coefficient 
at  the  open  end  of  our  slow-wave  resonator.  This  was  done 
by  measuring  the  spectral  width  of  each  resonance,  which 
broadens  as  the  end  of  the  slow-wave  resonator  is  opened.  We 
present  both  experimental  and  numerical  results  for  the  end 
reflection.  Section  VI  summarizes  our  work  and  describes  our 
latest  efforts  to  increase  the  accuracy  of  our  results.  Possible 
extensions  of  these  techniques  to  plasma  filled  systems  are 
also  discussed.  Two  appendices  are  also  attached  to  explain 
the  calculation  of  the  reflection  coefficient  and  the  analytic 
model  used  for  the  closed  cavity  quality  factor. 

II.  Axial  Modes  in  a  Finite  Length  Slow- Wave  Cavity 

The  dispersion  characteristics  of  electromagnetic  modes  in 
an  infinitely  long  spatially  periodic  structure  are  determined 


0  1  2  3  4  5  6 

R=Ro  +  h  cos(27tz/d) 

(a) 


Fig.  1.  (a)  Schematic  diagram  of  a  spatially  periodic  (corrugated  wall) 

structure  shorted  at  both  ends,  (b)  The  calculated  dispersion  diagram  for  six 
symmetric  transverse  magnetic  (TMo*,)  modes  in  an  infinitely  long,  spatially 
periodic  structure  of  Fig.  1(a). 


only  by  the  geometry  of  the  conducting  walls.  A  finite  length 
structure,  however,  can  be  described  by  a  simple  dispersion 
relation  only  when  it  is  well  matched  at  both  ends  (i.e.,  no 
reflections).  Structures  used  in  relativistic  BWO’s  are  usually 
not  well  matched.  First,  relativistic  BWO’s  utilize  a  strong 
reflection  at  the  entrance  of  the  structure  to  prevent  microwave 
propagation  into  the  diode  region  [7],  [8],  [10],  [11].  At  the 
output  end  of  the  structure,  part  of  the  microwave  radiation  is 
reflected  and  part  transmitted.  As  a  result  of  these  reflections 
a  standing  wave  pattern  is  created  leading  to  a  spectrum  of 
axial  modes.  This  effect  was  studied  in  [5]  and  [6]. 

Any  spatially  periodic  structure  with  end  reflections  cOii- 
taining  N  periods  will  support  N  +  1  different  axial  modes 
for  each  transverse  mode  [12].  Each  of  these  axial  modes  is 
characterized  by  a  discrete  frequency  (fr)  and  a  discrete  axial 
wave  number  ((3rX  which  are  located  on  the  dispersion  curve 
of  the  same  transverse  mode  in  an  otherwise  identical  structure 
of  infinite  length.  As  a  consequence  of  the  spatial  periodicity 
of  the  structure,  traveling  waves  can  be  presented  as  a  super 
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position  of  spatial  harmonics.  The  standing  wave  is  formed 
by  a  pair  of  such  traveling  waves  propagating  in  opposite 
directions.  As  an  example,  the  measured  frequencies  and  wave 
numbers  of  the  seven  axial  modes  associated  with  the  lowest 
order  symmetric  transverse  magnetic  mode  (TMoi)  of  the  six 
period  slow-wave  structure  of  Fig.  1  are  shown  in  Fig.  2(a).  In 
Fig.  2(a)  we  show,  for  simplicity,  only  axial  wave  numbers  for 
the  forward  wave  in  the  zero’s  Brillouin  zone.  Bear  in  mind 
that  for  a  given  frequency,  the  total  field  contains  a  set  of  axial 
wave  numbers  corresponding  to  different  spatial  harmonics. 
Throughout  this  paper,  resonance  plots  will  be  presented  as 
frequency  (in  GHz)  versus  the  normalized  wave  number,  0d, 
which  is  equal  to  the  phase  advance  per  structure  period.  Here 
0  is  the  axial  wave  number  and  d  is  the  length  of  the  structure 
period.  The  structure  wall  radius  is  given  by 

Ryj  =  Rq  +  hcos{2'irzfd)  (1) 

where  Ro  is  the  average  radius  and  h  is  the  amplitude  of  the 
wall  corrugations.  Even  though  all  seven  axial  modes  shown 
belong  to  the  same  transverse  magnetic  mode  (TMqi),  they  are 
characterized  by  completely  different  patterns  of  field  lines.  As 
an  example,  the  calculated  electric  field  line  pattern  of  two  of 
the  seven  axial  modes  in  the  slow-wave  structure  of  Fig.  1 
are  shown  in  Fig.  2(b).  The  technique  used  to  calculate  these 
patterns  will  be  described  in  Section  III. 

A  few  features  of  finite  length  stmctures  follow  from 
Fig.  2(a),  which  describes  an  X-band  slow-wave  structure 
having  a  passband  from  7.4  to  8.7  GHz  for  the  fundamental 
symmetric  transverse  magnetic  (TMoi)  mode.  First,  the  dis¬ 
crete  axial  modes  which  are  equally  spaced  in  wave  number 
are  not  equally  spaced  in  frequency  within  the  structure 
passband.  The  mode  separation  varies  between  0.05  and  0.25 
GHz  for  the  structure  studied  in  this  work.  Second,  the  spectral 
resonance  width  (and  thus  the  quality  factor)  of  each  of 
the  axial  modes  is  different  since  the  group  velocity  of  the 
electromagnetic  wave  varies.  This  feature  is  used  to  calculate 
the  end  reflection  of  shorted  and  open  slow-wave  structures, 
as  will  be  shown  in  Section  V. 

In  addition  to  the  dispersion  curve.  Fig.  2(a)  shows  all  seven 
TMoi  resonance  peaks  for  a  closed  six-period  slow-wave 
cavity  as  measured  with  a  microwave  network  analyzer.  For  a 
structure  with  N  periods  of  length  d  the  resonance  condition 
can  be  stated  as  [12] 

Nd  =  (n,./2)(27r//3)  (2) 


where  Ur  is  the  number  of  half  wavelengths  along  the  axis 
of  the  structure  and  {2it/l3)  is  the  axial  wave  length.  The 
resonant  axial  wave  number  fir  is  then  found  from  rir  using 
the  above  equation.  Points  on  the  dispersion  curve  are  found 
by  recording  the  resonant  frequencies  fr  and  associated  axial 
wave  numbers  for  the  set  of  axial  modes  associated  with 
each  transverse  magnetic  (TM)  mode.  The  complete  dispersion 
relation  can  then  be  constructed  from  these  discrete  points  [13] 
or  calculated,  as  will  be  shown  in  Section  III.  Measuring  the 
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Fig.  2.  (a)  The  dispersion  diagram  of  an  X-band  slow  wave  structure  having 

a  passband  of  7.4  to  8.7  GHz  for  the  TMoi  mode  as  calculated  by  Superfish, 
presented  as  a  smooth  line.  Only  half  of  the  first  Brillouin  zone  is  shown.  Also 
shown  is  the  measured  reflection  data  and  all  seven  TMoi  resonance  peaks 
of  the  closed  six  period  slow  wave  cavity  (open  circles),  (b)  The  electric  field 
line  pattern  associated  with  two  of  the  seven  axial  modes  (0,5^/6). 


dispersion  relation  is  therefore  reduced  to  finding  the  set  of 
resonant  frequencies  and  axial  wave  numbers  (/r,  Pr)  for 
desired  TM  mode.  In  practice  it  is  much  easier  to  find  the 
frequencies  than  the  wave  numbers.  For  a  spatially  periodic 
resonator  excited  by  an  ideal  mode  launcher,  the  N  1  axial 
modes  associated  with  a  single  TM  mode  should  be  equally 
spaced  in  wave  number.  Thus  if  the  dispersion  curve,  /(/?),  is 
known  to  be  an  increasing  or  decreasing  monotonic  function  of 
the  wave  number  over  one  Brillouin  zone,  then  it  is  necessary 
to  measure  only  the  resonance  frequencies.  If  the  dispersion 
relation  is  non-monotonic,  then  it  is  necessary  to  measure  both 
/  and  P  to  determine  /(/?).  In  this  study  we  consider  the 
fundamental  (TMoi)  and  the  next  two  higher  order  (TM02  and 
TM03)  transverse  modes,  which  are  monotonic;  however,  the 
techniques  presented  also  apply  to  nonmonotonic  modes.  In 
Section  IV  we  show  how  these  resonances  were  measured. 
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HI.  Numerical  Modeling  of  the 
Slow- Wave  Structure  Fields 

A.  Field  model 

The  first  step  in  the  simulation  of  overmoded  slow  wave 
strucmres  is  to  find  a  way  to  accurately  and  efficiently  calculate 
the  electromagnetic  fields  in  the  structure.  As  higher  modes 
^e  considered,  these  calculations  become  more  complex  and 
It  becomes  also  important  to  experimentally  verify  the  results 
Here  we  present  a  technique  which  shows  good  agreement 
with  expenmental  measurements  for  closed,  spatially^riodic 
structures  at  least  up  to  the  TM03  mode.  ^ 

We  use  a  model  [14H16]  in  which  the  fields  are  expanded 

SoS'^rf  Floquet’s  theorem. 

Solving  the  dispersion  relation  for  the  sinusoidal  boundary 

"  '’y  (')  gives  the  expansion 

coefficients.  Using  Maxwell’s  equation,  we  get  the  following 

expressions  for  the  electromagnetic  field  components,  E  E 
snu 


Z  =  0 


Z  =  L 
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Frequency,  f 
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Rn 
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Er(r,  z,  t)  —  ~iRo  ^2 


^  Ji  ( 

Xn  \  Rq 
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--ieouRo 

n=-oo  / 


{3nZ~ujt) 
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rrSer'o^anTl  f kind  of 

oraer  u  and  1,  respectively, 


xl  =  Rl{w'^/^  -  0l) 


(6) 

where  (ir^=  a  +  n^jd  is  an  integer  and  c  is  the  velocity  of 
ight  in  vacuum.  The  dispersion  relation  is  derived  from^the 
oundaiy  condition  [17]  requiring  that  the  tangential  electric 
field  IS  zero  at  the  wall  r  =  R[z), 

Et{R)  oc  E,{R)  +  Er{R)dR/dz  =  0.  (7) 

The  spatial  Fourier  transform  of  (7)  can  be  expressed  as: 
D-A  =  0 

where  A  is  a  coefficient  vector  in  (3)-(5)  and  is  a  matrix  of 

of  W  anf **  tke  non-trivial  solution 
or  (8)  and  is  given  by 

det  [£>(/?,  w)]  0.  (9) 

In  the  case  of  a  finite  length,  slow-wave  structure,  the  addi¬ 
tional  boundaiy  conditions  at  both  ends  (i:  =  0  and  ^  -  7t 

the  z-axis,  f ,  reflects  at  z  =  I  and  becomes  a  backward 
propagating  wave  5.  Fig.  3(a)  illustrates  this  process,  and 


fin  k  ■  ^  P'''’P=»§ating  forward  along  the  z-axis,  F  reflects  at  -  - 

5  propagating  wave.  <b)  The  location  of  the  F  , 

-o  waves  on  the  dispersion  diagram. 

difgram^  **^°'^*  location  of  these  waves  on  a  dispersii 

Both  forward  and  backward  waves  have  to  satisfy  tl 
dispersion  relation,  ^ 

det[D(/3r,u;)]  =  0.  (I, 

det[D(/3s,cj)]  =  0.  (1 

Electromagnetic  quantities  must  also  be  single  valued  at 
given  position  and  time  at  the  resonance, 

P(o{^(3f-3b)L}  _  2 

Here,  p  is  a  round  trip  reflection  coefficient.  With  the  assume 
ion  of  a  loss  less  cavity  that  is  completely  shorted  at  both  ends 
we  may  write  p  =  1  and  =  -/?e.  Hence,  (12)  becomes 

Pf  —  t^N/L  (^12 

where  N  is  an  integer.  Except  for  the  propagation  direction, 
the  backward  wave  is  identical  to  the  forward  wave;  that  is 
It  satisfies  the  same  radial  boundary  condiUon  and  has  the 
s^e  energy.  Therefore,  in  (3)-(5),  the  relationship  between 
ffie  coefficients  A^  of  backward  wave  and  those  of  the 
rorward  wave  may  be  written  as  ” 

>4^  —  aB 

-  A_^,  (24) 

By  summing  the  forward  and  backward  waves,  the  electro¬ 
magnetic  field  components  in  the  cavity  can  be  written 

«=  /  \ 

£.(r,.,()  =  c«s((3;..)  (15) 
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Er{r,  zA)  Roe  ^  f  sin(/3n^) 

°  (16) 

H0{r.  2,  t)  =  ~  *^1  ( cos(/3„2). 

n=-oo  \«0  / 

(17) 

In  our  calculations,  the  rank  of  the  matrix  was  truncated  at  9 
(-4  <  n  <  4)  [15],  [16].  Once  the  dispersion  relation,  (9),  is 
solved  numerically,  the  ratios  of  the  coefficients,  An/Ao^  are 
determined  from  (8).  With  a  known  set  of  An/Ao^  the  field 
components  Ez^  Er*  and  Htf,  are  calculated  from  (15)“(17). 

B.  Perturbation  Theory 

Perturbation  theory  indicates  that  slight  changes  in  the  shape 
of  a  resonant  cavity  can  affect  the  resonant  frequency.  Each 
resonance  (associated  with  an  axial  mode)  is  frequency  shifted 
by  a  perturbing  object  of  volume  AV^  by  a  different  amount. 
The  general  relation  for  the  relative  frequency  shift  is  given 
by  [2] 

"  jy  iitH  H’ +  cB  ■  E‘)dV 

where  the  integral  in  the  numerator  is  evaluated  over  the 
perturbing  object  volume  Al^  and  the  integral  in  the  denom¬ 
inator  is  evaluated  over  the  whole  cavity.  This  assumes  that 
the  perturbing  object  is  small  compared  to  the  amplitude  and 
length  of  the  wall  ripple.  For  the  case  of  a  sphere  [2]  the 
above  equation  is  geometrically  corrected  by  multiplying  the 
first  term  in  the  numerator  by  3/2  and  the  second  term  in  the 
numerator  by  3.  For  a  small  metal  bead  of  radius  ro  the  field 
is  almost  constant  over  so  the  integration  is  unnecessary. 
Thus  the  equation  commonly  used  is  [1] 

(19) 

where  En  and  Ho  are  the  field  amplitudes  normalized  so  that 
the  integral  of  or  E^  over  the  cavity  is  unity.  Since  the 
only  field  component  on  axis  is  the  simplest  measurement 
would  be  to  find  Pr  by  perturbing  the  field  on  axis.  In  our 
experiment,  a  bead  was  placed  0.54  cm  away  from  the  axis 
because  of  the  difficulty  associated  with  having  both  the  bead 
and  the  antenna  on  axis.  At  this  radius  it  was  still  possible 
to  identify  Pz  fot  the  TMoi  modes  by  inspection  of  the 
axial  profile  of  the  frequency  shift.  Once  the  cavity  fields  are 
calculated,  ( 19)  can  be  used  to  find  the  associated  perturbation. 
This  frequency  shift  can  then  be  directly  compared  with 
experimental  measurements. 

IV.  Results  of  the  Field  Model 
AND  Experimental  Verification 

There  are  two  levels  at  which  we  can  compare  the  field 
model  of  Section  III  with  experiment.  The  first  and  simplest 
level  is  to  compare  the  measured  and  calculated  resonant 
frequencies  associated  with  the  axial  modes.  The  second 
is  to  compare  the  field  distribution  (or  the  related  spatial 
distribution  of  frequency  shifts)  of  these  modes  using  the 
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Fig.  4.  Cold  test  system  used  for  measuring  the  electromagnetic  character¬ 
istics  of  spatially  periodic  structures.  The  perturbing  bead  could  be  translated 
both  axially  and  radially. 

spatial  perturbation  technique  described  in  Section  III.  In 
this  section  we  will  present  the  comparison  on  both  levels. 
Fig.  4  shows  the  measurement  system,  including  the  bead-pull 
apparatus  to  be  described  later. 

A.  Resonant  Frequencies 

Experimentally,  the  modes  of  our  slow-wave  structure  were 
excited  using  a  small  Hertzian  antenna  on  the  axis  which 
could  be  moved  axially  in  and  out  of  the  cavity  to  adjust 
the  degree  of  coupling  and  to  enable  calibration.  A  network 
analyzer  was  used  to  measure  the  microwave  reflection  from 
the  cavity  over  the  desired  frequency  range  (5u  single  port 
measurement).  The  resonances  appeared  as  narrow  spikes  at 
frequencies  where  the  magnitude  of  the  reflection  was  reduced. 

The  experimental  results  appear  as  plots  of  the  resonant 
frequency  (fr)  versus  the  normalized  axial  wave  number  (0d). 
All  measurements  were  compared  with  numerical  calculation. 
The  measurements  agree  very  well  with  the  boundary  specific 
calculation  of  Section  III  and  the  general  boundary  calculation 
using  the  2-D  electromagnetic  code  Superfish  [19].  A  graph¬ 
ical  comparison  of  the  experimental  results  with  Superfish 
calculations  is  shown  in  Fig.  2(a).  This  figure  shows  the 
dispersion  curve  and  reflection  data  with  a  common  frequency 
axis.  The  Superfish  results,  which  were  actually  calculated  at 
seven  points,  are  converted  to  form  a  smooth  curve  to  aid 
in  comparison  with  experiment.  Similar  plots  comparing  the 
measured  and  calculated  results  [16],  [20]  for  the  TM02  <ttid 
TMo3  modes  are  shown  in  Fig.  5(a)  and  5(b),  respectively.  The 
measured  frequencies  of  the  TMoi  resonances  are,  on  average, 
0.17%  lower  than  our  calculations  and  0.03%  higher  than 
Superfish.  For  the  TM02  TM03  modes  the  experimental 
data  are,  respectively,  0.05%  and  0.37%  lower  than  our 
calculations. 

B.  Field  Distribution 

To  perturb  the  cavity  a  small  object  had  to  be  placed  at 
known  locations  within  the  structure  and  the  resulting  shift  in 
the  frequency  of  a  given  resonant  mode  was  measured.  Fig.  4 
shows  a  schematic  diagram  of  the  apparatus  and  Table  I  gives 
a  few  critical  dimensions. 
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TABLE  I 

Dimensions  of  Slow- Wave  Structure. 

Structure  Period  1.67  cm 

Total  Length  10  cm  (6  periods) 

Radius  [cm]  1.5  +  .41  sin(3.7z) 

Bead  Diameter  0.24  cm 

Radial  Position  of  String _ 0.54  cm _ 


- Numerical 


(a) 


3  0  -3  -6  -9  -12 
Reflection  (dB) 


Normalized  Axial  Wavenumber,  ^ 
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Fig.  5.  (a)  The  dispersion  diagram  of  an  X-band  slow  wave  structure  having 

a  passband  of  7.4  to  8.7  GHz  for  the  TM02  mode  (smooth  line).  Also  shown 
is  the  measured  reflection  data  for  the  TM02  resonance  peaks  of  the  closed 
six  period  slow  wave  cavity  (open  circles),  (b)  Same  as  (a),  but  for  the 
TMosmode. 


°  Axial  Radial  ^ 
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Fig.  6.  Field  perturbation  results  in  the  shorted  structure  for  three  axial 
modes  of  the  TMoi  transverse  mode,  (a)  Plots  of  the  electric  field  lines, 
axial  and  radial  frequency  perturbation  for  7r/2  axial  mode;  solid  lines  - 
calculation  based  on  equation  18;  open  circles  -  experiment,  (b)  Same,  for 
27r/3.  (c)  Same,  for  57r/6. 


O  Momired 


Position  (cm)  Position  (cm) 


(b) 

Fig.  7.  Field  perturbation  results  of  the  shorted  structure  for  one  axial  mode 
of  high  order  transverse  modes,  (a)  Plots  of  the  electric  field  lines,  axial  and 
radial  frequency  perturbation  for  57r/6  axial  mode  of  the  TM02  mode;  solid 
lines  -  calculation  based  on  (18);  open  circles  -  experiment,  (b)  Same,  for  tt 
mode  of  the  TM03  transverse  mode. 


An  aluminum  bead  was  suspended  in  the  cavity  on  a  nylon 
thread  via  one  of  four  sets  of  access  holes:  one  set  parallel  to 
the  cavity  axis  and  three  along  the  diameter  at  axial  positions 
corresponding  to  the  maximum,  average,  and  minimum  radii 
of  the  slow  wave  structure.  Bead  movement  was  regulated 
by  tying  one  end  of  the  thread  to  a  dial  caliper  to  measure 
position  and  the  other  end  to  a  small  weight  to  provide  tension. 
An  HP  network  analyzer  measured  the  forward  and  reflected 
power  versus  frequency  (5ii)  from  a  coax  Hertzian  element 
inserted  on  axis  through  the  cavity  end  plate.  For  each  mode, 
the  resonant  frequency  was  recorded  as  a  function  of  bead 
position.  The  frequency  shift  at  a  given  point  is  then  the 
difference  between  this  frequency  and  the  frequency  in  an 
unperturbed  cavity.  It  should  be  noted  that  the  effect  of  the 
nylon  thread  was  negligible  (less  than  500  kHz  in  all  cases). 
Both  axial  and  radial  dependence  of  frequency  shifts  were 
measured.  The  first  shift  was  measured  with  the  bead  translated 
along  a  path  parallel  to  the  cavity  axis  at  a  radial  position  of 
0.54  cm.  The  second  shift  corresponded  to  a  radial  path  at  a 
position  of  maximum  cavity  diameter  (see  Fig.  4). 

In  general,  measurements  and  calculations  can  be  performed 
for  each  of  the  seven  axial  modes  associated  with  each  of  the 


three  transverse  modes  (TMqi,  TM02,  TM03).  In  Fig.  6  and 
7  we  present  selected  results.  All  seven  axial  modes  of  the 
TMoi  group  of  modes  and  some  of  the  axial  modes  of  the 
TMo2»  TMo3,  and  TM04  have  been  studied  numerically  and 
experimentally.  The  results  are  shown  in  Fig.  6(a),  (b),  (c) 
for  the  TMoi  7r/2,  TMoi  27r/3,  and  TMoi  57r/6  modes,  re¬ 
spectively.  The  experimental  results  appear  as  circles  overlaid 
on  solid  lines  that  represent  the  numerical  results.  To  help 
visualize  the  cavity  fields,  a  plot  of  the  electric  field  lines 
for  each  axial  mode  is  displayed  beside  the  corresponding 
data.  Similarly,  Fig.  7(a)  displays  the  57r/6  mode  of  the  TM02 
series  and  Fig.  7(b)  shows  the  tt  mode  of  the  TM03  series. 
The  normalized  wave  number,  (3d,  can  be  found  from  these 
results  by  simply  counting  the  number  of  local  maxima  along 
the  length  of  the  resonator  and  multiplying  the  number  by  7r/6. 

V,  Reflection  and  Q-Factor  of 
AN  Open  Slow-Wave  Structure 

In  most  low  power  linear  microwave  devices  the  output 
power  is  extracted  radially  at  one  end  of  the  interaction 
region.  Microwave  mode  converters  change  the  cylindrical 
mode  produced  in  the  interaction  region  into  the  fundamental 
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radiative  power  loss  through  the  output  aperture  of  the  cavity, 
whereas  the  power  loss  to  the  walls  is  related  to  the  ohmic 
quality  factor  of  the  cavity  (Qohm)- 

To  find  Qd  it  is  necessary  to  separate  the  ohmic  losses 
(<3ohm)  from  the  total  circuit  losses  (Qt)  which  include  the 
radiative  losses  of  the  cavity.  Since  the  Q-factor  is  inversely 
related  to  power  loss,  Qd  can  be  found  from  Qt  and  Qohm 
by  the  relation 


TMoi  7t/2Mode 


Fig.  8.  (a).  A  schematic  diagram  of  an  open  slow  wave  cavity  —  a  way  for 

axial  power  extraction  in  a  relativistic  backward  wave  oscillator.  In  this  work 
a  =  1.02  cm.  h  =  7.62  cm.  r  =  4.75  cm.  (b).  Demonstration  of  resonance 
broadening.  The  tr/a  resonance  of  the  TMoi  mode  is  broadened  when  the 
structure  is  open  (dotted  line)  as  compared  to  the  shorted  case  (solid).  This 
effect  is  used  to  calculate  the  reflection  at  the  end  of  the  structure. 


mode  of  a  waveguide  which  can  transport  the  energy  to  the 
load  or  antenna.  At  power  levels  of  tens  of  MW,  the  complex 
geometry  of  the  mode  converter  can  lead  to  breakdown  even 
in  vacuum.  To  avoid  this  problem  high  power  tube  designs 
eliminate  the  mode  converter  and  allow  the  output  power 
to  exit  the  slow-wave  structure  along  the  beam  axis.  Then 
the  beam  is  separated  from  the  microwave  by  tapering  the 
magnetic  guide  field  and  thus  dumping  the  beam  into  the  wall 
of  the  waveguide.  For  these  discussions  this  type  of  cavity  is 
labelled  “open”  [23].  A  diagram  of  a  practical  open  cavity 
used  in  high  power  BWO  experiments  [8],  [11]  is  shown  in 
Fig.  8(a). 

When  modeling  a  BWO  it  is  important  to  know  the  end  re¬ 
flection  of  the  structure.  End  reflections  affect  both  the  starting 
conditions  and  saturation  effects  in  the  device  [5],  [6].  Several 
numerical  codes  can  predict  the  end  reflection  coefficient  (p) 
for  a  given  geometry,  however,  it  is  often  difficult  to  accurately 
calculate  end  reflections  due  to  the  complex  geometry  of  such 
systems  and  the  often  unknown  surface  finish  conditions  used. 
For  these  reasons  it  is  important  to  perform  experimental 
measurements.  The  most  direct  method  of  measuring  p  is  to 
launch  the  desired  mode  down  the  structure  and  measure  the 
waves  which  are  transmitted  through  or  reflected  from  the 
end  of  the  structure.  This  technique  requires  both  launching 
the  desired  mode  and  completely  absorbing  the  reflected  wave 
inside  the  structure.  Either  of  these  would  be  very  difficult  in 
practice.  It  is  much  easier  to  find  p  indirectly  by  exploiting  its 
relation  to  the  diffractive  quality  factor  of  the  structure.  The 
diffractive  quality  factor  of  a  structure  (Qd)  is  related  to  the 


Q,  =  1/{1/Qt  -  1/Qohm).  (20) 

Both  Qt  and  Qohm  can  be  measured  directly.  The  measured 
quality  factor  of  the  system  is  Qt  when  the  structure  is  open, 
and  <3ohm  when  the  diffractive  loss  is  eliminated  by  placing  a 
short  at  the  output  side  of  the  structure.  Fig.  8(b)  demonstrates 
this  effect  by  showing  how  the  7r/2  resonance  of  the  TMoi 
mode  is  broadened  when  the  end  of  the  structure  is  opened. 


A.  Experimental  measurement  of  Q 

It  is  useful  to  give  some  details  of  our  Q  measurement  here. 
We  measured  the  voltage  standing  wave  ratio  (VSWR)  of 
the  cavity  across  the  desired  resonance  or  resonances  using 
a  HP8510C  vector  network  analyzer  (we  used  a  one-port 
measurement  for  simplicity).  From  this  data  we  found  the 
minimum  VSWR  (at  resonance)  and  the  maximum  VSWR 
(between  resonances).  From  these  two  values  we  used  a 
technique  [21]  to  find  the  VSWR  level  at  which  to  read  the 
upper  and  lower  frequencies,  fi  and  /2,  of  each  resonance. 
From  these  frequencies,  and  the  resonant  frequency  fr  we  can 
calculate  the  cavity  Q  using  the  relation 

Q  =  fr/{fl-f2)-  (21) 

With  the  additional  knowledge  of  whether  the  cavity  was 
under-coupled  (less  than  optimum  coupling)  or  over-coupled 
(more  than  optimum  coupling)  a  correction  factor  for  the  Q 
was  found  that  eliminated  the  loading  effect  of  our  mea¬ 
surement  apparatus  and  wave  launcher.  The  cavity  coupling 
could  be  adjusted  by  moving  the  launcher  in  and  out  of  the 
cavity.  Using  the  minimum  insertion  which  would  give  a  clear 
resonance  we  were  certain  to  be  in  the  under-coupled  regime. 
The  cavity  was  made  of  brass. 

Using  these  techniques  we  measured  the  Q’s  of  the  cavity 
modes  in  both  open-ended  and  close-ended  configurations. 
As  a  check  we  compared  the  measured  Q’s  for  the  shorted 
configuration  with  an  analytic  calculation  which  is  described 
in  detail  in  Appendix  II.  These  results  appear  in  Fig.  9. 
The  agreement  of  the  measured  and  calculated  Q’s  is  very 
good  when  we  use  a  surface  roughness  factor  of  1.7.  This 
factor  effectively  increases  the  skin  depth  in  the  calculation 
to  compensate  for  increased  field  penetration  due  to  surface 
imperfections  (i.e.  fabrication,  oxidation). 

We  use  (20)  to  calculate  Qd  from  the  measured  values  of 
Qt  and  Qohm.  The  only  remaining  step  is  to  calculate  the  end 
reflection,  p,  from  Qd-  Appendix  I  presents  a  derivation  of  the 
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Fig.  9.  Calculation  (open  squares)  and  measurement  (bars)  of  the  wall  losses 
for  various  axial  modes  associated  with  the  TMoi  mode  of  the  shorted  periodic 
cavity.  The  surface  finish  factor  is  1.7. 


relation  of  these  two  parameters.  The  result  is 


P  = 


1  -  {o^/Qd) 
1  -h  [a/Qd] 


(22) 


where  a  =  (u;/c)^(L//3).  Using  this  relation  the  measured 
values  of  p  were  found.  These  values  of  p  for  the  TMqi  axial 
modes  were  compared  with  a  set  of  calculations  performed 
by  A.  Bromborsky  [22].  The  results  appear  in  Fig.  10(b). 
The  measured  reflection  coefficients  are  up  to  15%  larger 
than  the  calculated  ones.  This  is  a  significant  difference  since 
reflections  are  important  for  both  the  linear  characteristics 
(starting  conditions)  and  for  simulations  of  non-stationary 
operation  of  pulsed  relativistic  backward  wave  tubes.  For 
example,  a  15%  change  in  the  reflection  coefficient  can  cause 
the  starting  current  to  double. 


VI.  Summary  and  Conclusions 

The  electromagnetic  properties  of  slow- wave  resonators 
for  high-power  relativistic  backward  wave  tubes  have  been 
studied  both  theoretically  and  experimentally.  The  agreement 
between  the  measured  and  calculated  frequencies  and  spatial 
field  distributions  is  very  good  (the  discrepancy  is  less  than 
0.2%)  for  the  first  three  symmetric  transverse  magnetic  modes 
(TMop^)  with  the  index  p  from  1  to  3  and  various  axial  indices 
i  (from  0  to  where  N  is  the  number  of  periods  of  a  slow- 
wave  resonator).  Preliminary  data  is  also  available  for  the 
TMo4  mode,  which  is  characterized  by  a  very  flat  dispersion 
curve  (see  Fig.  1). 

We  have  also  developed  a  simple  method  to  measure  the  end 
reflections  in  open  slow-wave  resonators  by  measuring  their 
Q-factors.  This  measurement  technique  gives  results  for  end 
reflections  which  are  up  to  15%  higher  than  those  that  follow 
from  computer  code  simulations  [22]. 

Applying  these  techniques  is  especially  important  to  high- 
power  relativistic  BWO’s  with  overmoded  structures  because 
their  dispersion  diagram  may  not  be  a  monotonic  function 
of  the  axial  wave  number.  For  this  reason  it  is  necessary  to 
measure  both  the  frequency  and  the  axial  wave  number,  as 
was  done  in  this  work,  in  order  to  identify  the  transverse  and 
axial  indices  of  the  operating  mode. 
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Fig.  10.  (a)  Measurement  of  the  diffractive  quality  factor  of  the  open 

periodic  cavity  for  various  axial  modes  associated  with  the  TMoi  mode,  (b) 
The  measured  reflection  coefficient  for  various  axial  modes,  derived  from  the 
diffraction  quality  factor  using  (22). 


The  results  obtained  are  important  in  the  design  of  future 
high-power  relativistic  backward  wave  oscillators,  relativistic 
travelling  wave  tubes,  overmoded  multiwave  Cherenkov  and 
multiwave  diffraction  generators.  We  also  plan  to  extend 
the  technique  described  here  to  the  study  of  plasma  filled 
relativistic  slow  wave  devices  [24]. 
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Appendix  I 

Dependence  of  Cavity  End-Reflectionon 
Diffractive  Quality  Factor 

We  consider  a  cavity  consisting  of  a  section  of  waveguide 
with  a  short  at  one  end  and  a  “leaky”  aperture  (antenna)  at 
the  other.  We  find  the  dependence  of  the  reflection  from  this 
aperture,  p,  on  the  diffractive  quality  factor  of  the  cavity  Qd 
in  two  steps.  First,  we  relate  Qd  to  the  ratio  of  circulating 
power  to  power  loss,  and  then  we  relate  this  ratio  to  the  end- 
reflection.  We  begin  with  a  description  of  Qd  and  how  it  relates 
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to  the  total  and  ohmic  quality  factors.  These  three  Q-factors 
are  defined  as  follows: 

Qohm  =  Q’factor  related  to  RF  power  lost  in  the  cavity  wall 
by  resistive  heating 

Qt  ~  QTactor  related  to  all  power  loss  mechanisms  (ohmic 
and  radiation) 

Qd  =  Q-factor  related  to  power  lost  by  diffraction  through 
an  aperture  in  the  cavity. 

By  definition,  the  Q-factor  is  related  to  the  cavity  power 
loss,  Floss,  hy  the  equation 

Q  =  a;(C//Fioss) 

where  uj  is  the  resonant  frequency  and  U  is  the  electromagnetic 
energy  stored  in  the  cavity.  Because  Q  is  inversely  related  to 
power  loss,  and  the  total  power  loss  is  equal  to  the  sum  of 
the  diffractive  and  ohmic  losses,  the  three  Q-factors  have  the 
relation 


l/Qt  =  1/Qohm  +  l/Qd.  (23) 

An  experimental  measurement  of  the  open  cavity  would  give 
the  value  for  Qt-  If  the  aperture  was  shorted  then  the  measure¬ 
ment  would  give  Qohm-  Thus,  using  this  equation,  all  three  Qs 
can  be  experimentally  determined. 

We  now  proceed  to  relate  Qd  to  the  end  reflection.  To 
accomplish  the  first  step  we  use  the  relation  between  power 
and  energy  in  an  empty  waveguide 


Thus,  substituting  into  the  equation  for  Qd 


Qd  =  oc 


i  +  \p? 

1  -  IpP 


where 

Inverting  this  expression  in  terms  of  p,  we  find 


P  = 


1  -  ja/Qd) 

1  +  {n/Qd)' 


(30) 


(31) 


(32) 


This  is  the  relation  we  used  to  find  the  cavity  end-reflection 
in  Section  V  of  this  paper  (22). 


Appendix  II 

Ohmic  Losses  in  Resonators  With  Corrugated  Walls 

Consider  a  resonator  formed  by  the  part  of  a  cylindrical 
waveguide  with  corrugated  walls  which  is  bounded  with  two 
end  walls.  In  general,  the  ohmic  quality  factor  (Q)  of  any 
resonator  can  be  defined  as 


W  IH^PdV 

Wohm  —  -f77  ^  .  — =: 


(33) 


P/u  =  {p/uj)c^  =  Vgr  (24) 

where  u  is  the  energy  per  unit  length  in  the  waveguide  (U/L) 
and  L’gr  the  group  velocity.  Using  this  relation  the  definition 
of  quality  factor  becomes 

Qd  -  {uj/c)HL/f3){P,irc/Pd)  (25) 

where  Pc\tc  now  represents  the  power  circulating  in  the  cavity. 
To  relate  Pc\rc/Pd  to  the  end-reflection  p,  we  consider  the 
forward  and  reverse  waves  in  the  cavity  with  associated  power 
F+  and  F_,  and  the  power  lost  through  the  cavity  aperture 
Pd.  By  definition 


P-/P+  =  IpI'  (26) 

and 

F+  -  F_  =  Pd  (27) 


so 


Here,  the  bars  imply  the  averaging  over  the  wave  time 
period  27r/u;,  W  is  the  microwave  energy  stored  in  the 
resonator,  Wohm  is  the  microwave  energy  concentrated  in  a 
skin-layer  of  the  metallic  walls.  Kkin  is  the  volume  associated 
with  the  skin  depth.  The  coefficient  2  reflects  the  fact  that 
inside  a  resonator  the  electric  field  amplitude  is  equal  to  the 
magnetic  field  amplitude  while  in  a  metallic  wall  (with  the 
finite  but  very  large  conductivity,  a)  only  the  magnetic  field 
is  significant  since  ~  (i/y/a)H^. 

We  will  restrict  our  consideration  to  the  simplest  but  most 
important  case  of  symmetric  TMop^  modes.  The  magnetic  field 
of  such  modes  has  only  one,  non-zero  component  directed 
along  azimuthal  coordinate  6  In  a  corrugated  structure 

with  strong  end  reflections  this  field,  excited  at  the  given 
frequency  to,  can  be  presented  as 

-h  complex  conjugate. 


where 


CX3  , 

Hd,^ieowRo  T  —.h{gr.r)P^-^ 
Xn 


(34) 


PdlP+  =  1  -  IpI'.  (28) 

By  eliminating  F+  and  F_  we  relate  the  power  ratio  to  the 
reflection  coefficient 

Fcirc  _  P+  P-  _  ^  ^ 

Pd  ■  Pd  "1-IpP* 


Pn  =  I'Tt/L  +  Tip 

Here  I'k/L  is  the  axial  wave  number  of  the  field,  the  wave 
number  P  is  determined  by  the  period  of  corrugation  d: 

=  27r/d,  the  transverse  wave  number  gn  is  equal  to  Qn  = 
[(a;/c)2-/32ji/2  jhis  is  the  form  given  by  (5).  For  slow  spatial 
harmonics  0)^  therefore, 
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where  Below  we  will  take  into  account  only  zero 

and  ±  first  spatial  harmonics  supposing  that  amplitudes  of  all 
other  harmonics  are  negligibly  small. 

Let  us  consider  periodical  corrugation  of  the  resonator  wall, 

R  =  Rq  ^  h  '  cosfiz  (35) 


expressions: 

^ohm,corr  =  7r6gjfeL  RoJiiXp) 

1  +  alql  +  +  a-ig-i)  +  |, 

(40) 


and  suppose  that  the  height  of  these  ripples,  h,  is  much  smaller 
than  the  wavelength,  A.  In  such  a  case  of  small  ripples  it  is 
possible  to  replace  the  boundary  condition  for  the  field  at  the 
real  corrugated  wall  by  an  approximate  boundary  condition 
at  the  wall  of  the  cylindrical  waveguide  of  radius  Ro  [25]. 
From  this  condition  one  can  find  the  relation  between  spatial 
harmonic  amplitudes  [26]: 

=  _.j!l .  go  T  ■  p  _  MgoRo) 

^0  2  p±i  Io{p±iRo) 

In  the  case  of  arbitrary  relations  between  h  and  A  the  ratio  of 
harmonic  amplitudes  can  be  found  only  numerically. 

In  addition  to  the  assumptions  given  above  let  us  also 
assume  that  the  structure  consists  of  an  integer  number  of 
periods:  L  =  N  '  d  where  N  is  an  integer.  Then,  for  the 
microwave  energy  stored  in  the  resonator  one  can  obtain  the 
following  formula: 


2 

(l  +  OLi<Pl  -1- 

(37) 

where  Xp  =  90^0  is  the  pth  root  of  the  equation  Jo(x)  =  0^ 
which  is  the  boundary  condition  for  the  TMop^  mode  in  a 
waveguide  of  the  constant  radius  Rq, 


W  =  irLRi 


£0Wi?0— )  Jl(Xp) 
Xo  / 


a±i 
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fh 
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0±i  [a  7o(p±ii?o)72(p±i7?o)]  • 

(38) 

The  last  two  terms  in  (II. 5)  correspond  to  the  ±  first  spatial 
harmonics. 

In  a  similar  manner  one  can  find  the  microwave  energy 
stored  in  a  skin-layer  of  the  depth  5sk^  Note  that  due  to  a 
certain  roughness  of  the  wall,  a  realistic  skin  depth  is  1.5-2 
times  larger  [27]  than  the  theoretical  value.  Correspondingly, 
the  ohmic  Q-factor  of  the  closed  cavity 


^closed 

^ohm 


w 


^  ohm  ,corr 


(39) 


where  l^ohm,corr  is  the  average  microwave  energy  concen¬ 
trated  in  a  skin-  layer  of  the  corrugated  wall  of  a  waveguide 
section,  and  I^ohm,e.w.  is  the  microwave  energy  concentrated 
in  a  skin-layer  of  end  walls.  Denoting  the  skin-depth  by 
Sgk  one  can  derive  for  FFohm,corr^  ^ohm,e.w.  the  following 


W^ohm,e.w.  =  ‘2'K6sk(^^O^Ro^^  ^0‘^l(Xp) 

|l  4*  al<f>i  +  +  Aa\  0_^2li7v/L)  ' 

+  4q:_i  ^_2(i7^/L)  ^  P-lQl)^ 

Here  q±i  denotes  Iiip±iRo)/Io{p±iRo)- 
Substituting  equations  (11,5),  (II.6),  (IL8)  and  (II.9)  into 
(IL7)  one  can  find  the  ohmic  Q-factor  of  the  closed  cavity.  In 
microwave  generation  experiments  [3,4, 6, 7]  structures  without 
end  walls  are  used.  Correspondingly,  when  reflections  of  the 
microwave  power  at  both  ends  (^  =  0  and  2  =  L)  are  large 
enough,  the  ohmic  Q-factor  of  the  resonator  can  be  estimated 
by  the  following  equation 


gopen  ^ 
^  ohm 


w 


ohm 


,corr 


(42) 


where  W  and  H^ohm.corr  are  defined  by  Eqs.  (II.5)  and  (II.8), 
respectively.  Equation  (II.7)  was  used  in  Section  V  (and  Fig.  9) 
to  calculate  the  ohmic  losses  in  a  shorted  cavity. 

Two  limiting  cases  should  be  considered  separately.  In 
case  of  the  0-mode,  we  have  /  =  =  P^nPi  — 

p-i,qi  =  q-i^ai  —  a^i  and  —  0-i,  which  means 
that  the  first  and  minus  first  harmonics  degenerate  to  each 
other  except  for  difference  in  propagating  direction.  Thus  we 
have  to  recalculate  our  formulas  from  the  first  step  by  taking 
the  degeneracy  into  consideration.  This  leads  to  simplified 
expressions  for  terms  defined  by  eqs.(II.8),(II.9).  We  do  not 
present  here  these  formulas  because  it  was  found  that  for 
parameters  of  our  slow-wave  structure  the  amplitudes  of  ± 
first  harmonics  are  small  (Ai  —  A^i  —  0.02Ao),  and  therefore 
these  harmonics  can  be  ignored  in  calculation.  This  fact  has 
been  used  to  calculate  the  Q-value  of  the  0-mode  in  Fig.  6. 

In  case  of  7r-mode,  a  zero  harmonic  of  forward  wave  and 
minus  first  harmonic  of  backward  wave  have  the  same  axial 
wave  number  as  well  as  a  zero  harmonic  of  backward  wave 
and  plus  first  harmonic  of  forward  wave.  They  propagate  in 
opposite  directions  and  thus  form  a  standing  wave  together. 
In  other  words,  these  harmonics  degenerate  to  each  other  in 
terms  of  dependence  of  the  radial  coordinate  and  therefore  all 
formulas  must  account  for  this  effect.  We  do  not  present  here 
these  formulas  recalculated,  because  it  was  found  that  in  our 
experiment  the  intersection  between  the  light-line,  a;  = 
and  the  dispersion  curve  is  close  to  the  point  of  the  7r-mode, 
which  means  \  go  0,  where 

2  2  /u;\2  ^2 

=  +->0. 

This  condition  \  go  P  simplifies  the  consideration.  For 
experimental  parameters  the  following  relation  between  har- 
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monic  amplitudes  was  found, 

^--1.57xl0-^  ^  =  1.025. 

/lo  ^0 

Thus  we  were  able  to  ignore  the  plus  first  harmonic  and  to 
approximate  the  magnetic  field  as 


(^eouRo^^h{poRo)cos{^z/2)  (43) 

This  expression  has  been  used  to  calculate  the  Q-value  of  the 
TT-mode  in  Fig.  6. 
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Theory  of  relativistic  backward  wave  oscillators  operating  near  cutoff 
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A  numerical  model  for  analyzing  backward-wave  oscillators  (BWOs)  operating  near  the  upper 
edge  of  the  transmission  band  is  presented.  The  model  is  used  to  calculate  starting  currents  for 
two  finite  length  devices,  an  X-band  BWO  (/=8.4  GHz)  and  a  J-band  BWO  (/=5.5  GHz). 

The  operating  frequency  and  efficiency  predicted  by  the  nonlinear  numerical  simulations  are 
compared  with  experimental  data  for  each  device. 


I.  INTRODUCTION 

Backward- wave  oscillators  (BWOs)  utilize  a  high- 
current  electron  beam  to  produce  high-power,  coherent  ra¬ 
diation  in  the  centimeter  and  millimeter  wavelength  re¬ 
gime,  Under  certain  voltage  and  beam  current  operating 
conditions,  a  backward-wave  oscillator  can  operate  near 
the  upper  edge  of  the  transmission  band,  where  the  group 
velocity  of  the  electromagnetic  wave  goes  to  zero.  Most 
frequently,  this  corresponds  to  operation  at  the  so-called 
‘‘77’*  point,  where  the  electromagnetic  fields  have  a  period 
twice  that  of  the  structure.  There  are  a  number  of  practical 
reasons  for  operation  near  this  point.  First,  the  start  cur¬ 
rent  usually  has  its  lowest  value  in  this  case.  Thus,  if  the 
length  of  the  interaction  region  and  the  available  beam 
current  are  constrained,  one  can  still  design  a  device  that 
will  operate  efficiently  by  placing  the  operating  point  at  the 
minimum  of  the  group  velocity.  Second,  once  a  particular 
structure  is  chosen,. its  operating  point  can  be  varied  by 
varying  beam  voltage.  The  highest-frequency  operation 
will  occur  at  the  ‘V”  point.  Thus,  it  is  natural  to  consider 
operation  in  this  regime. 

For  frequencies  near  the  “tt”  point  the  cold  structure 
dispersion  relation  can  be  approximated  as  a  quadratic 
function  of  the  wave  number.  A  theoretical  model,  similar 
to  those  presented  in  Refs.  1-3,  is  developed  in  this  paper 
to  describe  the  operation  of  the  device  in  this  regime.  The 
model  consists  of  the  self-consistent  set  of  equations  that 
describe  the  slow  evolution  of  the  envelope  of  the  radiation 
field  and  the  relativistic  motion  of  the  electrons  along  a 
strong  magnetic  field.  Numerical  calculations  of  the  start¬ 
ing  current  are  performed  and  compared  with  an  analytic 
expression  for  the  starting  current  derived  by  assuming  a 
fixed  field  profile. 

To  compare  with  experimental  results,  the  effects  of  a 
time-varying  voltage  and  beam  current  pulse  are  included 
in  the  theoretical  modeling.  We  consider  two  BWOs,  an 
X-band  BWO  (/=8.4  GHz)  and  a  J-band  BWO  (/=5.5 
GHz).  We  compare  predictions  of  peak  efficiencies  and 
operating  frequencies  from  the  theoretical  model  with  ex¬ 
perimental  results."^’^  From  studies  performed  on  the 
X-band  BWO,  agreement  is  shown  between  experiment 
and  theory.  However,  as  will  be  seen,  some  aspects  of  the 
J-band  BWO  experiment  are  not  explained  by  the  model. 

Earlier  theoretical  modeling  of  the  X-band  device^  was 
developed,  with  the  assumption  that  the  interaction  takes 


place  at  frequencies  far  away  from  the  zero  group  velocity 
point.  However,  at  high  voltages  this  assumption  breaks 
down.  The  present  model  allows  us  to  accurately  model  the 
behavior  of  the  device  at  operating  points  near  the  upper 
edge  of  the  transmission  band. 


11.  THEORETICAL  MODEL 

In  this  section  we  develop  the  theoretical  model  to 
describe  the  BWO  interaction  near  the  upper  edge  of  the 
transmission  band,  where  the  group  velocity  vanishes.  In 
our  model,  we  assume  the  interaction  between  the  electron 
beam  and  the  electromagnetic  waves  of  the  cold  slow-wave 
structure  is  weak  in  the  following  sense.  Over  a  length 
equal  to  a  period  of  the  structure,  the  electromagnetic 
fields  have  the  same  temporal  and  spatial  dependence  as 
the  cold-structure  fields.  The  effect  of  the  beam  is  to  cause 
the  envelope  of  the  electromagnetic  fields  to  vary  slowly  in 
axial  distance  and  time.  In  addition,  the  beam  produces  a 
small  space-charge  field  that  is  proportional  to  the  beam 
density. 

Our  derivation  will  begin  by  reviewing  the  approach  of 
Refs.  4  and  6.  The  electromagnetic  fields  in  the  structure 
can  be  written  as  follows: 

E(x,/)  =  [E/x,^)e(z,r)  +c.c.,  ( 1 ) 

B(x.r)  =  [B^(x,A:)e(z,/)  +  B'^>(x,A)]e"*"-''">+c.c.,  (2) 

where  6(2, r)  is  the  slowly  varying  amplitude  of  the  elec¬ 
tromagnetic  waves,  E^(x,/:)  and  B^(x,/c)  are  the  periodic 
eigenfunctions  of  the  empty  slow-wave  structure  and 
E^^^(x,/:)  and  B^^^(x,A:)  are  the  “space-charge”  fields.  The 
space-charge  fields  represent  the  first-order  correction  (in 
beam  density)  to  the  vacuum  fields  and  contain  both  a 
periodic  dependence  due  to  the  corrugations  and  a  slow 
dependence  [analogous  to  that  imparted  by  the  envelope 
e(Zyt)  to  the  vacuum  fields],  reflecting  the  degree  of  bunch¬ 
ing  at  a  particular  point  in  the  interaction  region.  Note  that 
the  space-charge  fields  are  also  modulated  by  an  exponen¬ 
tial  with  wave  number  k  and  frequency  o).  This  is  due  to 
the  fact  the  these  fields  are  the  responses  to  the  bunched 
beam,  which  is  moving  with  speed  v^^r^oj/k,  A  more  de¬ 
tailed  treatment  of  the  space-charge  field  will  be  given  in 
Appendices  A  and  B. 
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(5) 


FIG.  1.  A  schematic  picture  of  a  cold-structure  dispersion  coik)  of  a 
TMoi  mode  of  an  infinitely  long  pcnodic  structure.  The  straight  lines  that 
intersect  the  dispersion  curve  are  the  electron  beam  Doppler  lines, 
(i)  =  kv^  where  ‘S  the  initial  beam  velocity.  The  two  Doppler  beamlines 
correspond  to  beam  voltages  of  550  and  350  kV. 


The  periodic  eigenfunctions  of  the  slow-wave  structure 
satisfy  Maxwell’s  equations  for  a  wave  of  frequency  co  and 
Floquet  wave  number  k,  namely 

E^  =  VXB^i  ikzXBp^  (3a) 


-^B^=VXE^+/^iXE^, 


(3b) 


where  f  is  the  unit  vector  in  the  axial  direction.  The  solu¬ 
tion  to  (3a)  and  (3b)  with  the  boundary  condition  that  the 
tangential  component  of  vanishes  at  the  metal  wall  gives 
the  cold-structure  dispersion  relation,  a){k).  For  example, 
consider  an  empty  corrugated  structure  with  period  d  and 
wall  radius  r(z)  =  ray[l— A  cos(2A:oz)],  where  is  the  av¬ 
erage  wall  radius  and  k^^^rr/d  is  the  wave  number  at  the 
‘V’  point.  The  dispersion  curve  for  the  lowest  axisymmet- 
ric  mode  of  this  structure  is  plotted  in  Fig.  1  for  the  case 
where  r^^—\.5  cm,  A=0,1,  and  t/=1.67  cm.  Also  shown 
are  the  beam  dispersion  curves  for  two  values  of  the  beam’s 
axial  velocity,  corresponding  to  energies  of  350  and  550 
kV.  The  significance  of  these  curves  will  be  discussed  sub¬ 
sequently. 

To  determine  the  equations  governing  the  slowly  vary¬ 
ing  amplitude  of  the  electromagnetic  fields,  one  substitutes 
(1)  and  (2)  into  Maxwell’s  equations. 


1  dE  47r 

-^=VXB  — j, 
c  dt  c 


(4a) 


1  dB 

VXE.  (4b) 

c  dt 

Then  one  subtracts  the  dot  product  of  (4a)  with  B*  from 
the  dot  product  of  (4b)  with  E*  and  integrates  over  the 
volume  contained  in  one  period  of  the  structure.  Using 
(3a)  and  (3b),  the  periodicity  of  E^,  E^^\  B^  and  B^^\  and 
the  boundary  condition  that  the  tangential  component  of 
the  electric  field  vanishes  on  the  metal  wall,  the  following 
slowly  varying  amplitude  equation  is  obtained:'^’^ 


U 


where  the  group  velocity  was  shown  to  be 


4,6 


■X,  z- 


(e,xb;+e*xb,) 

u 


(6) 


the  normalized  energy  stored  in  one  period  of  the  structure 
is  given  by 


C/=T  Jrfz' (|Ep|'+|B,|'),  (7) 

and  J  is  the  component  of  the  current  oscillating  at  (o.  The 
integrals  in  (5),  (6),  and  (7)  are  carried  out  over  the 
volume  contained  in  one  period  of  the  structure.  The  inte¬ 
grands  in  (6)  and  (7)  are  strictly  periodic  in  z,  and  there 
is  no  ambiguity  in  the  evaluation  of  the  integral.  In  the  case 
of  ( 5 )  the  integral  is  carried  out  over  the  range  of  z'  values 
of  length  dy  the  period  of  the  structure,  centered  about  the 
point  z.  The  formal  basis  for  this  type  averaging  is  pre¬ 
sented  in  Appendix  A.  To  be  specific,  we  note  that  in  the 
case  of  a  symmetric  structure  the  volume  integral  may  be 
written  as 


When  operating  far  away  from  the  “tt”  point  (for  ex¬ 
ample,  with  the  lower  beam  voltage  of  Fig.  I ),  the  disper¬ 
sion  relation  of  the  cold  structure  can  be  approximated  by 
a  straight  line  with  slope  which  is  constant.  This  is  the 
situation  that  has  been  treated  previously,  and  is  well  de¬ 
scribed  by  (5).  As  the  operating  point  moves  closer  to  the 
zero  group  velocity  point  as  the  voltage  is  raised  (for  ex¬ 
ample,  with  the  higher  beam  voltage  of  Fig.  1 ),  the  straight 
line  approximation  breaks  down  and  (5)  requires  modifi¬ 
cation.  A  rigorous  derivation  of  the  amplitude  equations 
for  this  case  is  presented  in  Appendix  A.  Here  we  present 
a  heuristic  derivation  that  yields  the  same  result.  Near  the 
upper  edge  of  the  transmission  band  the  dispersion  relation 
is  approximated  by  the  parabola 

I  , 

oy(k) {k  k^)  ,  (8) 


where  a>o  and  k^  are  the  frequency  and  wave  number  at  the 
“tt”  point  (zero  group  velocity),  and  dvg/dk  is  a  negative 
constant.  Thus,  in  this  range  of  wave  numbers  we  have 
Vg(k)  =  {k—kQ)dUg/dk.  We  now  regard  k  as  an  operator, 
which  near  the  “tt”  point,  is  written  as  k  =  kQ—id/dz, 
where  the  derivative  with  respect  to  z  operates  on  the 
slowly  varying  envelope.  The  correct  generalization  of  ( 5 ) 
is  obtained  by  replacing  v^/dz  in  ( 5 )  by  the  following: 


d  -idVgiki^)  d^ 


(9) 


Substituting  (9)  into  (5),  we  obtain  the  following  slowly 
varying  amplitude  equation: 
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de  i  dvg  d^€ 
Tt~2dk^ 


E** 

U  ’ 

(10) 


where 


dk 


_  d  (e,xb;+e*xbp 

dk  U 

(11) 


Equation  ( 10)  is  the  envelope  equation  for  an  electro¬ 
magnetic  wave  driven  by  a  source  term  that  represents  the 
interaction  of  the  cold-structure  fields  with  an  electron 
beam  current.  The  envelope  varies  slowly  in  space  com¬ 
pared  to  the  period  of  the  structure  d  and  the  wavelength 
Ztt/A:,  and  varies  slowly  in  time  compared  to  Itt/cu.  The 
boundary  conditions  for  the  slowly  varying  amplitude  that 
are  to  be  applied  at  both  ends  of  the  structure,  are  of  the 
form 


d€ 

—  (2=0)  ( I2a) 


de 

—  (2=  L)  =  -//c^,e(z=  L),  ( I2b) 

where  and  are  parameters  that  depend  on  the  cou¬ 
pling  of  the  fields  inside  and  outside  the  structure.  By  anal¬ 
ogy  to  the  theory  of  transmission  lines,  these  quantities  can 
be  thought  of  as  admittances  describing  the  fields  outside 
the  structure.  These  boundary  admittances  can  be  written 
in  terms  of  the  reflection  coefficients  at  the  ends  of  the 
periodic  structure.  Consider  the  following  slowly  varying 
amplitude: 

(13) 

where  6^  represent  the  forward  and  backward  propagating 
waves,  and  A/c  and  Ao;  represent  the  shift  in  the  wave 
number  and  frequency  away  from  the  point.  Inserting 
the  above  expression  into  the  boundary  conditions  (12a) 
and  (1 2b),  one  obtains  the  following  expression  for  the 
reflection  coefficient: 


^k-ku 

^~^k  +  k„' 


(14) 


for  /=0,1,  where  Ro=e+/€_  and  cor- 

responding  to  the  reflection  at  2=0  and  2=  L,  respectively. 
As  the  wave  number  approaches  /cq,  ^k  goes  to  zero  and 
the  reflection  /?,  goes  to  —  1.  Now  suppose  that  the  reflec¬ 
tion  coefficient  Rf{k)  has  been  determined  by  some  other 
means.  For  example,  the  reflection  coefficient  for  sinusoi¬ 
dally  rippled  waveguides  connected  to  waveguides  of  con¬ 
stant  wall  radius  has  been  determined  numerically  in  Ref. 
7.  It  isfound  that  Ri(k)  —  —  I  as  k^k^,  the  wave  number 
at  the  “tt”  point.  The  boundary  condition  parameter 
can  then  be  determined  by  solving  ( 14),  assuming  Ri{k)  is 
known  and  taking  the  limit  as  k-*k^. 


kf,i=  Urn  (k-ko) 


\-Ri(k) 

\^Ri(k)‘ 


(15) 


By  doing  this  one  obtains  a  unique  value  for  the  boundary 
admittances  valid  for  operation  near  cutoff. 

From  (15)  ones  sees  that  the  boundary  admittance  is  a 
measure  of  how  the  magnitude  of  reflection  coefficient  for 
a  particular  structure  goes  to  unity  as  the  operating  wave 
number  approaches  k^,  the  wave  number  at  cutoff.  Large 
values  of  the  boundary  admittance  correspond  to  struc¬ 
tures  with  large  values  of  the  magnitude  of  the  reflection 
coefficients  for  ail  operating  wave  numbers  under  consid¬ 
eration,  whereas  small  values  of  the  boundary  admittance 
correspond  to  structures  with  low  reflection  coefficients, 
except  when  operating  near  the  cutoff. 

To  evaluate  the  source  term  in  (10)  we  need  to  char¬ 
acterize  the  electron  beam  current.  Consider  an  annular 
electron  beam  immersed  in  a  strong  magnetic  field.  Assum¬ 
ing  the  beam  has  initially  only  an  axial  component  of  mo¬ 
mentum  and  the  magnetic  field  confines  the  motion  of  the 
electrons  along  the  2  axis,  then  only  the  interaction  of  the 
axial  current  with  the  axial  electric  field  needs  to  be  con¬ 
sidered.  For  synchronism  between  the  electromagnetic 
fields  and  the  electron  beam,  we  pick  out  the  component  of 
the  current  that  varies  as  Therefore,  we  average  the 
current  as  follows: 

J,/(Uof=z^  f  (16) 

ZTT  Jo 

The  current  due  to  an  annular  electron  beam  with  radius 
can  be  written  as  follows: 


Pz 

yiPz)^n 


Ijrr 


dPzy 


(17) 


where  /{z^p^d)  is  the  nonthermal  electron  distribution 
function  defined  as 


/(2,/7,,/)=«0l^^  J:  dtQd[z~z(t,to)\3{p,-p,(UtQ)], 


(18) 


where  Wq  is  the  initial  electron  density  per  unit  length,  ^20  is 
the  initial  axial  electron  velocity  and2(r,/o)  and/7^(r,ro)  are 
the  position  and  momentum,  respectively,  of  a  particle  at 
time  t  that  entered  the  interaction  region  of  the  structure  at 
time^  fo- 

Substituting  the  above  expression  for  the  current  into 
(10)  and  then  normalizing  the  equation,  we  arrive  at  the 
following  slowly  varying  amplitude  equation: 


da  d^a 


(19) 


with  the  following  normalized  boundary  conditions: 


^(^=0)=/7:«,Lfl(^=0),  (20a) 

^  (^=  1)  =  ,La(^=l),  (20b) 

where  xlj=k(^—o)Qt  is  the  electron  phase,  ^=z/ L  is  the 
normalized  axial  distance,  r=^\dPg/dk  L\  ct/L  is  the  nor¬ 
malized  time,  and  0(2, r)  =qLE^pi^ri^)e{z,t)/m<?'  is  the  nor¬ 
malized  amplitude.  Here  we  have  assumed  that  dvg/dk  <  0, 
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appropriate  to  the  upper  edge  of  the  transmission  band. 
The  angled  bracket  denotes  the  average  over  the  initial 
phase  of  the  electrons  with  respect  to  the  high-frequency 
field.  In  addition,  we  have  introduced  a  normalized  current 
defined  as 


(I\(W  Cir,) 

\dpg/dkL\ 


(21) 


where  I \=mc^/q=  1.7X  10'*  A  and  C(rj,)  is  the  coupling 
coefficient.  The  coupling  coefficient  is  a  measure  of  how 
well  the  electromagnetic  fields  of  the  cold-structure  couple 
to  the  electron  beam.  The  coupling  coefficient  is  defined  as 


_ d\  Sjdz  E^p{r,,^)\^ _ 

”  /o/;“‘'’27rr  dr  dz[  \  E,(r^)  | '+  i  Bp(r^)  | ']  ’ 

(22) 


where  the  numerator  is  the  square  of  the  average  of  the 
axial  electric  field  over  one  period  of  the  structure  and  the 
denominator  is  the  total  electromagnetic  field  energy 
stored  in  one  period  of  the  structure.  For  our  model  we  use 
the  value  of  the  coupling  coefficient  at  the  radius  of  the 
electron  beam  for  the  cold  structure  electromagnetic  fields 
with  frequency  coq  and  wave  number  k^. 

An  electron  in  the  beam  is  characterized  by  its  phase 
xl)  =  k^  —  ii)Qt  and  its  energy  mc^y,  where 
y  =  ( 1  Since  the  electron  beam  is  assumed  to 

be  immersed  in  a  strong  axial  magnetic  field,  only  motion 
in  the  axial  direction  is  considered.  As  the  beam  travels 
through  the  interaction  region,  its  phase  and  energy  evolve 
according  to 


dxb 

-  +  ...^=^,-0,0, 


(23) 


dy  dy  u.  q 
dt  3z  m 


(24) 


where  q  is  the  electron  charge.  The  particles  respond 
strongly  only  to  the  synchronous  spatial  harmonic  of  the 
axial  electric  field.  Therefore,  in  (24),  we  average  the  axial 
electnc  field  over  one  period  of  the  structure.  The  quantity 
E^p  is  the  axial  component  of  the  space-charge  field.  In 
Appendix  B  this  field  is  shown  to  be  proportional  to  the 
alternating  beam  current  density, 

—  El^’oc  — (25) 

C  -p  c 


Consistent  with  our  previous  assumption  of  slow  time  ev¬ 
olution,  the  dominant  terms  in  the  above  equations  of  mo¬ 
tion  are  the  derivatives  with  respect  to  the  axial  coordinate. 
Here  we  are  basically  assuming  that  the  slowly  varying 
field  amplitude  changes  little  during  the  time  in  which  an 
electron  is  in  the  structure.  As  an  estimate,  this  requires 
dPg/d{kL)  4p2-  With  these  assumptions  the  equations  of 
motion  reduce  to  the  following: 


dih  coqL 

—  —k  T  — 

dr ""  cp,  ^ 


(26) 


|=2Re[(a(^.r)-/^(.-''^>,,)e^ 

The  space-charge  parameter  is  defined  as 

where  C2  is  the  coefficient  describing  the  reduction  of  the 
plasma  frequency  due  to  the  metal  wall  of  the  periodic 
structure.  See  Appendix  B  for  a  precise  definition  of  C2. 
For  our  calculation  we  use  the  coefficient  C2  calculated 
assuming  an  annular  beam  in  a  smooth  waveguide  of 
radius*’^  r^=r^^.  The  analytic  expression  for  Cj  is  given  by 

1  . Ko{ki_  rjloiki 

'■*V“^o(*i  ^»)/o(^i  Oj’ 
(29) 

where  Iq  and  Kq  are  modified  Bessel  functions  and 
kl  ={kl-(Oo/c^).  This  concludes  the  derivation  of  our 
theoretical  model.  Equations  (19),  (20),  (26),  and  (27) 
describe  mathematically  our  model  and  will  be  solved  nu¬ 
merically  to  determine  the  operating  frequency  and  effi¬ 
ciency  for  two  different  devices.  In  the  next  section  we  use 
this  model  to  calculate  the  starting  current  for  the  two 
devices.  In  Sec.  IV,  results  from  nonlinear  simulations  will 
be  presented. 


(27) 


(28) 


III.  STARTING  CURRENT 

The  start  current  is  defined  as  the  current  above  which 
the  linear  gain  exceeds  losses,  and  signals  are  self-excited. 
Several  authors  have  investigated  conditions  for  self¬ 
oscillation  near  the  edge  of  the  transmission  band.^  Pre¬ 
sented  here  is  both  an  analytic  and  numerical  calculation 
for  the  starting  current  for  a  BWO  operating  near  the  zero 
group  velocity  point.  Starting  current  calculations  are  pre¬ 
sented  in  Ref.  6  for  operation  far  from  the  ‘V’  point.  We 
will  compare  this  starting  current  calculation  based  on  the 
linear  dispersion  curve  valid  far  from  the  ‘V”  point  with 
our  starting  current  calculations  based  on  a  quadratic 
curve  valid  near  the  ‘V”  point. 

For  the  analytic  calculation  of  the  starting  current  we 
assume  a  fixed  field  profile  for  the  slowly  varying  envelope 
and  use  perturbation  theory  to  solve  the  equations  of  mo¬ 
tion  for  the  system.  To  do  this,  let  us  choose  the  normal¬ 
ized  slowly  varying  amplitude  of  the  electromagnetic 
waves  to  have  the  following  form: 

a(^)  =aosin(A/:  L^),  (3(^) 

where  A*  for  /2  =  1,2,...,A^,  where  N  is  the  number 

of  periods  of  the  slow-wave  structure.  Such  a  field  profile 
will  result  when  the  boundary  condition  parameters 
are  very  large.  Using  energy  balance  for  the  amplitude 
equation  (19)  along  with  the  appropriate  boundary  condi¬ 
tions  (20a)  and  (20b),  the  analytic  expression  for  the 
starting  current  is 
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TABLE  I.  Simulation  parameters. 


Parameters 

Backward  wave  oscillator 
X-band  BWO 

J-band  BWO 

/o  (GHz) 

8.7395 

5.8480 

d  (cm) 

1.67 

2.15 

N 

8 

10  14  35 

r^y  (cm) 

1.500 

8.440 

rfy  (cm) 

0.775 

7.375 

C(r,) 

7.24X10-^ 

2.91X10"^ 

00 

1.169 

5.80 

1.169 

d0/dkd 

^0.4613 

-2.6834 

V  (kV) 

300-^ 

200-900 

I  (kA) 

0.10-0.25 

2.0-10.0 

r^(dycplorl  \d0/dkd\ 

4ir(z.j  oJoL  C(r^)F(Ak  L) 

(fcjoL  kjj^L) 

where  the  function  F(A/c  L)  is  defined  as  follows: 


F{Ak  L)-- 


G{k_)-G{k^) 


K  _K  )sin“(K-^  -j-K^) 

where  k  ^  =  [A/:  ±  (  /cq  ^zo ) ]  ^ 


G(/f)=K-  — 


sin^  K 


(32) 


(33) 


In  the  above  calculation  we  have  allowed  for  outgoing 
waves  (finite  /c^  and  /c^i)  at  both  ends  of  the  structure.  If 
at  the  input  of  the  structure  (z<0)  modes  are  cut  off  [a(z 
=  0)— 0]  then,  according  to  (20a)  this  case  is  treated,  by 
allowing  to  become  infinite. 

Equation  (31)  shows  the  starting  current  is  inversely 
proportional  to  the  boundary  admittance.  This  is  reason¬ 
able,  since  by  (15)  a  large  value  of  implies  the  reflection 
coefficient  is  close  to  unity.  As  the  boundary  admittance 
increases,  the  starting  current  of  the  device  decreases.  The 
boundary  admittance  at  the  “tt”  point  is  inversely  propor¬ 
tional  to  the  slope  of  the  reflection  coefficient  therefore  as 
the  slope  approaches  zero  [corresponding  to  a  constant 
reflection  /?(/:)  =  !  for  all  k]  more  energy  is  stored  in  the 
device.  Less  energy  radiates  out  of  the  device,  making  it 
easier  to  start  electromagnetic  waves  oscillating. 

Using  our  model  we  calculate  numerically  the  starting 
current  for  two  devices  and  compare  our  analytic  and  nu¬ 
merical  results.  Table  I  shows  the  operating  parameters  for 
the  two  devices.  The  values  of  quantities  such  as  frequency, 
coupling  coefficient,  boundary  condition  parameters,  and 
group  velocities  are  obtained  using  the  computer  code  de¬ 
scribed  in  Ref.  7,  The  first  device  considered  is  an  X-band 
BWO.  This  device  has  been  extensively  studied  in  Refs.  4 
and  6,  using  the  model  based  on  a  linear  dispersion  curve 
(constant  Ug).  However,  the  assumption  that  the  group 
velocity  is  constant  breaks  down  for  voltages  above  400 
kV.  As  the  voltage  is  increased,  the  beamline  moves  toward 
the  ‘V’  point,  where  the  dispersion  curve  rolls  over. 
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FIG.  2.  Starting  current  versus  voltage  for  the  X-band  BWO  with  period 
length  d=  1.67  cm,  number  of  periods  A  =  8,  cutoff  frequency  /o= 8.7395 
GHz,  average  wail  radius  1.50  cm,  and  beam  radius  /'^=0.775  cm. 


Therefore,  for  small  changes  in  voltage  the  group  velocity 
of  the  electromagnetic  waves  is  no  longer  constant.  Figure 
2  shows  the  starting  current  calculation  based  on  the  zero 
group  velocity  model  (both  analytic  and  numerical)  and 
the  start  current  from  the  model  of  Refs.  4  and  6.  Clearly, 
the  analytic  and  numerical  calculations  based  on  the  qua¬ 
dratic  dispersion  agree  very  well.  As  expected,  the  slowly 
varying  amplitude  is  very  nearly  equal  to  the  assumed  sine 
profile.  The  cusp  in  the  starting  current  corresponds  to  a 
change  in  the  axial  mode  number.  For  voltages  between 
200  and  400  kV  the  field  profile  has  A/:  £  =  37r,  and  for 
voltages  greater  than  400  kV  the  field  profile  has  A/c  L 
=  2v".  For  the  range  of  voltage  from  300  to  4(X)  kV,  the  two 
models  ( linear  and  quadratic  dispersion )  agree  well,  how¬ 
ever,  for  voltages  greater  than  400  kV  our  model  near  cut¬ 
off  (quadratic  dispersion)  is  more  reliable. 

The  second  device  considered  here  is  a  J-band  BWO. 
As  detailed  in  Table  I,  this  device  operates  in  the  range  of 
voltage,  where  the  dispersion  relation  for  the  cold- 
structure  fields  is  quadratic.  The  group  velocity  of  the  elec¬ 
tromagnetic  wave  is  very  nearly  equal  to  zero.  The  previ¬ 
ous  linear  dispersion  model  is  not  capable  of  treating  this 
device.  Therefore  for  this  device  our  “rr”  point  model  is 
preferable.  For  this  device  three  lengths,  iV=14,  16,  35, 
periods  are  studied.  Figure  3  shows  the  analytic  and  nu¬ 
merical  starting  currents  for  the  three  structures.  A  dis¬ 
crepancy  exists  between  the  two  curves  because  the  actual 
slowly  varying  field  profile  obtained  from  the  numerical 
code  differs  from  the  sine  wave  assumed  by  the  analytic 
calculation.  This  discrepancy  is  shown  in  the  inset  figure  of 
Fig.  3  for  the  16  period  structure.  We  note  that  the  bound¬ 
ary  admittances  for  this  device  are  not  nearly  as  large  as  for 
the  X-band  BWO.  Thus,  the  assumption  of  a  sinusoidally 
varying  field  profile  in  the  analytic  calculation  is  not  ex¬ 
pected  to  be  as  good. 


IV.  NONLINEAR  SIMULATIONS 

In  this  section  we  present  the  modifications  to  our  the¬ 
oretical  model  in  order  to  simulate  actual  BWO  experi- 
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FIG.  3.  Starting  current  versus  voltage  for  the  J-band  BWO  with  penod 
length  d=l.\5  cm,  number  of  penods  A^=14,  16,  35,  cutoff  frequency 
/q  =  5.8488  GHz,  beam  radius  1.115  cm,  and  average  wall  radius 
/-^,  =  8.44  cm.  The  profile  in  the  upper  right  corner  is  the  normalized 
slowly  varying  amplitude  versus  normalized  length  for  voltage  V 
=  750  kV. 


ments.  Comparisons  between  the  nonlinear  simulations 
and  experimental  results  for  operating  frequency  and  peak 
efficiency  are  shown. 

To  simulate  the  BWO  experiments,  we  include  in  our 
theoretical  model  a  time-dependent  voltage  and  beam  cur¬ 
rent  pulse.  We  assume  a  rising  exponential  for  the  beam 
current  and  voltage,  namely 

/(/)  =  /p(l-e''^V/'0P  (34) 


K(0= 


(35) 


where  is  the  rise  time  and  tp,  Vp,  and  Ip  are  parameters 
chosen  to  best  model  the  characteristics  of  the  beam.  We 
assume  the  Child-Langmuir  law  is  approximately  valid, 
and  let  p=  1.5.  To  model  the  end  of  the  pulse,  we  assume 
the  voltage  and  current  drop  exponentially  at  some  time 
t=tf.  In  addition,  we  include  the  effects  of  direct  current 
space  charge.  When  the  electron  beam  enters  the  electro¬ 
magnetic  structure  the  beam  slows  down  due  to  a  buildup 
of  an  electrostatic  potential.  This  voltage  depression  re¬ 
duces  the  electron  beam  energy.  The  reduction  in  the  beam 
energy  is  calculated  using  the  following  formula:*^ 


/ 

A7=2  “  In 

f  A 


(36) 


where  r^v  is  the  average  wall  radius  of  the  corrugated  struc¬ 
ture  and  is  the  radius  of  the  electron  beam.  Since  the 
voltage  and  current  are  time  dependent,  the  initial  value  of 
the  relativistic  factor  y,  is  given  by  7^=70— ^7*  where 

The  two  BWO  devices  considered  here  do  not  operate 
exactly  at  the  “tt'’  point,  therefore,  we  made  two  additional 
modifications  to  our  theoretical  model  to  more  accurately 
simulate  these  devices.  The  first  modification  to  our  model 
was  to  change  the  equations  of  motion  to  included  the 
frequency  shift  away  from  the  “tt”  point.  Equation  (26) 
then  becomes 


FIG.  4.  (a)  Operating  frequency  versus  voltage  for  the  X-band  BWO 
with  period  length  4/=  1.67  cm,  number  of  penods  A^  =  8,  cutoff  frequency 
7^^8  7395  GHz,  average  wail  radius  r,,  =  1.50  cm,  and  beam  radius 
0.775  cm.  (b)  Peak  efficiency  versus  voltage. 


dtp  ^  ^  (cL)o+Aa))L 

where  Aco  is  the  average  frequency  shift  obtained  from 
(19).  The  slowly  varying  envelope  is  a  complex  quantity 
representing  the  amplitude  and  phase  of  the  wave.  The 
frequency  shift  Aco  is  the  rate  of  change  of  the  phase  of  the 
amplitude,  a.  To  find  Aco  we  multiply  ( 19)  by  a*,  integrate 
over  the  length,  and  take  the  imaginary  part  to  find  the 
average  frequency  shift.  In  addition,  since  the  coupling 
coefficient  varies  very  sharply  with  wave  number  near  the 
“tt”  point,  we  included  the  coupling  coefficient  as  a  func¬ 
tion  of  operating  wavelength.  Therefore,  after  each  correc¬ 
tion  of  the  operating  frequency  we  updated  the  value  of  the 
coupling  coefficient. 

Using  our  modified  code  for  a  finite-duration  current 
and  voltage  pulse,  we  compare  numerical  simulation  with 
experimental  results.^’^  First,  consider  the  X-band  BWO 
with  cutoff  frequency  /o=8.7395  GHz,  period  ^=1.67 
cm,  and  length  L  =  13.36  cm.  Figure  4(a)  shows  the  com¬ 
parison  of  the  operating  frequency  versus  voltage  and  Fig. 
4(b)  shows  the  comparison  of  the  peak  efficiency  versus 
voltage.  With  the  above  improvements  to  our  model  we 
obtain  very  good  agreement  between  the  calculated  and 
measured  frequencies,  as  shown  in  Fig.  4(a).  The  small 
discrepancy  between  the  measured  and  calculated  frequen- 
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FIG.  5.  Axial  modes  for  the  X-band  BWO  with  period  length  d—\.tl 
cm,  number  of  periods  A  =  8,  cutoff  frequency  /o  =  8.7395  GHz,  average 
wall  radius  /*„,=  1.50  cm,  and  beam  radius  /■ft  =  0.775  cm.  The  solid  line  is 
the  cold-structure  dispersion  curve,  where  the  open  circles  represent  dif¬ 
ferent  axial  modes.  The  dashed  line  is  the  parabolic  dispersion  curve  used 
in  the  theoretical  model  for  operation  near  the  "tt"  point,  where  the 
closed  circles  represent  different  axial  modes.  The  frequencies  for  the  axial 
modes  of  the  cold-structure  dispersion  are  listed  on  the  right  and  the 
frequencies  for  the  axial  modes  of  the  parabolic  dispersion  arc  listed  on 
the  left. 


cies  is  explained  by  the  difference  between  the  actual  dis¬ 
persion  relation  and  its  parabolic  approximation.  Figure  5 
shows  the  cold  structure  dispersion  relation  with  the  axial 
modes  for  the  structure,  assuming  a  shorted  boundary  con¬ 
dition  [/2(/c)  =  l]  at  the  end.  The  dotted  line  shows  the 
assumed  parabolic  dispersion  relation  for  the  structure 
when  operating  near  the  ‘V”  point  with  the  corresponding 
axial  modes.  Fof  the  first  three  axial  modes  the  quadratic 
dispersion  agrees  with  the  cold-structure  dispersion  to 
within  65  MHz. 

The  simulated  frequency  for  diode  voltages  below  400 
kV  is  not  shown  in  Fig.  4(a),  since  in  the  simulations  the 
device  did  not  start  in  this  range  of  voltages.  It  should  be 
noted,  however,  that  if  the  diode  current  in  the  simulations 
is  increased  the  device  will  operate  with  a  diode  voltage  of 
365  kV  in  the  third  axial  mode  (/=8.19  GHz,  which  is 
not  the  observed  frequency).  Additionally,  if  the  duration 
of  the  pulse  is  sufficiently  increased  (nearly  doubled)  the 
device  will  operate  in  the  second  axial  mode  (/=8.48 
GHz)  at  365  kV,  however,  the  efficiency  is  less  than  \%. 
With  the  diode  voltage  set  to  305  kV,  neither  raising  the 
current  nor  lengthening  the  pulse  leads  to  efficient  opera¬ 
tion  in  the  second  axial  mode  {/=8.19  GHz).  This  is  a 
consequence  of  the  fact  that  for  the  measured  frequency 
(8.5  GHz)  and  the  wave  number  inferred  from  Fig.  5 
(k^ir/d-^ln/ L),  the  resonant  energy  is  211  kV.  Since 
efficient  operation  requires  the  energy  of  the  entering  beam 
exceed  the  resonant  energy  by  some  margin,  high  efficiency 
operation  at  this  frequency  K=305  kV  is  not  possible. 
Thus,  while  the  calculated  efficiency  shows  the  same  basic 
trend  as  the  experimentally  measured  efficiency,  the  calcu¬ 
lated  efficiency  is  shifted  in  voltage  by  approximately  100 
kV.  The  origin  of  this  shift  is  under  current  investigation. 

We  consider  now  the  J-band  BWO  with  cutoff  fre¬ 
quency  /=5.848  GHz  and  period  c/=2.15  cm.  For  this 
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FIG.  6.  (a)  Operating  frequency  versus  voltage  for  the  J-band  BWO  with 
period  length  d=2. 15  cm,  number  of  penods  A  =10,  cutoff  frequency 
5.8488  GHz,  beam  radius  r.,=7.375  cm,  and  average  wall  radius 
-8.44  cm.  The  solid  line  is  the  best-fit  curve  through  the  simulation 
points,  (b)  Peak  theoretical  efficiency  versus  voltage. 


device  we  used  the  experimentally  measured  voltage  and 
current  versus  time  traces  instead  of  the  exponential  pulse 
described  by  (34)  and  (35).  The  operating  parameters  for 
this  device  are  presented  in  Table  I.  Figures  6(a),  7(a), 
and  8(a)  show  the  comparison  between  the  measured  and 
calculated  frequencies  for  each  of  the  three  structure 
lengths:  ^=10,  14,  and  35,  respectively. 

The  simulated  and  measured  frequencies  of  operation 
show  the  same  general  dependence  on  beam  voltage.  Fur¬ 
ther,  both  the  simulation  and  experiment  show  roughly  the 
same  frequency  shift  from  the  '‘tt”  point  (/= 5.848  GHz). 
The  agreement  appears  to  be  best  for  the  35  period  struc¬ 
ture.  While  the  shorter  structures  exhibit  less  variation  of 
frequency  with  voltage  in  the  experiment  than  the  simula¬ 
tion.  This  suggests  that  the  effective  cavity  Q  in  the  exper¬ 
iment  is  higher  than  in  the  simulation. 

The  peak  theoretical  efficiency  versus  voltage  for  each 
of  the  structures  is  shown  in  Figs.  6(b),  7(b),  and  8(b), 
respectively.  Results  from  the  experimental  study^  of  this 
device  predict  peak  efficiencies  near  5%,  based  on  output 
power  measurements.  Numerically  calculated  efficiencies 
are  between  10%-15%. 

Two  possible  explanations  for  the  relatively  large  dis¬ 
parity  in  predicted  and  measured  efficiencies  are  the  fol¬ 
lowing.  First,  there  is  a  large  uncertainty  in  the  measure- 
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FIG.  7.  (a)  Operating  frequency  versus  voltage  for  the  J-band  BWO  with 
penod  length  d  =  2.I5  cm,  number  of  periods  iV=14,  cutoff  frequency 
/o=  5.8488  GHz,  beam  radius  375  cm,  and  average  wail  radius 

/'^,^8.44  cm.  The  solid  line  is  ilie  heslTit  curve  through  the  simulation 
points,  (b)  Peak  theoretical  efficiency  versus  voltage. 


FIG.  8.  (a)  Operating  frequency  versus  voltage  for  the  J-band  BWO  with 
period  length  d=2. 15  cm,  number  of  periods  A^—35,  cutoff  frequency 
/o=  5.8488  GHz,  beam  radius  rft= 7.375  cm.  and  average  wail  radius 
r,^=8.44  cm.  The  solid  line  is  the  best-tit  curve  through  the  simulation 
points,  (b)  Peak  theoretical  efficiency  versus  voltage. 


ment  of  the  experimental  efficiencies,  as  will  be  described  in 
a  forthcoming  paper. Second,  the  theory  treats  only  the 
lowest-order  symmetric,  transverse  mode  of  the  structure. 
While  this  mode  is  the  one  that  interacts  most  strongly 
with  the  electron  beam,  the  next  two  symmetric  modes 
maybe  important  as  well  since  they  are  above  cutoff  at  the 
operating  frequency.  These  modes  will  be  excited  at  the 
discontinuity  at  the  ends  of  the  structure  and  can  affect  the 
transverse  field  profile  of  the  radiated  power.  A  more  com¬ 
plete  theory  would  follow  the  evolution  of  these  modes  as 
well  as  the  lowest-order  mode  and  include  the  coupling 
among  the  modes  at  the  ends  of  the  device. 

V.  SUMMARY 

We  have  used  the  theoretical  model  for  the  operation 
of  a  BWO  near  the  edge  of  the  transmission  band  to  cal¬ 
culate  starting  current  for  two  structures.  This  model 
agrees  well  with  the  previous  model  fvirmulated  with  non¬ 
zero  group  velocity  assumptions  for  operating  voltages, 
where  the  cold-structure  dispersion  makes  a  transition 
from  a  linear  to  a  quadratic  function  of  wavelength.  After 
modifying  the  model  near  cutoff  to  include  a  time- 
dependent  current  and  voltage,  we  have  simulated  two 
BWO  experiments,  and  compared  the  predicted  and  mea¬ 
sured  results  for  peak  efficiency  and  operating  frequency. 


Predicted  and  measured  frequencies  generally  agree  to 
within  4%,  and  show  the  same  dependence  on  beam  volt¬ 
age.  However,  predicted  and  measured  efficiencies  show 
various  forms  of  disagreement.  In  the  case  of  the  X-band 
BWO,  the  magnitude  of  the  measured  and  predicted  effi¬ 
ciencies  are  in  agreement.  However,  the  peak  efficiencies  in 
the  two  cases  occur  at  different  voltages.  Various  theoret¬ 
ical  attempts  to  induce  the  “model”  BWO  to  operate  at 
lower  voltages  failed.  Thus,  the  efficient  experimental  op¬ 
eration  at  the  reported  voltages  is  unexplained.  In  the  case 
of  the  J-band  BWO,  the  predicted  efficiency  is  substantially 
higher  than  that  which  is  measured.  However,  there  are 
large  uncertainties  in  this  measurement.^  Further,  the 
model  treats  only  the  lowest-order  symmetric  mode  of  the 
structure.  An  improved  theory  that  is  currently  under  de¬ 
velopment  will  employ  a  “full- wave”  solution  for  the  sym¬ 
metric  electromagnetic  fields,  including  coupling  at  the 
output  and  input  ends  of  the  cavity. 
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J^PeiDfX  A:  DERIVATION  OF  SLOWLY  VARYING 
AMPLITUDE  EQUATION 

In  this  appendix  we  will  present  a  rigorous  derivation 
of  the  slowly  varying  amplitude  equation  used  in  our  nu¬ 
merical  simulations.  We  begin  by  assuming  that  the  oper¬ 
ating  point  is  near  the  “tt”  point  and  write  the  electric  and 
magnetic  fields  in  the  structure  as 

E(x,f)  =  E(x,r)e'<*°*-“""+c.c.,  (Al) 

B(x,f)=B(x,r)e'<*o*-‘^'+c.c.,  (A2) 

where  kQ=ir/d  and  coq  is  the  “tt”  point  frequency.  These 
expressions  are  then  substituted  in  Maxwell’s  equations 
(4a)  and  (4b),  yielding  evolution  equations  for  the 
complex  amplitudes  E  and  B, 

1  ICOq  ^  .  477  ^ 

E  =  VXB+/VXB— 

c  dt  c  c 

(A3) 

1  ^  ^ 

— B=VXE  +  /A:ofXE.  (A4) 

c  ot  c 

We  now  introduce  a  multiple  length  and  time  scale 
analysis  to  describe  the  slow  evolution  in  time  and  space  of 
the  amplitude  of  the  field.  Primarily  this  involves  separat¬ 
ing  the  axial  dependence  of  the  fields  into  a  fast  dependence 
on  the  variable  Zo  and  a  slow  dependence  on  the  variable  . 
Here,  variations  with  Zq  will  be  periodic,  with  a  period 
equal  to  that  of  the  structure  describing  the  cold  cavity 
fields,  and  variations  with  will  describe  the  profile  of  the 
fields  over  a  distance  equal  to  the  length  of  the  structure. 
Thus,  our  formal  expansion  is  based  on  the  small  param¬ 
eter  6  =  d/ L,  where  d  and  L  are  the  period  and  length  of 


where  e{zxd)  is  the  amplitude  envelope,  which  is  to  be 
determined  in  higher  order. 

At  this  point  it  is  necessary  to  decide  the  ordering  with 
respect  to  the  parameter  8  of  all  the  terms,  which  we  have 
dropped  from  the  lowest-order  version  of  (A3)  and  (A4). 
For  example,  the  theory  described  in  Refs.  4  and  6,  which 
is  valid  far  away  from  the  ‘V”  point  assumes  that  the  slow 
time  derivative  and  the  beam  current  are  first  order  in  6,  as 
is  the  axial  derivative  with  respect  to  the  variable  Zj .  This 
results  in  (5),  which,  as  we  have  seen,  cannot  be  used  at 
the  “tt”  point  due  to  the  vanishing  of  the  group  velocity. 
The  resolution  is  that  one  must  consider  the  time  deriva¬ 
tive  and  the  beam  current  to  be  second  order  in  8.  Physi¬ 
cally,  this  is  a  consequence  of  the  fact  that  the  dispersion 
curve  is  parabolic  near  the  “rr”  point,  and  the  frequency 
depends  quadratically  on  the  deviations  of  the  wave  num¬ 
ber  from  the  “tt”  point  value.  From  a  formal  point  of  view, 
we  will  now  have  to  go  to  second  order  in  the  expansion 
parameter  to  obtain  the  envelope  equation. 

Proceeding  now  with  our  expansion,  the  first-order 
version  of  (A3)  and  (A4)  are  written  as 

— —  E|  =  TqX B I -|- //cqzX B I -f-zX B^(x,/co)  , 

’  (A12) 

io}c\  ^  ^  ^  dc 

—  B,-VoXEj-f/VXE,-hfXE.(x,^o)  T'-  (A13) 
c 

Equations  (A  12)  and  (A  13)  represent  a  pair  of  coupled, 
inhomogeneous  partial  differential  equations  for  the  first- 
order  fields  Ej  and  Bj.  It  can  be  verified  by  differentiation 
with  respect  to  A:  of  (3a)  and  (3b)  that  the  particular 
solution  of  (A12)  and  (A13)  is 


the  structure,  respectively.  We  then  expand  the  field  and 


axial  derivatives  in  powers  of  6, 

E  =  E^  +  Ei-i-E2+’*‘  , 

(A5) 

B  =  Bf)  -h  B  j  +  B2  +  ■  *  ■  , 

(A6) 

a  a  a 

dz  dz^~^ 

(A7) 

The  beam  current  and  the  time  derivative  of  the  fields 
in  (A3)  and  (A4)  will  be  small  quantities;  exactly  how 
small  will  be  discussed  in  a  moment.  Thus,  the  lowest- 

(9E. 

de 

E,=  -- 

— 

(A14) 

dk 

/  dzi  ’ 

dB, 

d€ 

B,  =  -, 

i  ^ 

(A15) 

dk 

k 

The  homogeneous  solution  has  the  same  form  as  the 
lowest-order  solution  of  (A  10)  and  (All)  and  thus  can  be 
absorbed  into  it. 

Proceeding  to  second  order,  we  find 


order  equations  resulting  from  (A3)  and  (A4)  are 


icon  ^  ^ 

Eo=VoXBo+/MXBo, 

C 

(A8) 

-^Bo=VoXEo+/MXEo, 

(A9) 

which  are  the  empty  structure  equations  (3a)  and  (3b)  for 
(o=o}q  and  which  is  at  the  “tt”  point.  The  solution 

of  these  equations  is 


I  d€  io)f\  A  ^ 

*0 

477 

_ — (A16) 


1  d€  icOfi  ^  u€ 

(A17) 


Eo=E^(x,/:o)e(zi,0. 

Bo=B^(x,A:o)«Ui.'). 


(AlO) 

(All) 


where  we  have  made  use  of  (AlO),  (All),  (A14),  and 
(A15).  Taking  the  dot  product  of  (A16)  with  E*  and 
(A  16)  with  B*  and  then  subtracting,  we  obtain 
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<3B„ 


dk 


XB* 


*0 


a^e 


=  V ’  (B2XE* -E.XBp  -— 

(A18) 

In  obtaining  (A  18)  we  have  made  use  of  the  complex 
conjugate  of  (3a)  and  {3b).  Integrating  (A18)  over  the 
volume  in  one  period  of  the  structure  and  dividing  through 
by  the  electromagnetic  energy  per  |6|^  in  one  period,  we 
obtain  the  slowly  varying  amplitude  equation. 


/  dVg  d^€  r  f  , 


Se  i  dv^  d^€ 
Jt 


u 


(A19) 


where 


U  =  —  ^  dzoj  dh,(.\Ep\^+\B^\^)  (A20) 

is  the  electromagnetic  energy  per  |e|^  in  one  period  and 


dv,  c  r  j-  ,  d  (/t.xb;  -  k:xb,j 

s=rj  -ot - V' — 

(A21) 


is  the  derivative  of  the  group  velocity  with  respect  to  k  at 
the  “tt”  point  for  the  vacuum  electromagnetic  wave.  In 
deriving  (A  19),  we  have  used  the  periodicity  of  E^,  £3, 
and  and  the  boundary  condition  that  the  tangential 
component  of  the  electric  field  vanishes  on  the  metal  wall. 
Equation  (A  19)  is  our  desired  result  and  is  clearly  identi¬ 
cal  to  ( 10),  which  was  derived  heuristically  in  the  text. 


APPENDIX  B:  DEVELOPMENT  OF  SPACE-CHARGE 
FIELDS 

In  Appendix  A  we  introduced  an  ordering  scheme  and 
perturbation  expansion  for  calculation  of  the  fields  in  a 
corrugated  structure  in  the  presence  of  an  electron  beam. 
The  fields  are  represented  as  a  sum  of  terms,  (A5)  and 
(A6),  which  constitutes  a  series  expansion  in  the  small 
parameter  S  =  d/L,  where  d  and  L  are  the  period  and 
len  of  the  structure,  respectively.  As  discussed  in  /  .p- 
pe  A  A,  in  solving  these  equations  we  assumed  that  the 
beai '  current  and  density  was  second  order  in  the  I^aram- 
eter  Thus,  the  lowest-order  fields  (Eo,Bo)  and  (E|,Bi) 
[givr  )y  Eqs.  (Al(  All),  (A14),  and  (A15)]  repre¬ 
sent,  essentially,  vacuun'i  fields  with  slowly  varying  enve¬ 
lopes.  The  second-order  fields  (£2,82)  will  contain  contri¬ 
butions  proportional  directly  to  the  local  charge  and 
current  density.  It  is  these  terms  that  constitute  the  space- 
charge  fields  included  in  our  model. 

To  obtain  these  space  charge  fields  we  must  solve  Eqs. 
(A  16)  and  (A  17).  To  do  this,  we  exploit  the  variation  of 
parameters  technique  used  to  solve  the  first-order  system  in 


Eqs.  (A  12)  and  (A  13).  We  write  the  second-order  fields 
as  the  sum  of  a  space-charge  field  and  a  second-order  cor¬ 
rection  to  the  vacuum  field, 

1  d^E„\ 


2^^ 


(Bl) 


(B2) 


and  insert  these  expressions  in  Eqs.  ( A16)  and  ( A17).  We 
then  use  Eq.  (A  19)  to  eliminate  the  term  proportional  to 
d€/dt  in  Eqs.  (A  16)  and  (A  17).  We  note  that  terms  pro¬ 
portional  to  d^€/dz\  cancel,  since  they  multiply  the  second 
derivative  of  Eqs.  (3a)  and  (3b).  That  is,  by  differentiating 
Eqs.  (3a)  and  (3b)  twice  with  respect  to  k,  we  have  the 
identitv 


id^w 

cdk^ 


E 


/a»o  d~E^ 


c  dk- 


d^B, 

=  VX  ^ 


dk- 


(3-B,  <3Bj 
+  iVX^  +  2;iX  — I 


(B3) 


/  d^o) 
cW- 


co 


+  c  dk- 


d% 
:VX  ^ 


dk- 


d% 


dE, 


(B4) 


After  noting  the  cancellations,  the  space-charge  fields  are 
found  to  satisfy 
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c  ~  ^  ~  c 


+  d\x 
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"<■'>  --VnXEV'*  -t-/A:ofXEi’' 


+  — BV 
c 


d"xi 


_ 


(B6) 

Thus,  the  space-charge  fields  Ej^^Zq^Zi)  and 
are  linearly  proportional  to  the  beam  current. 
These  space-charge  fields  are  referred  to  in  the  main  text  as 
E),^^  and  B^^^  For  the  model  considered  in  our  paper  in 
which  the  beam  is  annular,  and  the  magnetic  field  restricts 
the  current  to  point  in  the  axial  direction  we  have,  follow¬ 
ing  the  steps  leading  to  Eq.  (19), 


Irrr^ 

Note  that  the  ‘‘bunching  function”  (exp[— )]) 
varies  on  the  slow  scale,  Z\,  that  is,  the  bunching  varies 
slowly  over  a  distance  equal  to  one  period  of  the  structure. 
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vary  by  a  large  amount  over  the  length  of  the 
structure.  Thus,  we  may  write  the  space  charge  fields  as  a 
product  of  the  beam  current,  the  bunching  factor,  and  a 
field  profile,  Qtix,  which  depends  on  the  geometry 

of  the  structure,  the  operating  point,  and  the  location  ot 

the  beam, 

(B8) 

The  axial  component  of  this  field  when  evaluated  at  the 
beam  radius  and  averaged  over  axial  distance  enters  the 

equations  of  motion  in  Eq.  (24). 

It  is  reasonable  to  ask  the  basis  for  retaining  the  space- 
charge  field  in  the  interaction,  but  not  retaining  the  higher- 
order  corrections  to  the  vacuum  field.  There  are  two  justi¬ 
fications.  The  first  is  that  the  space-charge  field  produces  a 
qualitatively  different  interaction  with  the  beam  than  the 
vacuum  electromagnetic  field.  This  is  due  to  the  fact  that  it 
is  locally  proportional  to  the  beam  current,  unlike  the  vac¬ 
uum  field,  whose  effect  is  integrated  over  the  length  of  the 
structure.  Thus,  including  higher-order  corrections  to  the 
electromagnetic  field  will  not  change  the  results  as  much  as 
including  the  space-charge  field.  The  second  reason  is  that 
the  space-charge  field,  though  formally  small  in  the  small 
parameter  S  =  d/L,  tends  to  be  peaked  at  the  location  of 
the  beam,  whereas  the  synchronous  component  of  the  elec¬ 


tromagnetic  field  tends  to  be  larger  at  the  wall,  and  de¬ 
pends  also  on  the  ripple  size.  In  a  device  with  smaU  nppl« 
the  space-charge  field  can  dominate.  Thus,  we  retain  only 
the  space-charge  field  and  the  lowest-order  electromagnetic 
field. 
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In  microwave  sources  of  coherent  Cerenkov  radiation  the  electrons  usually  propagate  near  the 
rippled  wall  of  a  slow-wave  structure.  These  ripples  cause  the  periodic  modulation  of  elecuon 
potential  depression,  and  therefore,  lead  to  periodic  modulation  of  electron  axial  velocities.  Since 
the  period  of  this  electrostatic  pumping  is  the  period  of  the  slow-wave  structure  the  parametric 
coupling  of  electrons  to  originally  nonsynchronous  spatial  harmonics  of  the  microwave  field  may 
occur.  This  effect  can  be  especially  important  for  backward-wave  oscillators  (BWO’s)  dnven  by 
high  current,  relativistic  electron  beams.  In  the  paper  both  linear  and  nonlinear  theories  of  the 
relativistic  resonant  BWO  with  periodic  modulation  of  electron  axial  velocities  are  developed 
and  results  illustrating  the  evolution  of  the  linear  gain  function  and  the  efficiency  of  operation  in 
the  large-signal  regime  are  presented. 


I.  INTRODUCTION 

Operation  of  Cerenkov  sources  of  electromagnetic  ra¬ 
diation  is  based  on  the  synchronism  between  the  electron 
axial  velocity,  v^,  and  the  phase  velocity  of  the  electromag¬ 
netic  wave,  Vp^.  To  provide  this  synchronism  slow-wave 
structures  are  used.  Usually  these  structures  are  made  in 
the  form  of  metallic  waveguides  with  periodically  corru¬ 
gated  walls  (see  Fig.  1).  According  to  the  Floquet  theo¬ 
rem,  the  electromagnetic  fields  of  any  mode  in  such  a  struc¬ 
ture  can  be  represented  as  a  superposition  of  spatial 
harmonics 

E=  ^  [fcj— j(2ir/</) )i} ^ Q  (  1  ) 

s 

(A  similar  expression  can  be  written  for  H.)  In  this  equa¬ 
tion  d  is  the  structure  period,  co  is  the  field  frequency,  is 
the  axial  wave  number  of  the  zero  spatial  harmonic,  and  E, 
is  the  harmonic  amplitude.  The  above  mentioned  condition 
of  synchronism 

v,ssVp^=a)/k^,  (2) 

where  k^=ko-s2ir/d,  may  be  fulfilled  either  for  the  zero 
spatial  harmonic  (in  traveling  wave  tubes,  TWT’s)  or  for 
the  minus  first  spatial  harmonic  (in  backward  wave  oscil¬ 
lators,  BWO’s).  Usually  the  theory  of  these  devices  con¬ 
siders  only  the  interaction  of  electrons  with  the  synchro¬ 
nous  harmonic  of  the  electromagnetic  field. 

During  the  last  20  years  high-power  microwave 
sources  based  on  Cerenkov  radiation  of  relativistic  electron 
beams  have  been  actively  studied  both  theoretically  and 
experimentally  (an  extensive  review  of  these  studies  has 
recently  been  done  by  Benford  and  Swegle')-  In  spite  of  a 
large  number  of  publications  there  is  still  at  least  one  effect 
which  may  be  important  for  operation  of  these  devices  but 
which,  as  far  as  we  know,  has  not  been  analyzed  before. 


•’Pennanent  address:  Department  of  Physics.  Moscow  State  University, 
Moscow  119899,  Russia- 


This  effect  is  based  on  the  periodic  axial  modulation  of  the 
velocities  of  electrons  moving  linearly  along  the  metallic 
structure  with  rippled  walls  (see  Fig.  1). 

In  principle,  depression  of  electron  potential  with  re¬ 
spect  to  the  applied  voltage,  which  is  caused  by  the  finite 
distance  between  the  electron  beam  and  the  accelerating 
structure  wall,  is  well  known  and  its  effect  on  the  operation 
of  relativistic  Cerenkov  devices  has  been  discussed 
elsewhere.^  ’  However,  usually  this  effect  is  treated  only  as 
a  reduction  of  the  mean  value  of  electron  energy. 

In  fact,  for  electron  beams  propagating  in  a  metallic 
structure  with  rippled  walls,  the  voltage  depression  de¬ 
pends  also  on  the  electron  axial  position  which  leads  to 
axial  modulation  of  electron  velocities.  The  period  of  this 
modulation  is,  obviously,  the  period  of  the  slow-wave 
structure  d.  Therefore,  this  modulation  may  cause  the  res¬ 
onance  parametric  coupling  between  electrons  and  other 
originally  nonsynchronous  spatial  harmonics.  Of  course,  as 
the  current  grows,  the  voltage  depression  effect  becomes 
more  pronounced.  There  is  also  another  not  so  obvious 
reason  which  makes  this  effect  potentially  important  for 
relativistic  BWO’s.  From  the  scaling  laws  for  Cerenkov 
microwave  sources  with  ultra-relativistic  electron  beams  it 
follows  that  with  the  increase  in  operating  voltage,  K*,  the 
optimal  amplitude  of  the  synchronous  spatial  harmonic  of 
the  field  scales  as  yj"' ,  where  yo=  1  +eV,/mc^  is  the  elec¬ 
tron  Lorentz  factor.  This  tendency  becomes  distinctive 
enough  at  voltages  above  300—400  kV  In  the  case  of 
BWO’s  the  minus  first  harmonic  is  synchronous.  So,  its 
amplitude  becomes  smaller  as  the  voltage  grows  while  the 
amplitude  of  the  zero  spatial  harmonic  becomes  larger  be¬ 
cause  the  radiated  power  grows  proportionally  to  the  beam 
power.  From  here  it  follows  that  at  high  voltages  the  role 
of  the  electrons’  interaction  with  the  zero  spatial  harmonic 
becomes  more  important,  and  this  interaction  may  partly 
be  caused  by  the  above  discussed  axial  modulation  of  elec¬ 
tron  velocities. 

In  a  general  form  this  effect  was  taken  into  account  in 
Ref.  6  where  the  dispersion  equation  for  the  waves  hia 
waveguide  with  corrugated  walls  was  derived  and  the  cor^ 
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Uniform  axial  magnetic  field 


FIG.  1.  A  schematic  picture  of  a  relativistic  backward  wave  oscillator.  A 
thin  annular  relativistic  electron  beam  of  radius  propagates  inside  of  a 
metallic  structure  with  rippled  walls. 


responding  terms  proportional  to  dv^dz  were  presented. 
However,  no  analysis  of  the  effect  has  been  carried  out. 

The  present  paper  is  aimed  at  studying  the  effect  of 
parametric  coupling  between  electrons  and  the  zero  spatial 
harmonic  of  the  field  in  the  simplest  model  of  a  relativistic 
Cerenkov  device,  the  backward  wave  oscillator  with  large 
end  reflections.  The  paper  is  organized  as  follows.  Section 

II  contains  general  equations  describing  the  effect.  In  Sec. 

III  the  linear  theory  is  presented.  In  Sec.  IV  the  stationary 
nonlinear  theory  is  given.  Finally,  Sec.  V  contains  the  dis¬ 
cussion  and  summary  of  our  results. 

il.  GENERAL  EQUATIONS 


eh  1 


A-/ 


Here  the  initial  electron  velocity  normalized  to  the  speed  of 
light  is  equal  to  As  follows  from  Eq.  (6h 

the  effect  of  energy  depression  is  rather  small.  However,  it 
may  play  an  important  role  in  phase  relations  between 
electrons  and  synchronous  spatial  harmonic  of  the  wave. 
Indeed,  by  introducing,  according  to  Eq.  (1),  the  phase  of 
an  electron  with  respect  to  the  minus  first  harmonic  as 
(f>  =  cot- [ito+  (2ir/d)]z,  where  electrons  z  and  t  are  related 
as  dz/dt=v^,  one  can  obtain  the  standard  equation  for 
d4>/dz: 


d<f>  1  ko+liT/d 
dz  Pi  (o/C 

Here  the  axial  coordinate  z  is  normalized  to  co/c. 

Correspondingly,  by  introducing  the  mean  value  of 
electron  energy,  which  accounts  for  energy  depression 


rd.m.=ro-2 


1  A 


(7) 


(cf  Eq.  6),  the  mean  value  of  the  electron  unperturbed 
normalized  velocity, 

P]0~  Yd-m.’ 

and  the  mean  mismatch  of  synchronism. 


First  of  all,  let  us  consider  the  effect  of  voltage  depres¬ 
sion  on  an  electron  beam  propagating  near  a  corrugated 
metallic  wall.  Let  a  thin  annular  electron  beam  propagate 
at  a  distance  from  the  corrugated  wall  of  the  cylindrical 
waveguide  that  for  ail  z  is  much  smaller  than  the  period  of 
corrugation,  d.  In  such  a  case,  which  corresponds  to  the 
condition 

A,  l<d,  (3) 

where  A  is  the  distance  between  the  average  radius  of  the 
wall,  Tjv,  and  the  electron  beam  radius,  r*,  and  /  is  the 
height  of  corrugation  (see  Fig.  1),  the  axial  modulation  of 
voltage  depression  can  be  considered  locally,  which  leads 
to  the  standard  definition  of  electron  energy  depression 


"  ri, ' 

In  this  equation  is  the  electron  current,  and  the  wave¬ 
guide  radius  is  a  slowly  varying  function  of  z.  For  small 
voltage  depression  we  can  assume  that  in  the  right-hand 
side  of  Eq.  (4),  y=;7o-  Also,  supposing  that  the  sizes  A  and 
/  are  small  not  only  in  comparison  with  d,  as  given  by  Eq. 
(3),  but  also  with  the  beam  radius,  r*,  we  can  expand 
InCr^r*)  in  the  form 


A  / 

sa - —  cos 

rb  n 


(5) 


The  latter  expression  is  given  for  the  wall  corrugation  in 
the  form  r„,=rfr-|- A-/cos(2iTz/<f).  Correspondingly,  the 
depression  of  the  beam  potential  can  be  described  by  a 
simple  equation  for  Ay=yo~7‘ 


1  kQ+ltr/d 

5=’z — - » 

Pa 

one  can  rewrite  the  equation  for  the  phase  <(>  in  the  form 


dz  Pi  Pa 

The  equation  for  the  electron  energy  should  include 
the  periodic  modulation  of  energy  depression  [see  Eq.  (6)] 
and  the  field  of  the  nonsynchronous,  zero  spatial  harmonic. 
The  latter  could  be  important  because,  as  was  mentioned  in 
the  Introduction,  in  relativistic  BWO’s  the  amplitude  of 
this  harmonic  can  be  much  larger  than  the  amplitude  of 
the  synchronous,  minus  first  harmonic.  Therefore,  in  spite 
of  the  lack  of  synchronism,  the  effect  of  the  nonsynchro¬ 
nous  interaction  with  the  large  amplitude  harmonic  can  be 
comparable  to  that  of  the  synchronous  interaction  with  the 
small  amplitude  harmonic.  Correspondingly,  by  represent¬ 
ing  the  axial  component  of  the  electric  field  acting  on  elec¬ 
trons  in  the  form 


£,=Re 


(here  z  is  normalized  to  co/c,  A  is  the  wavelen^h),  one  can 
obtain  for  the  normalized  energy  y  the  following  equation: 

-/(yg- 1  Y sin(Az)  -Re(^e'^(  1  +hLge‘'^)]. 
dz  ^ 

(9) 


Here 
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/= 


eh 


/ 

-L 


is  the  normalized  beam  current  parameter  responsible  for 
axial  modulation  of  electron  energies  and  velocities  near 
the  ripple  structure,  L  is  the  normalized  length  of  a  struc¬ 
ture,  h=2irc/<od  is  the  normalized  axial  wave  number  of 
the  slow-wave  structure,  hL = 2irA^,  where  N  is  the  number 
of  structure  periods,  A = eE_  x  j/mca  xs  the  normalized  am¬ 
plitude  of  the  synchronous  spatial  harmonic,  and  the  ratio 
of  harmonic  amplitudes  is  presented  as  EQjE_\_^=hLq. 
Note  that  since  both  the  ratio  £(>_,/£_  1,2  in  relativistic 
BWO’s  and  hL  are  large,  the  parameter  q—EQjlTrNE_\^p 
which  describes  the  effect  of  a  large  amplitude,  nonsyn- 
chronous  zero  spatial  harmonic,  may  have  an  arbitrary 
value.  Also  note  that  when  the  condition  given  by  Eq.  (  3) 
is  fulfilled,  the  ratio  of  harmonic  amplitudes  Eq^E in¬ 
troduced  by  Eq.  ( 1 )  may  be  found  from  formulas  given  in 
Ref.  7: 


Eo^  2go _ h^P-\’’av) 

J\{gorav)  hiP-\rb)’ 

Here  go  is  the  transverse  wave  number  for  the  fast-wave 
zero  spatial  harmonic  and  p_i  is  related  to  the  transverse 
wave  number  of  the  minus  first,  slow  wave  spatial  har¬ 
monic,  g_„  by  g.i=ip-,.  From  Eq.  (10),  it  follows  that 
the  ratio  of  harmonic  amplitudes  can  be  positive  or  nega¬ 
tive. 

In  comparison  with  the  standard  equation  for  electron 
energy  in  backward  wave  oscillators  with  negligibly  small 
space  charge  forces  (see,  e.g..  Refs.  4  and  5),  Eq.  (9) 
conuins  two  additional  terms  responsible,  respectively,  for 
spatial  modulation  of  electron  energy  depression  and  for 
interaction  with  the  zero  spatial  harmonic  of  the  wave. 

The  boundary  conditions  for  equations  of  electron  mo¬ 
tion  (8)  and  (9)  are  0(0)  =i^o^[*^;2ir]  and  y(0)=!7o-  The 
efficiency  of  interaction  between  the  electrons  and  the  elec¬ 
tromagnetic  field  is 

In  principle,  the  equations  for  the  electron  motion 
should  be  supplemented  with  the  equations  defining  the 
amplitude  of  the  electromagnetic  field,  which  varies  along 
z  when  at  least  one  end  of  the  structure  is  opened  (as 
shown  in  Fig.  1 ) .  However,  below  we  shall  restrict  our 
consideration  to  the  simplest  model  of  a  resonant  Cerenkov 
device  with  strong  reflections  at  both  ends  of  the  interac¬ 
tion  region.  In  such  a  system  the  amplitude  of  the  field 
acting  on  electrons  is  practically  independent  of  z  and  can 
be  determined  by  a  standard  balance  equation  between  the 
microwave  power  radiated  by  an  electron  beam  and  the 
power  of  microwave  diffractive  and  Ohmic  losses. 


III.  LINEAR  THEORY 

To  develop  the  linear  theory  one  has  to  linearize  the 
equations  for  electron  motion  (8)  and  (9)  with  respect  to 


the  field  ampUtude  A.  Zeroth-order  terms  lead  us  to  the 
expression  for  the  electron  energy  given  by  Eq.  (6)  and  for 
the  phase  <f>: 

2eh  ^  I  ■  ('^J\ 

Then,  linearizing  Eqs.  (8)  and  (9)  one  can  find  the  first- 
order  terms  in  y  and  0.  However,  after  averaging  over 
initial  phase  0o  in  Eq.  ( 1 1 )  these  terms  will  not  contribute 
to  energy  exchange  between  the  electron  beam  and  the 
microwave  field.  Therefore,  it  is  necessary  to  take  into  ac¬ 
count  the  second-order  terms  in  perturbations  in  electron 
energy.  The  corresponding  formalism  has  been  described 
in  detail  elsewhere,*  so  we  will  present  here  only  the  final 
expression  for  the  electron  efficiency  defined  in  the  frame  of 
the  linear  theory; 


A^L^ 

’^^'^""2(yo-l)(ro-l)''' 


s>-{e,q,i). 


(13) 


Here  the  function  ^  depends  on  the  electron  transit  angle 
through  the  interaction  region,  6=5L,  on  the  parameter  q 
responsible  for  the  interaction  with  the  zero  spatial  har¬ 
monic  of  the  microwave  field,  and  on  the  normalized  beam 
current  parameter  /  describing  the  effect  of  axial  modula¬ 
tion  of  electron  velocities  by  rippled  walls: 


S^{e,  q,  I)  =  sin  d[  1  ^qd  +  0)  ] 

—  2(1— COS  0)(l-h^/)}* 

This  expression  is  valid  for  any  structure  containing  an 
integer  number  of  periods,  d.  Note  that  among  the  terms 
accounting  for  parametric  resonance  with  nonsynchronous 
spatial  harmonics  only  the  terms  responsible  for  interac¬ 
tion  with  the  large-amplitude  zero  spatial  harmonic  are 
accounted  for  here.  When  both  new  effects,  interaction 
with  the  zero  spatial  harmonic  and  axial  modulation  of 
electron  velocities,  are  negligibly  small  Eq. 

(14)  is  reduced  to  the  standard  gain  function  well  known 
in  the  theory  of  the  “O”  type,  linear  beam  devices.  When 
we  neglect  only  the  effect  of  axial  modulation  of  electron 
velocities,  the  effect  of  interaction  with  the  zero  spatial 
harmonic  still  exists  [see  Eq.  (14)  for  7  =  0,  How¬ 

ever,  if  we  neglect  the  effect  of  interaction  with  the  zero 
spatial  harmonic  (^=0),  the  parameter  /  does  not  play 
any  role  in  Eq.  (14);  that  demonstrates  the  parametric 
nature  of  this  effect. 

In  order  to  illustrate  the  effect  of  the  two  new  param¬ 
eters  presented  in  Eq.  (14),  let  us  give  some  results  of  the 
linear  analysis.  In  Fig.  2(a)  the  standard  linear  gain  func¬ 
tion  is  shown  for  the  case  ^=7=0  as  the  function  of  the 
transit  angle  <9.  In  Figs.  2(b)  and  2(c)  the  lines  of  equal 
ratios  of  gain  functions  ^/^{q—l—^)  are  given  in  the 
plane  of  these  parameters  q  and  7.  Figure  2(b)  corresponds 
to  the  transit  angle  0=  —2.62  which  gives  the  maximum 
linear  gain  function  for  ^=7=0.  Figure  2(c)  corresponds 
to  0=  -5,  where  operation  with  the  high  efficiency  can  be 
realized,  as  will  be  shown  below.  In  the  case  shown  in 
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FIG.  2.  (a)  Standard  linear  gain  function  in  the  absence  of  new  effects. 

(b)  and  (c)  The  lines  of  equal  ratios  of  linear  gain  functions  y{gj)/ 

/— 0)  in  the  plane  of  normalized  parameters  proportional  to  the 
axial  modulation  of  electron  potential  (/)  and  to  the  coupling  with  the 
zero  spatial  harmonic  of  the  microwave  field  {g).  In  (b)  electron  transit 
angle  0=  —2.62,  which  corresponds  to  the  minimum  starting  current;  in 

(c)  0=— 5,  which  corresponds  to  highly  efficient  operation. 


Fig.  2(b)  the  changes  in  the  linear  gain  function  are  small 
while  in  the  case  shown  in  Fig.  2(c),  the  effect  is  much 
more  pronounced.  As  follows  from  Fig.  2(c),  the  mam 
role  is  played  by  the  interaction  with  the  zero  spatial  har¬ 
monic.  At  negative  q  this  effect  reduces  starting  currents 
and  vice  versa. 


FIG.  3.  Contours  of  equal  efficiency  for  a  resonant  BWO  in  the  plane  of 
normalized  amplitude  of  the  microwave  field  and  detuning  of  Cerenkov 
synchronism. 


IV.  NONLINEAR  THEORY 

Numerical  simulations  of  the  nonlinear  set  of  equa¬ 
tions  (8)  and  (9)  have  been  carried  out  for  the  specific 
value  of  electron  initial  energy,  ro=2.  The  electron  effi¬ 
ciency  was  first  studied  in  the  absence  of  the  two  new 


fig.  4.  Lines  of  equal  efficiency  in  the  plane  of  normalized  param«ers  ^ 
and  /  for  the  normalized  length  of  a  slow-wave  structure  L  =  20.  iV-8. 
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FIG.  5.  Axial  dependences  of  electron  velocities  with  and  without  ac¬ 
count  for  axial  modulation  of  electron  potential  and  interaction  with  zero 
spatial  harmonics  of  the  field.  All  figures  correspond  to  L=20, 
5= —0.25.  Other  parameters  are  (a)  >4=0.05,  7=^=0;  (b)  >4=0.05, 
/=0,  ^=-0.1,  iV=8;  (c)  >4=0.05,  7=0.1,  ^=-0.1,  Ar=8;  (d) 
>4=0.07,  7=^=0;  (c)  >4=0.07.  7=0.15,  ^=-0.03. 


effects  (at  ^=7=0)  as  the  function  of  the  mismatch  of 
synchronism,  5,  and  the  normalized  field  amplitude,  A,  at 
different  values  of  the  normalized  length,  L.  It  was  found 
that  when  L  varies  from  15  to  25  the  maximum  efficiency 
is  a  rather  weak  function  of  L,  The  typical  dependence  of 
efficiency  on  6  and  A  is  shown  in  Fig.  3  for  L  =  20  as  lines 
of  equal  efficiency  in  the  plane  of  frequency  detuning,  5, 
and  the  normalized  amplitude,  A. 

Then,  the  effect  of  new  parameters,  q  and  /,  was  stud¬ 
ied.  It  was  found  that  these  parameters  may  give  relatively 
weak  enhancement  in  the  efficiency  of  operation  when  the 
detuning,  6,  and  amplitude,  A,  correspond  to  the  maximum 
efficiency  in  the  absence  of  the  new  effects.  However,  these 
effects  may  widen  the  region  of  parameters  where  efficient 


operation  may  take  place,  especially  at  small  amplitudes. 
Corresponding  results  are  presented  in  Fig.  4,  where 
>4=0.05  and  L  =  20.  In  Fig.  4(a)  the  lines  of  equal  effi¬ 
ciency  are  shown  in  the  plane  of  parameters  q,  7  for  8  — 
—  0.25  and  an  8-period  structure  (results  of  the  study  of  a 
6-period  structure  for  the  same  values  of  A,  L,  and  8  are 
very  similar).  In  Fig.  4(b)  the  same  lines  are  shown  for 
<5=  —0.2  and  an  8-period  structure.  In  both  cases  the  max¬ 
imum  efficiency  is  above  26%  which  is  close  to  the 
77max— 29%  shown  in  Fig.  3  for  the  system  without  these 
new  effects. 

To  illustrate  the  effect  of  velocity  modulation  by  inter¬ 
action  with  the  zero  spatial  harmonic  of  the  wave  and  a 
rippled  wall  structure,  it  is  expedient  to  present  the  axial 
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dependencies  of  electron  velocities  with  and  without  these 
effects.  Such  dependencies  are  shown  in  Fig.  5.  Figures 
5(a)-5(c)  correspond  to  the  field  amplitude^  =0.05;  Figs. 
5(d)  and  5(e)  to  v4=0.07.  All  other  parameters  are 
L=20,  6= -0.25,  N=&.  Figures  5(a)-5(d)  correspond 
to  g=I=0-,  Fig.  5(b)  to  1=0,  9= -0.1;  Fig.  5(c)  to 
7=0.1,  9= —0.1;  and  Fig.  5(e)  to  7=0.15,  ^=— 0.03. 
Comparing  these  figures,  one  can  find  that  the  interaction 
of  electrons  with  the  zero  spatial  harmonic  and  the  axial 
modulation  of  the  electron  potential  have  a  significant  ef¬ 
fect  on  the  reduction  of  electron  velocities.  In  particular, 
Fig.  5(e)  shows  uiat  the  operation  at  given  parameters  is 
close  to  the  appearance  of  reflected  particles  that  can  also 
be  stimulated  by  an  initial  electron  velocity  spread,  which 
was  not  taken  into  account  here. 


V.  DISCUSSION 

To  illustrate  the  results  of  the  theory  developed  in  the 
normalized  parameters  above,  let  us  express  these  param¬ 
eters  in  real  numbers  corresponding  to  a  typical  experi¬ 
mental  setup.  Consider  the  8-period  slow-wave  structure 
described  in  Ref.  9  and  used  in  a  number  of  experiments 
with  relativistic  BWO’s  at  the  Institute  for  Plasma  Re¬ 
search,  University  of  Maryland.  This  structure  operated  in 
the  8  GHz  frequency  range.  The  period  of  the  structure  is 
1.67  cm,  its  average  radius  is  1.5  cm  and  the  height  of 
corrugation  is  0.4  cm.  Correspondingly,  if  we  take  the  axial 
wave  number  for  the  zero  spatial  harmonic  as  equal  to 
—  -rr/2d,  then  for  an  electron  beam  with  500  k V  operating 
voltage  and  1  cm  beam  radius  we  get  that  the  value  of 
normalized  current  parameter  0.05  corresponds  to  —0.5 
kA  of  the  real  electron  beam  current. 

In  the  same  manner  one  can  find  from  Eq.  (10)  that 
parameters  of  the  slow-wave  structure  and  the  electron 
beam  given  above  correspond  to  the  ratio  of  harmonic  am¬ 


plitudes  Eq/E^^  ,=^-\0.6,  which  in  turn  corresponds  to 
the  value  of  the  normalized  parameter  9=— 0.236. 

This  example  shows  us  that  the  results  presented  in  the 
paper  can  be  used  for  interpretation  of  a  number  of  exper¬ 
iments. 

Note  in  conclusion  that  the  effects  studied  may  lead, 
first,  to  achieving  a  highly  efficient  operation  at  low  beam 
currents  and,  second,  to  a  more  smoother  dependence  of 
the  BWO  efficiency  on  the  beam  current  than  that  pre¬ 
dicted  by  the  theory  developed  without  these  effects.  Be¬ 
sides  the  efficiency  increase,  these  effects  may  also  cause 
significant  reduction  in  efficiency  at  certain  parameters. 
Also  important  are  pulsations  in  electron  axial  velocities 
caused  by  these  effects  because  such  pulsations  may  lead  to 
the  appearance  of  reflected  electrons  in  high  efficiency  op¬ 
eration  regimes. 
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Plasma  Source 

Plasma  is  formed  when  a  high  voltage  pulse  is 
applied  between  the  center  and  outer  conductor 
of  gun.  Ionized  gas  is  ejected  from  titanium  hydride 
impregnated  rings  concentric  with  the  center 
conductor. 
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End  View: 


Graphite  — 


Stainless  steel  outer 
jacket/  ground  connection 


TiH2in  Sodium 
Silicate 


Center  conductor 


A  system  for  mapping  the  plasma  density  inside  the  RF 
interaction  region  of  a  plasma  loaded  relativistic  BWO. 
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Probe  Bias  Circuits  for  Plasma  Measurement: 


Pulse  Bias  Circuit: 


To  probe 


D.  C.  Bias  Circuit: 


Ground  on  Vacuum  Vessel 
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Measuring  the  Plasma  Density 
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Radial  Profile  of  Probe  Current 


Probe  Position:  25  cm  Gun  Voltage:  12  kV 
Gun  Position:  32  cm  Magnetic  Field:  6.0  kC 


Resonant  Cavity 

Plasma  density  may  also  be  measured  by  determining  the  amount 
of  shift  in  the  resonant  frequency  of  a  simple  cavity  due  to  the 
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numbers  observed  to  the  predictions  from  a  numerical  model  for  the 
system.  The  details  of  this  method  were  presented  in  session  1D03, 
by  S.  Kobayshi. 


Schematic  Diagram  of  Resonant  Cavity 
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Measuring  the  Electromasnetic  Characteristics 
of  Plasma  Loaded  Cavities 
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Langmuir  probe  was  used  to  find  the  peak  plasma  density 
and  the  radial  profile  of  the  plasma  column. 

Cavity  technique  provided  an  independent  calibration 
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Future  experiments  with  the  plasma-loaded  backward  wave 
oscillator  can  now  select  a  particular  plasma  density  through 
timing  control  of  the  plasma  gun. 
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Starting  Energy  and  Current  for  a  Large 
Diameter  Finite  Length  Backward  Wave 
Oscillator  Operated  at  the  Fundamental  Mode 

K-  Minami,  K.  Ogura,  Y.  Aiba,  M.  R.  Amin,  X.  D.  Zheng,  T.  Watanabc,  Y.  Cannel,  Senior  Member,  lEEE^ 

W.  W.  Destler,  Fellow,  IEEE,  and  V.  L.  Granatstein,  Fellow,  IEEE 


Abstract — ^We  study  the  startiiig  condhloiis  for  a  large  diameter 
(diameterAraveleiigth  »  4.8)  finite  length  backward  wave  oscil* 
lator  designed  for  24*GHz  operation  at  the  fundamental  TMoi 
mode.  This  geometry  is  very  promising  for  high  power  handling 
capability.  We  analyze  two  separate  threshold  conditions.  First, 
finite  length  effects  give  rise  to  a  threshold  in  electron  beam 
energy  below  which  oscillations  cannot  be  sustained  at  any  beam 
current  The  second  is  the  more  familiar  current  threshold  known 
as  a  start  current  It  is  also  found  that  the  growth  rate  for  the 
fundamental  mode  can  be  much  hurgo'  than  those  of  other  higher 
order  modes  thus  leading  to  coherent  operation  of  large  diameter 
sources  free  from  mode  competition. 

I.  Introduction 

IGH-POWER  MICROWAVE  sources  are  important  for 
a  number  of  advanced  applications  ranging  from  current 
drive  and  RF  heating  of  magnetically  confined  plasmas  in 
fusion  devices  to  high  resolution  nanosecond  radars  (l]-[3]. 
Pulsed  high  power  microwave  oscillators  utilizing  intense 
relativistic  electron  beam&have  been  extensively  studied  in 
recent  years  [2]-[6].  Among  various  microwave  sources,  the 
multiwave  Cerenkov  generators  (MWCG’s)  developed  by  [7] 
at  the  High  Current  Electronics  Institute,  Tomsk,  Russia,  have 
recently  attained  record  outputs  [7],  [8].  Radiation  powers  of 
7.5  GW  at  a  wavelength  A  =  9.7  mm  with  an  electronic 
efficiency  of  20%  were  reported.  The  MWCG  involves  a 
slow  wave  structure  (SWS)  with  an  average  diameter  D  much 
larger  than  A,  In  the  above  example,  D/\  =  13.  This  large 
diameter  SWS  enables  a  larger  output  power  for  a  given 
RF  energy  density  inside  the  SWS  before  breakdown  occurs. 
Despite  a  possibility  of  mode  competitions,  efficient,  high 
power  single  mode  oscillation  was  attained  in  a  device  utilizing 
overmoded  SWS  [9].  Twenty  years  ago,  high  average  power 
millimeter  microwave  sources  (gyrotrons)  were  invented  by 
combining  an  electron  cyclotron  maser  and  an  overmoded 
open  barrel  cavity  [2],  [3].  The  large  diameter  (overmoded) 
SWS  in  MWCG’s  may  be  an  innovative  key  technology 
which  corresponds  to  the  overmoded  open  cavity  employed 
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in  gyrotrons,  although  the  final  evaluation  of  performance  and 
the  feasibility  of  using  MWCG’s  in  actual  applications  are 
yet  to  be  established.  The  MWCG  has  a  close  relationship  to 
conventional  backward  wave  oscillators  (BWO’s),  since  both 
employ  SWS’s  in  common  [4],  [5],  [7H14]. 

Although  the  physical  processes  involved  in  the  MWCG’s 
[7],  [8]  are  complicated,  we  point  out  here  that  high-power 
millimeter  microwaves  can  be  generated  without  decreasing 
the  mean  diameter  of  the  SWS.  Both  BWO’s  and  large 
diameter  BWO’s  can  be  operated  in  the  fundamental  TMoi 
mode  as  will  be  shown  in  this  paper.  The  difference  in  SWS’s 
between  the  two  devices  is  as  follows.  The  inner  radius  of 
the  metal  surface  of  the  SWS,  R{z),  is  assumed  to  vary 
sinusoidally  around  the  mean  radius  Rq,  i.e.,  R{z)  =  i2o  + 
h  cos  Kqz,  Ko  =  2'ir/zo.  Here,  zq  and  h  are,  respectively,  the 
axial  length  of  periodicity  and  the  amplitude  of  corrugation.  In 
BWO’s  [4],  [5],  [10H14],  the  mean  radii,  Rq  =  D/2,  were 
chosen  such  that  D/X  =  2Ro/X  ^  1,  whereas,  in  the  present 
large  diameter  BWO,  the  condition  D/X  ^  1  is  explored. 
The  key  point  is  that,  in  the  latter,  oscillation  frequency  is 
raised  by  carefully  choosing  small  values  of  h  and  zq,  keeping 
Rq  larger  than  A.  This  situation  is  quite  similar  to  the  SWS 
design  of  MWCG’s  [7]-[9].  In  the  present  paper,  numerical 
studies  are  made  within  the  scope  of  linear  analysis  for  a 
large  diameter  SWS  with  D/X  -  4.8.  Specifically,  the  starting 
energy  and  the  starting  current  of  the  beams  for  initiating 
microwave  oscillation  in  such  a  finite  length  SWS  are  analyzed 
in  detail.  Previously,  the  starting  current  conditions  for  finite 
length  BWO’s  have  been  studied  by  various  authors  [5],  [13], 
[15],  [16].  Sweglc  analyzed  the  starting  condition  for  a  BWO 
[15].  In  his  treatment,  100%  reflection  at  the  beam  entrance 
boundary  and  0%  reflection  at  the  RF  output  boundary  were 
assumed.  This  assumption  was  valid  for  the  case  of  matched 
termination  at  the  output  end  of  the  SWS.  We  consider  here 
an  alternative  case  of  nonzero  round  trip  reflection  which  may 
correspond  to  the  experimental  situations  in  [4]  and  [10],  Such 
nonzero  round  trip  reflection  produces  discrete  spectra  of  axial 
wavenumber  of  eigenmodes  of  an  axisymmetric  TM  mode 
[13],  [14],  [17]. 

In  Section  n,  we  describe  the  model  of  our  linear  analysis. 
The  difficulties  to  be  overcome  in  the  practical  numerical 
computation  in  the  case  of  D/X  ^  1  are  discussed.  Numerical 
results  are  shown  in  Section  HI.  Discussion  and  conclusions 
are  given  in  Section  IV. 
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H  PMtMUiAnnrAMD  CMmcuiaiEs  ^ 

IN  Numerical  Analysis 

The  model  ot  iit  exi^nuBeiiie  laige-  diemetef  SWS  with 
length  L  shown  in  Rg.  1  is  analyzeil  The  radiiis  iJ(2!)  of 
perfectly  conducting  metal  wall  varies  with  R{z)  =  Ro  + 
hcoeKoz.  The  entire  system  is  immersed  in  a  strong  longitudi¬ 
nal  external  magnetic  field.  The  transverse  motion  of  the  beam 
electrons  is  to  be  negligible.  We  also  assume  that  an 

infinitely  thin  annular  beam  With  radius  Ri,,  monoenergetic 
beam  energy  14  and  current  A  is  incident  to  the  SWS.  The 
Mnrtar  dispersioD  relation  for  the  system  in  Rg.  1  was  derived 
previously  [15],  [18],  [19]  and  the  results  are  summarized  here. 
The  boundary  condition  that  the  tangential  components  of  the 
RF  electric  field  at  the  corrugated  metal  surface  of  the  SWS 
must  be  zero  gives 

[1  +  (n  -  +  CnC^n) 

m,  n=— oo 


CORRUGATED  WALL 
WAVEGUIDE 
R(2)-Ro+h  cos  KoZ 


ANNULAR  BEAM 
ENERGY  Vb 
CURfOIT  lb 


—  Dmn  ■ 

=  0  (1) 

Cmn  =  /  (*(”  “  rn)Koz] 

J-k/Ko 

■  A(yn(l  +  Q  cos  Koz)]  dz  (2) 

2^  g!(g  -h  |n  —  m|)! 

a  ^  h/Ro  (3) 

4 

where  Qn  =  Koknliui'^lc^  -  ^n)>  =  k  +  nKo,  Vn  = 

!<?  -  kn)  and  c  is  the  vacuum  velocity  of  light.  The  RF 
field  amplitudes,  A„,  5„,  and  Cn  of  the  Roquet  harmonics  are 
the  same  as  those  in  (7)  in  [18].  The  beam  quantities  such  as  A 
and  14  are  involved  in  Bn  and  C„.  In  (1),  C7^„  is  obtained  by 
replacing  the  Bessel  function  of  the  first  kind,  Jo,  in  by 
that  of  the  second  kind.  No-  In  the  above  expressions,  temporal 
and  spatial  phase  factors  exp  -wt)]  have  been  assumed 
for  every  RF  quantity.  The  Aspersion  relation  combining  w 
and  k  is  given  from  0)  as  the  determinant  equation, 

D{k,  u)  =  det  [Bmn]  =  0  (4) 


some  difficulties  for  computing  (2)  and  (3).  In  order  to  get 
a  large  oscillation  frequrmcy  a»/2ir  without  decreasing  Ro 
we  must  choose  a  small  h  and  zo  in  Pi?-  I  1^***  results  ir 
a  large  Ko  and  that  makes  <  q  in  (2).  For  negativf 
yl,  the  Bessel  function  Jo  becomes  modified  Bessel  toctioi 
A  consequently,  and  the  integrand  of  (2)  is  rewritten  a 
exp  [i(n  — Tn)/ifo'2]A(yn(^4'‘^  -Ko^)]! '^^lere  y^  —  — y„ 

The  modified  Bessel  function  A  in  the  integrand  become 
extremely  large,  if  y(,  is  large.  This  effect  causes  an  overflow  i; 
the  process  of  numerical  computation  of  (4).  Such  a  difficult 
arises  often  in  the  calculation  of  the  large  diameter  SWS  witi 
small  zo-  To  avoid  the  difficulty,  we  express  Iv{z)  as 

r 

~  (27rz)i/2  8z  2!(8z)2  J 

|arg{z)l<|,  = 

Using  this  asymptotic  expansion,  (1)  is  expressed  as 

mt  n*-*oo 

G„  =  A„  exp  (y^). 

D'{k,  oj)  =  det  [IC„1 


which  is  required  to  have  nontrivial  An  in  (1).  As  is  well 
known.  No  goes  to  infinity  at  the  origin  which  occurs  in  the 
case  of  light  line  in  vacuum,  w  =  ckn.  Because  the  roots  of 
(4)  are  slow  waves,  they  are  located  in  the  complex  w  and  k 
planes  considerably  seps’^  ad  from  where  w  =  ckn-  Difficulty 
of  No  infinity  does  not  3ur  for  finding  the  roots. 

We  discuss  here  some  peculiar  difficulties  that  arise  in 
the  computation  for  the  large  diameter  SWS.  For  the  time 
being,  we  consider  the  simple  case  without  the  beam,  i.e., 
A  =  0,  and  An  =  Bn,  Cn  =  0  in  (1)  accordingly  [18]. 
Because  no  input  energy  is  present  in  the  SWS,  oscillations 
caimot  be  expected.  For  a  real  w,  real  k’s  are  obtained  in 
(4)  and  vice  versa.  Equations  (2)  and  (3)  were  useful  to 
compute  the  disposion  relation  of  a  BWO  [19],  [20].  In  the 
present  case  of  large  diameter  BWO’s,  however,  fliete  exist 


=  0.  (f 

The  dispersion  relation  is  given  by  (5)  instead  of  (4).  Equado 
(5)  does  not  have  problems  with  numerical  overflow. 

The  next  problem  concerns  the  adequacy  of  the  Taylc 
used  in  (3).  In  general,  numerical  integration  of  C 
takes  much  computation  time  and  the  expanded  form  of  (3) : 
preferable,  as  long  as  the  sinusoidally  corrugated  SWS  inFij 
1  is  considered.  A  large  Ko  again  may  result  in  |yT>a|  > 
in  (3),  and  the  summation  does  not  converge  rapidly.  We  fin 
that  tens  of  terms  in  (3)  must  be  calculated  to  get  a  converge 
value  of  summation  because  of  increased  denOTiinators.  1 
avoid  this  difficulty,  we  truncate  the  rank  of  die  determinant  L 
(5)  to  as  small  a  value  as  possible.  This  is  because  the  value 
of  kn  and  accordingly  Vn  will  not  be  very  large  for  sma 
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Arirnmnam  Por  exmpte^  \jf^a\  >  1  may 
iimita  *-1  <n,m<0(2x2  detemrinant) 
qr  <n,m<l(3x3  detensmant)!  Hg.  2  shows  the 
i^^  iesallof  die  design  study;  die  oadlladoo  frequency  /  =  u)/2ir 
versus  the  corrugation  wavenumber  Kq  for  several  corrugation 
amplitudes  h.  Here,  iZo  3.0  cm  has  been  assumed  as  a  given 
parameter.  Curves  in  the  figure  are  obtained  by  computing  the 
frequencies  at  the  crossing  point  of  the  structure  mode  and 
the  beam  space  charge  line  with  =  0  for  various  values  of 
Ko  and  h.  The  straight  line,  w  =  Vi,Kol2,  in  Fig.  2  is  the 
boundary  that  divides  the  BWO  (absolute  instability)  from  the 
TWT  (convective  instability)  for  the  case  of  7r-mode  operation 
{k  =  K(s/2)  with  14  =  100  keV.  Microwave  oscillation  can  be 
expected  in  the  cases  of  Kq  to  the  left  of  the  straight  boundary 
line.  Roughly  speaking,  BWO’s  and  large  diameter  BWO’s 
arc  expected,  respectively,  in  cases  of  Ko  <  10  cm~^  and 
Kq  >  10  cm~^  The  black  circle  represents  the  parameters 
we  have  chosen  for  the  present  design  study  of  a  large 
diameter  BWO  with  the  following  parameters:  Kq  =  18.5 
cm"^  (zq  =  0.34  cm),  /i  =  0.17  cm  and  /  =  24  GHz. 
The  solid  and  dashed  curves  are,  respectively,  the  calculations 
using  (3)  and  (5)  of  the  9  x  9  and  2x2  approximate  truncated 
determinants.  The  differences  between  both  curves  are  less 
than  4%.  For  the  chosen  size  parameters,  the  dispersion  curves 
arc  calculated  for  various  truncated  ranks  from  9  (9  x  9)  to 
2  (2  X  2)  of  the  determinant  in  (5).  For  ranks  larger  than 
6.  the  oscillation  frequencies  are  almost  unchanged,  and  the 
results  are  considered  to  be  exact.  Even  in  the  case  of  rank  2, 
the  deviations  from  exact  values  are  less  than  several  percent. 
These  facts  suggest  that  the  truncation  to  the  rank  2  of  the 
determinant  in  (5)  is  almost  correct  for  our  purpose.  Hence, 
in  order  to  save  computation  time,  we  adopt  the  truncated 
determinant  of  rank  2  in  the  subsequent  numerical  analysis 
of  complex  u  and  k  with  incident  beams,  at  the  sacrifice  of 
accuracy  within  a  few  percent. 

in.  Numerical  Results 

A.  Infinitely  Long  Slow  Wave  Structure 

The  parameters  of  the  large  diameter  SWS  used  in  the 
following  numerical  analysis  arc:  Rq  =  3.0,  h  =  0.17, 
Zq  =  0.34,  and  Rh  =  2.63  cm.  Beam  energies  of  14  =  100 
and  65  keV  are  used  as  typical  values.  The  beam  current  h  is 
mostly  assumed  to  be  0.4  kA,  unless  specified  otherwise. 

The  dispersion  relation  (5)  for  complex  a;/27r  versus  real 
wavenumber  k  is  shown  in  Fig.  3  for  two  values  of  14, 100  and 
63  keV  represented,  respectively,  by  solid  and  dashed  curves. 
Real  and  imaginary  parts  of  u;/27r  are  shown,  respectively, 
in  (a)  and  (b).  Here,  the  beam  current  of  h  =  0.4  kA  has 
been  assumed.  Because  of  FIoquet*s  theorem,  the  dispersion 
curves  have  periodicity  Kq  for  wavenumber  fc,  in  other  words, 
u){k)  =  uj{k  -h  uKq),  Usually,  figures  are  shown  for  the  first 
Brillouin  zone  —KqI2  <  k  <  KqI2,  however,  we  depict 
hereafter  the  curves  for  0  <  fc  <  /iTo,  because  the  important 
parts  of  the  curves  arc  located  near  k  —  Kq  /2.  As  were  pointed 
out  previously,  there  exist  generally  4  independent  roots  for 
the  fundamental  TMoi  mode  [18],  [20].  They  are  the  fast 
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Fig.  2.  Design  study  of  the  slow  wave  stnicture.  Oscillation  frequency  for 
100  keV  beam  versus  corrugation  wavenumber  Kq  for  various  cocnigation 
amplitudes  h.  The  mean  radius  is  Ro  =  3.0  cm.  The  solid  and  dashed  curves 
are.  respectively,  9x9  and  2x2  approxinoate  results  of  the  detennmant  of 
the  dispersion  relation.  The  black  circle  is  the  chosen  parameters  of  our  large 
diameter  backward  wave  oscillator. 


and  slow  space  charge  waves  and  the  backward  and  forward 
structure  waves.  The  first  wave  is  always  heavily  damped  and 
it  is  ignored  entirely  in  the  subsequent  analysis.  The  second 
wave  is  the  energy  source  for  the  third  wave,  i.e.,  the  output 
radiation.  The  fourth  wave  can  serve  as  a  positive  feedback 
mechanism  in  the  case  of  a  finite  length  SWS.  For  simplicity, 
the  fourth  wave  is  also  ignored  hereafter.  For  a  real  A,  (5) 
can  have  a  pair  of  complex  conjugate  roots  of  cj.  They  are  a 
growing  slow  space  charge  wave  and  a  decaying  fast  space 
charge  wave.  The  single  beam  line  for  infinitesimal  currents 
is  split  into  the  fast  (shown  by  F)  and  slow  (shown  by  S) 
beam  space  charge  waves  in  Fig.  3(a),  because  of  nonzero 
beam  current  Absolute  instability  can  be  found  around  the 
crossing  point  of  the  backward  structure  wave  (shown  by  B) 
and  the  line  S.  It  is  noted  in  Fig.  3  that  the  ranges  of  k  for 
instability,  i.e.,  complex  u)  shown  by  thick  solid  and  dashed 
lines  in  Fig.  3(b),  are  much  more  limited  than  the  previous 
case  of  BWO’s  [18],  [19].  This  is  because  the  corrugation 
parameter  a  =  h/Ro  is  presently  smaller  than  those  in  the 
conventional  BWO’s.  In  other  words,  the  oscillation  condition 
in  laige  diameter  BWO’s  with  D/A  >  1  is  more  restrictive 
or  stringent  than  those  in  BWO’s  with  D/A  ^  1. 

In  our  large  diameter  (overmoded)  SWS,  there  exists  a 
possibility  of  mode  competition  between  various  candidates 
including  the  fundamental  and  higher  modes.  Steady  oscilla¬ 
tion  of  the  fundamental  mode  may  not  be  realized  because 
the  beam  energy  may  be  fed  to  various  modes  with  higher 
frequencies.  In  order  to  clarify  the  situation,  we  extended  the 
analysis  for  the  fundamental  mode  shown  in  Fig.  3  up  to  sixth 
higher  modes.  The  oscillation  frequency  and  the  maximum 
temporal  growth  rate  of  the  fundamental  TMoi  mode  are 
compared  with  those  of  higher  modes.  We  computed  the 
dispersion  curves  similar  to  Fig.  3  for  each  higher  mocte.  Fig.  4 
shows  the  dispersion  curves  for  an  infinitesimal  beam  cunent. 
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Fig.  3.  Dispersion  relation  for  complex  uj/25r  versus  real  wavenumber  fc. 
4  =  0.4  kA,  Solid  curves  are  for  Vi  =  100  keV  and  dashed  curves  are  for 
V(,  =  63  keV.  (a)  /  =  Re  (u))/2x  versus  k.  (b)  /m(vj)/2)r  versus  fc. 

from  which  the  oscillation  frequencies  can  be  estimated.  It  is 
interesting  to  note  that  the  dashed  beam  lines,  100  and  65 
keV,  intersect  the  curve  of  the  fundamental  TMqi  mode  at  the 
positions  where  the  group  velocities  are  much  smaller  than 
those  for  higher  modes.  The  results  from  complex  analysis 
like  Fig.  3  are  summarized  in  Table  1,  where  the  maximum 
temporal  growth  rate,  corresponding  oscillation  frequency  and 
oscillation  wavenumber  are  listed  for  six  modes.  Here,  the 
beam  current  0.5  kA  has  been  assumed.  As  is  clearly  shown 
in  the  table,  the  maximum  temporal  growth  rates  of  the 
higher  modes  are  much  smaller  than  that  of  the  fundamental 
TMoi  mode.  Accordingly,  oscillations  of  higher  modes  may 
be  ignored  in  the  scope  of  the  present  linear  analysis.  Of 
course,  mode  competitions  may  be  expected  in  large  diameter 
BWO’s  at  the  saturated  level  of  oscillation  in  general.  At 
the  linear  stage,  however,  mode  competitions  are  insignificant 
in  the  large  diameter  SWS,  if  the  size  and  beam  parameters 
are  carefully  chosen.  For  this  reason,  we  limit  ourselves  the 
subsequent  analysis  to  the  fundamental  TMoi  mode. 

The  dashed  curves  in  Fig.  3  are  for  the  case  of  Vb  =  65  keV . 
The  temporal  growth  rate  Im  (w)  for  Vi  =  65  keV  shown  by 
the  dashed  curve  in  Fig.  3(b)  is  considerably  smaller  than  that 
for  Vb  =  100  keV  shown  by  the  solid  curve.  This  difference 
suggests  that,  like  conventional  BWO’s,  there  exists  a  starting 
energy  threshold  for  oscillation  in  a  finite  length  large  diameter 
SWS  as  will  be  analyzed  in  the  next  subsection. 

In  Fig.  3,  the  real  wavenumber  k  has  been  assumed.  The 
real  k  means  that  we  have  assumed  a  sinusoidal  small  origin 
of  oscillation  with  infinite  extent  in  the  z  direction  and 
have  calculated  its  ten^xwal  evolution.  On  the  other  hand,  in 
the  case  of  a  localized  small  origin  of  disturbance  for  the 
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Fig  4  The  dispersion  curves  for  TMoi ,  •  ■  TMos  modes.  The  two  dashec 
lines  are  the  beam  lines  for  100  and  65  keV  with  infinitesimal  currents 
respectively.  The  crossing  points  of  the  beam  lines  and  smictuie  curves  give 
possible  oscillation  ficequeacies. 


TABLE  I 

Comparison  of  Oscillation  of  Six  Modes  Including  the 
Fundamental  TMoi  Mode.  The  Osciu-ation  Wavenumber  and 
Frequency  at  the  Maximum  Temporal  Growth  Rate  are  Listed  for 
Our  Designed  Large  Diameter  Slow  Wave  Structure.  Beam 
Energy  and  Current  are.  Respecttvely  100  kbV  and  0.5kA 


instability,  the  origin  develops  asymptotically  with  tempon 
and  spatial  factors  exp  -  k,z)],  as  was  show 

in  (2.22)  in  [21].  Here,  u,  and  fc,  are,  respectively,  th 
complex  angular  frequency  and  complex  wavenumber  at  th 
saddle  point  du./dk,  =  0  in  (5).  The  localized  origin  ca 
monotonically  grow  up  in  time  at  every  point  in  z.  This  ^ 
therefore  an  absolute  instability.  Once  an  unstable  root  of  (! 
in  Fig.  3  is  found,  it  is  not  very  difficult  to  access  the  sadd! 
point,  w,  and  fc„  which  exists  uniquely  for  a  given  set  of  a 
parameters,  by  using  the  Newton-Raphson  technique.  In  th 
saddle  point  analysis  of  the  large  diameter  BWO,  it  is  four, 
that  the  oscillation  frequency  Re  (w,)/2ir  and  the  temper; 
growth  rate  Im  (w.)  are,  respectively,  sU^tly  decreasing  an 
increasing  functions  of  /s,  and  the  spatial  growth  rate  Im  (fc 
is  almost  unchanged  with  increase  in  h- 

B.  Finite  Length  Slow  Wave  Structure 

In  a  large  diameter  SWS  with  finite  lengffi  L  as  shown  i 
Fig.  1,  the  wave  reflections  (or  leakages)  at  both  ends  ai 
taken  into  account.  The  end  reflections  result  in  a  fe^bac 
mechanism  by  the  backward  structure  wave,  and  the  distinc 
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dm  betwm  abtcdutB  and  cmanctive  iMtabili^  becomes 
somewbat  andngnous  in  the  £r  finite  case*  The  complex 
wavenumbers  (rf  die  slow  space  diaige  wave  and  the  backward 
stnictuxe  wave  are  denoted,  respectively,  by  a^.  and  a..  In  the 
limit  of  infinite  L,  a+  (=  a.)  coincides  with  the  saddle  point 
of  (5).  For  L  finite  cases,  the  following  equations  must  be 
satisfied  instead  of  (S),  [20]. 


D(a+,u)  =  0 

(6) 

D{a-,  w)  =  0 

(7) 

R  exp  {— t(a_  -  a+)L]  =  1. 

(8) 

Equation  (8)  comes  from  the  requirement  that  the  electromag¬ 
netic  field  must  be  a  single  value  at  every  z  point,  when  the 
propagating  wave  comes  back  after  one  round  trip  [20].  In 
(8),  R  is  the  one  round  trip  reflection  coefficient  of  wave 
fields  at  the  both  ends  of  the  SWS.  and  assumed  to  be  a  given 
parameter.  The  length  L  is  assumed  to  be  L  =  70zo  =  23.8 
cm.  In  general,  the  reflection  coefficient  R  in  (8)  must  be  a 
complex  number,  however,  for  simplicity  it  is  approximately 
replaced  by  a  real  number  0  <  <  1  in  the  subsequent 

analysis.  This  is  because  the  argument  of  the  complex  number 
causes  only  changes  in  equivalent  structure  length  less  than 
zo  much  smaller  than  L.  The  additional  requirement  given  by 
(8)  results  in  the  existence  of  the  starting  energy  in  addition 
to  a  nonzero  starting  current  in  finite  L  BWO’s.  Equation  (8) 
is  rewritten  as 

Re  (a_  -  a+)  =  >  0  (9) 

Im  (a_  -  a^)  =  -  In  {R)IL  >  0  (10) 

where  is  an  integer.  The  following  analysis  is  limited  to 
the  case  =  -1  in  (9),  because  A/  =  -1  is  the  easiest  case 
for  initiating  oscillation  among  various  choices  of  N  for  a 
given  L. 

To  solve  (6)-{8)  correctly,  it  is  necessary  to  watch  the 
movements  of  roots  a+  and  a_  on  the  complex  k  plane  for 
various  u;’s  with  positive  Im  (u;)  approaching  to  zero  [20], 
[21].  This  is  because  the  physical  meaning  (stable  or  unstable) 
of  the  each  root  is  lost  in  the  course  of  numerical  calculation 
of  (6)-(8).  We  must  reconsider  the  physical  meaning  through 
watching  the  movement  of  the  root  on  the  complex  k  plane. 
If  the  panicular  root  traverses  the  real  axis  during  the  change 
in  Im  (a;)  from  a  large  positive  number  to  zero,  the  root  is 
unstable,  otherwise  it  is  stable.  A  pair  of  solutions  a+  and  a- 
is  shown  on  the  complex  k  plane  in  Fig.  5(a)  for  the  case  of 
Vh  =  100  keV.  The  locations  of  the  roots  a+  and  are  found 
at  the  centers  of  the  contour  mapping  of  |i9'|  given  by  (6)-(8) 
on  the  complex  k  plane.  The  corresponding  common  value  of 
complex  u;/27r  =  24.103  +  z0.063  is  found  in  Fig.  5(b)  for 
both  a+  and  a.  shown  by  black  circles  on  the  same  complex 
k  plane  as  Fig  5(a).  The  white  circles  in  Fig.  5(b)  arc  the 
locations  of  Im  (u;)  =  0  that  give  the  boundary  lines  between 
BWO’s  (Im  (u;)  >  0),  and  TWT’s  (Ira  (w)  <  0).  The  arrows 
on  the  solid  lines  show  the  directions  of  decrease  in  Im  (a;)/27r 
toward  zero  for  constant  Re  (u;)/27r.  The  diamond  is  the 
saddle  point  of  (5).  In  L  infinite  case,  oscillation  (absolute 
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Fig.  5.  Locations  of  the  pair  of  roots,  a+  and  a^,  on  the  complex  k  plane 
for  Vb  =  100  keV  and  h  =  0.4  kA,  L  =  23.8  cm,  iV  =  -1  and 
=  24.103+10.063.  (a)  Contour  mapping  of  \D\  where  D  is  defined  by 
(5).  (b)  Locations  of  the  roots  (black  circles)  that  arc  in  the  oscillation  region 
within  border  lines  of  Im(u;)  >  0  shown  by  white  circles.  The  diamond  is 
the  saddle  point. 


instability)  is  expected  at  the  saddle  point  as  was  stated  in  the 
previous  subsection,  if  it  is  located  in  Im  (u^)  >  0  region  in 
complex  k  plane.  In  this  case,  oscillation  is  expected  even  for 
infinitesimal  beam  currents  and  no  limitations  exist  for  starting 
energy.  In  L  finite  case,  however,  the  oscillation  happens  not  at 
the  saddle  point  but  at  a+  obtained  from  (6)-(8).  Oscillation  is 
impossible,  if  the  pair  of  the  roots  are  located  in  Im  (w)  <  0 
region.  Because  the  black  circles  are  located  in  the  region 
Im  (cj)  >  0  in  Fig.  5(b),  the  oscillation  can  be  expected  in 
the  present  case  of  Vt,  =  100  keV.  Since  we  have  chosen 
L  =  23.8  cm  and  AT  =  ~1,  Re  (a-  -  a+)  =  0.264  cm”^ 
from  (9)  is  a  fixed  value,  which  is  independent  of  the  value 
of  reflection  coefficient  R.  The  horizontal  distance  between 
a+  and  a-  is  a  constant  0.264  cm"^  on  the  complex  k  plane 
shown  in  Fig.  5(b). 

A  pair  of  solutions  a+  and  a»  in  (6)-(8)  in  the  case  of 
Vf,  =  65  keV  is  shown  on  the  complex  k  plane  in  Fig.  6.  Other 
parameters  are  the  same  as  those  in  Fig.  5.  The  corresponding 
common  value  of  complex  a;/27r  =  23.295  —  i0.092  is  found 
in  Fig.  6  for  both  a+  and  a_  shown  by  black  circles.  In 
this  case,  and  a-  are  located  in  the  region  of  TWT’s, 
Im  (w)  <  0,  and  oscillation  cannot  occur  for  VJ,  =  65  keV, 
although  the  beam  current  h  =  0.4  kA  is  identical  to  that  for 
Vb  =  100  keV  in  Fig.  5.  Comparing  Figs.  5(b)  and  6,  one  finds 
that  there  may  exist  a  threshold  value  in  the  beam  energy  14 
below  which  oscillation  in  the  BWO’s  stops,  no  matter  how 
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Fig.  6.  Locations  of  the  pair  of  roots,  a-f  and  a-,  on  the  complex  k 
plane  for  H  =  65  kcV  and  A  -  0.4  kA,  L  =  23.8  cm,  iV  =  -1  and 
u;/2t  =  23,295  -  t0.092.  Locations  of  the  roots  (black  circles)  are  outside 
the  oscillation  region  Im  (w)  >  0  shown  by  white  circles.  The  diamond  is 
the  saddle  point 


Fig.  7.  Allowable  oscillation  regions  of  two  solid  curves  for  =  100 
and  65  keV,  i.e.,  horizontal  separation  of  open  circles  as  shown  in  Figs.  5(b) 
and  6  versus  beam  current  /&.  The  dashed  line  is  the  horizontal  separation 
Re  (a-  -  a-f.)  =  0.264  cm*“^  of  the  roots  for  h  =  0,4  kA,  L  =  23.8 
cm,  andiV  =  — 1  in  (9). 


large  the  beam  current  h-  This  statement  is  clarified  in  Fig. 
7,  where  the  fixed  horizontal  separation  between  two  black 
circles  Re  (a-  -  a+)  =  0.264  cm"*^  given  from  (9)  is  shown 
by  a  horizontal  dashed  line.  Two  solid  curves  are  the  horizontal 
separations  Re  (Aib)  between  two  white  circles,  respectively, 
as  shown  in  Hgs.  S(b)  and  6  as  functions  of  h-  Here,  Re  ( A/s) 
is  the  allowable  range  of  oscillatioa  for  a  given  h.  In  the  case 


Fig.  8.  Oscillation  region  versus  beam  energy  \  1.  The  dashed  line  is  the 
horizontal  separation  Re  (a-  -  a+)  =  0.264  for  h  =  0.4  kA. 

L  =  23.8  cm,  and  iV  =  -1.  The  black  circle  is  the  starting  energy. 


of  14  =  100  keV,  the  solid  curve  always  lies  above  the  dashec 
line.  This  fact  suggests  that  the  oscillation  is  possible.  On  the 
other  hand,  in  the  case  of  H  =  65  keV,  the  curve  stays  below 
the  dashed  line,  and  oscillation  does  not  occur. 

For  a  given  current  h  =  0.4  kA,  the  starting  energy  for 
oscillation  is  found  in  Fig.  8,  where  the  horizontal  separations 
Re  (Afc)  of  white  circles  such  as  those  in  Fig.  5(b)  or  t 
versus  14  plotted.  It  is  clearly  shown  that  the  startini 
energy  denoted  by  the  black  circle  is  76.5  keV  for  the  giver 
parameters  of  the  present  large  diameter  SWS.  The  large' 
length  L,  the  smaller  starting  energy,  because  the  dashed  line 
becomes  lower.  It  must  be  emphasized,  however,  that  the 
starting  energy  is  not  a  sufficient  condition,  but  a  necessar^ 
condition  for  initiating  oscillation  in  the  finite  length  SWS.  Ii 
order  to  have  oscillation,  the  pair  of  roots  a+  and  a-  must  alsc 
satisfy  (10)  in  addition  to  (9),  This  additional  condition  yield 
the  starting  current  for  oscillation.  Equation  (10)  includes  th' 
reflection  coefficient  i?  as  a  given  parameter,  and  the  startini 
current  ht  is  affected  sensitively  by  R.  The  results  are  showi 
in  Fig.  9,  where  lat  versus  R  is  calculated  for  14  =  10^ 
keV  and  L  =  23.8  cm.  It  is  reasonably  shown  that  ht  is 
decreasing  function  of  R.  If  we  take  into  account  the  effec 
of  the  fourth  wave  (forward  structure  wave)  as  an  additiona 
feedback  mechanism  in  (8),  an  oscillatory  nature  arises  in  th 
curve  lat  versus  L,  as  was  pointed  out  in  [13],  [16]. 

rv.  Discussion  and  Conclusions 

A  large  diameter  backward  wave  oscillator  of  D/A  4. 
has  been  design^  and  analyzed  numerically  in  detail.  The  ke; 
point  of  the  design  study  is  to  raise  the  oscillation  frequence 
without  decreasing  the  mean  diameter  of  the  SWS.  In  ou 
linear  analysis,  parameter  selection  for  such  purposes  has  beer 
readily  performed. 

It  was  shown  in  the  previous  section  that  there  existed  : 
starting  current  of  the  beam  required  to  withstand  the  leakage 
of  radiation  at  both  ends  of  SWS.  Moreover,  there  existed  ^ 
starting  energy  of  the  incident  beam,  below  which  oscillatioi 
could  not  occur,  because  of  the  finite  length  of  the  SWS.  These 
two  statements  are  never  trivial  and  the  distinction  betweer 
both  conditions  has  not  been  clarified  in  the  literatures  ir 
the  past  The  fonner  condition  for  the  starting  current  giver 
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Fig.  9.  Starting  current  1st  versus  reflection  coefficient  R  for  beam  energy 
=  100  kcV,  L  =  23.8  cm,  and  N  =  -1. 


by  (10)  arises  from  the  requirement  that  the  radiation  loss 
depending  on  h  at  the  ends  of  the  SWS  must  be  smaller  than 
a  threshold  value.  In  other  words»  the  external  Q  value  of 
the  corrugated  wail  resonator  shown  in  Fig.  1  must  be  large 
enough  for  the  occurrence  of  oscillations.  The  resonator  must 
store  a  minimum  amount  of  RF  energy  for  oscillation.  On  the 
other  hand,  the  latter  condition  for  the  starting  energy  given 
by  (9)  arises  from  the  requirement  that  the  coupling  condition 
of  real  wavenumbers.  Re  (a+)  and  Re  (a-.),  between  the 
donor  of  RF  energy  (slow  space  charge  wave)  and  the 
acceptor  a_  (backward  structure  wave)  becomes  stringent, 
when  L  is  small.  The  locations  of  the  a+  and  on  the 
complex  k  plane  are  considerably  separated  in  relatively  short 
L  case,  and  the  coupling  condition,  namely,  phase  matching 
condition  of  the  waves  at  the  ends  may  not  be  satisfied.  This 
restriction.  Re  {Ak),  is  relaxed  in  the  case  of  large  14  as 
was  shown  in  Fig.  8.  This  is  the  reason  for  the  existence  of 
the  minimum  starting  energy  for  oscillation.  This  oscillation 
condition  for  14  is  peculiar  to  the  finite  L  case,  and  such 
additional  requirement  for  the  coupling  of  waves  did  not  occur 
in  conventional  infinite  L  case  shown  in  Fig.  3.  It  is  well 
known  that  beam  wave  interaction  is  strong  near  tt  or  27r 
mode  operation,  where  the  group  velocity  of  the  backward 
structure  wave  is  small.  In  L  finite  case,  the  oscillation  may 
stop  somewhere  between  two  mode  operations  because  of 
increased  group  velocity.  This  qualitatively  well  known  fact  is 
for  the  first  lime  analyzed  quantitatively  in  the  present  paper, 
and  it  is  the  significance  of  the  starting  energy  analyzed  in 
this  paper. 

The  diffraction  of  the  radiation  at  the  ends  of  large  diameter 
SWS  relating  the  value  of  reflection  coefficient  i?  is  a  com¬ 
plicated  problem  to  analyze,  and  it  is  not  treated  here.  In  our 
numerical  analysis,  L  =  23.8,  i?o  =  3.0,  and  A  =  1.25  cm. 
The  structure  length  L  may  be  considered  to  be  longer  enough 
than  /?o  and  A,  and  our  analysis  ignoring  end  diffraction 
problem  may  be  qualitatively  valid.  The  practical  value  of 
the  reflection  coefficient  R  of  the  large  diameter  SWS  can 
be  measured  by  means  of  vector  network  analyzers  for  a  real 
fabricated  sample.  The  coefficient  R  is  estimated  to  be  0.2 
at  20  GHz  from  a  straight  cylindrical  waveguide  assuming 
that  the  corrugation  parameter  a  =  h/R^  =  0.0567  is  small 
enough.  Then,  the  starting  current  =  1  A  is  found  in 
Fig.  9.  In  the  previous  section,  we  assumed  the  typical  beam 


cunent  h  -  0.4  kA,  which  is  400  times  larger  duut  In 
such  cases,  overbunch  instability  may  occur  in  die  oscillation 
and  degradation  of  the  microwave  output  may  result  as  was 
predicted  and  observed  in  [5]  and  [13].  In  order  to  suppress 
the  instability,  we  may  raise  lat  by  decreasing  the  length  L 
under  the  restriction  that  the  pair  of  roots  in  (9)  remain  in 
the  oscillation  region  shown  in  Fig.  5(b).  On  the  other  hand, 
space  charge  limiting  current  in  a  cylindrical  pipe  for  the 
average  radius  i2o  =  3  cm  and  beam  radius  R^  ^  2.68 
cm  is  calculated  to  be  17(7^/^  -  l)^/^/2  In  (iZ^/iio)  = 
2.9  kA,  which  is  much  larger  than  the  present  A,  where 
7  is  the  relativistic  factor.  The  nonlinear  analysis  [22]  or 
numerical  simulation  [13]  is  required  to  predict  the  power 
level  and  the  performances  of  the  designed  large  diameter 
BWO. 

The  energy  sources  in  gyrotrons  are  the  beam  velocityvx 
perpendicular  to  the  axial  magnetic  field.  Electron  beams  with 
=  L5  ~  2  are  usually  required  for  high  efficiency 
operation  in  gyrotrons  [2],  [3].  For  that  purpose,  a  sophisticated 
technology  to  generate  such  spiral  beams  was  devised  by 
means  of  magnetron  injection  gim,  and  analytically  and  em¬ 
pirically  adjusted  nonuniform  axial  magnetic  field  profiles  for 
optimized  performance.  In  contrast  to  gyrotrons,  BWO’s  are 
easy  to  operate,  because  beams  with  vx/t^n  1  are  available, 
and  uniform  axial  magnetic  fields  are  usually  applicable. 

The  growth  rates  of  the  large  diameter  BWO  analyzed  in 
the  present  paper  are  smaller  than  those  in  the  conventional 
BWO’s.  This  is  mainly  because  we  have  chosen  small  a  = 
h/Ro.  In  fact,  we  carried  out  a  preliminary  measurement 
of  BWO  operation  of  our  designed  large  diameter  SWS.  It 
was  found  that  the  microwave  output  at  21  GHz  was  quite 
small  and  that  the  operation  was  made  in  linfear  regime.  It 
may  be  important  to  explore  a  possibility  of  enhancing  the 
growth  rates  of  radiation  from  the  large  diameter  BWO.  Our 
analysis  has  been  confined  to  the  case  of  very  strong  applied 
magnetic  field.  An  extension  of  the  present  linear  analysis  to 
the  finite  magnetic  field  case,  especially  the  case  u  ^  Q,  is 
very  important  [6].  Here,  Q  is  the  relativistic  electron  cyclotron 
frequency. 

Reference  [23]  pointed  out  that,  in  slow  wave  cyclotron 
devices  with  <  w/k  <  c,  the  beam  energy  can  be  converted 
to  wave  energy  through  a  transformation  from  Vjj  to  ui.  This 
device  was  called  the  slow  wave  electron  cyclotron  maser. 
Although  they  assumed  a  dielectric  loaded  smooth  cylindrical 
configuration,  their  statement  can  be  applied  to  a  metal  wall 
SWS  [24],  [25],  In  fast  wave  devices  such  as  gyrotrons,  the 
normal  Doppler  shifted  beam  cyclotron  wave,  u;  =  /cvn  -f  f2, 
is  used  to  interact  with  structure  TE  modes.  On  the  other 
hand,  the  anomalous  Doppler  shifted  beam  cyclotron  wave, 
uj  =  -  f2,  may  be  available  to  interact  with  structure 

TM  modes  in  our  large  diameter  SWS.  The  growth  rates 
in  the  large  diameter  BWO  in  the  present  paper  may  be 
resonantly  enhanced  by  the  novel  effect  suggested  by  Kho 
and  Lin  [23].  Moreover,  they  showed  in  their  Fig.  6  that  the 
slow  wave  electron  cyclotron  maser  was  more  tolerant  of  beam 
momentum  spread  than  a  fast  wave  device  (CARM)  for  high 
efficiency  operation.  Their  results  are  especially  encouraging 
in  the  case  of  high  beam  current  in  which  beam  qualities  such 
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as  momentum  spread  and  emittance  arc  greatly  inferior  to  the 
beams  generated  from  thermoionic  cathodes  [26]. 

In  conclusion,  the  slow  wave  devices  such  as  MWCG’s  will 
be  a  hopeful  candidate  for  the  purpose  of  generating  muld-MW 
millimeter  microwaves  for  a  variety  of  advanced  applications. 
The  large  diameter  BWO  studied  in  the  present  paper  may  be 
helpful  for  that  purpose. 


References 

[1]  A.  V,  L.  Gaponov-Grekhov  and  V.  L.  Granaistein.  Ed.,  Hi^h  Power 
Microwave  Applications  for  the  2 1st  Century.  Dedham,  MA:  Artcch 
House,  1994. 

[2]  V.  L.  Granaistein  and  I.  Alexeff,  Ed.,  Hi^h-Power  Microwave  Sources. 
Dedham,  MA:  Artcch  House,  1987. 

[3]  J.  N.  Benford  and  J.  A.  Sweglc,  High  Power  Microwaves.  Boston, 
MA:  Artech,  1992. 

[4]  R.  A.  Kchs  et  at,  “A  high  power  backward- wave  oscillator  driven  by 
a  relativistic  electron  beam. ’’  IEEE  Trans.  Plasma  ScL,  voi.  PS-13,  pp. 
559-562,  Dec.  1985. 

[5]  Y.  Carmel  et  ai,  "From  linearity  toward  chaos:  Basic  studies  of 
relativistic  backward-wave  oscillators,"  Phys.  Rev.  Lett.,  vol.  69,  no. 
11,  pp.  1652-1655,  Sept.  1992. 

[6]  A.  Vlasov  et  oL,  "Relativistic  backward-wave  oscillators  operating  near 
cyclotron  resonance,"  Phys.  Fluids  B,  vol.  5,  no.  5,  pp.  1625-1638, 
May  1993. 

[7]  S.  P.  Bugaev  et  ai,  "Investigation  of  a  millimeter-wavelength-range 
relativistic  diffraction  generator.”  IEEE  Trans.  Plasma  Sci..  vol,  18,  pp. 
518-524,  June  1990. 

[S] _ "Relativistic  muitiwave  Cerenkov  generators."  IEEE  Trans. 

Plasma  Sci,  vol.  18,  pp.  525-536,  June  1990, 

[9]  D.  K.  Abe,  "Experimental  studies  of  overmoded  high  power  microwave 
generators,”  Ph.D.  dissertation,  Univ.  Maryland,  College  Park,  1992. 

[10]  Y.  Carmel  et  ai,  "High-power  microwave  generation  by  excitation  of 
a  plasma-filled  rippled  boundary  resonator,"  IEEE  Trans.  Plasma  Sci, 
vol.  18,  pp.  497-506,  June  19W. 

[11]  X.  Zhai  et  at,  “Observation  of  Trivelpiece-Gould  modes  in  a  plasma 
filled  backward  wave  oscillator,”  Phy.s.  Rev.  A,  vol.  45,  no.  12  ,  pp. 
R8336-R8339.-’June  1992. 

[12]  J.  A.  Swegle  et  ai,  "Scaling  studies  and  time-resolved  microwave 
measurements  on  a  relativistic  backward-wave  oscillator,"  IEEE  Trans. 
Plasma  Sci,  vol.  21,  pp.  714-724,  Dec.  1993. 

[13]  B.  Levush  et  ai,  "Theory  of  relativistic  backward-wave  oscillators  with 
end  reflections,”  IEEE  Trans.  Plasma  Sci,  vol.  20,  pp.  263-280,  June 
1992. 

[14]  Y.  Carmel  et  ai.  "Relativistic  plasma  microwave  electronics;  Studies  ot 
high-power  plasma-filled  backward- wave  oscillators,  '  Phvs.  Fluids  B. 
vol.  4,  no.  7,  pp.  2286-2292,  July  1992. 

[15]  J.  A.  Swegle,  "Starting  conditions  for  relativistic  backward  wave  os¬ 
cillators  at  low  currents,"  Phys.  Fluids,  vol.  30.  no.  4,  pp.  1201-1211, 
Apr.  1987. 

[16]  K.  Minami  it  ai,  "Analysis  of  starting  currents  in  a  backward  wave 
oscillator  with  finite  structure  length,"  J.  Phys.  Soc.  Jpn.,  vol.  61,  no. 
10,  pp.  3566-3575,  Oct  1992, 

[17]  M.  R.  Amin  et  oL,  "Analysis  of  the  electromagnetic  waves  in  an 
overmoded  finite  length  slow  wave  structure,"  IEEE  Trans.  Microwave 
Theory  Tech.,  vol.  43,  Apr.  1995. 

[18]  J.  A.  Swegle  et  ai,  "Backward  wave  oscillators  with  rippled  wall 
resonators:  Analytic  theory  and  numerical  simulation,"  Phys.  Fluids, 
vol.  28,  no.  9,  pp.  2882-2894,  Sept  1985. 

[19]  K.  Minami  et  ai,  "Linear  theory  of  electromagnetic  wave  generation  in 
a  plasma-loaded  conugated-wall  resonator,”  IEEE  Trans.  Plasma  Sci, 
vol.  18,  pp.  537-545,  June  1990. 

[20]  M.  M.  Ali  et  ai,  "Linear  analysis  of  a  finite  length  plasma-filled 
backward  wave  oscillators,”  Phy.s.  Fluids  B,  vol.  4,  no.  4,  pp.  1023-1032, 
Apr.  1992. 

[21]  R.  J.  Briggs,  Electron-Stream  Interaction  with  Plasmas.  Cambridge, 
MA:  MIT  Press,  1964. 

[22]  T.  Watanabe  et  ai,  "Effect  of  beam  energy  spread  on  radiation  intensity 
in  a  high-power  backward  wave  oscillator  with  finite  length,"  J.  Phys. 
Soc.  Jpn.  Lett.,  vol.  61,  no.  4,  pp.  1136-1140,  Apr.  1992. 

[23]  T.  H.  Kho  and  A.  T.  Lin,  "Slow-wave  electron  cyclotron  maser,”  Phys. 
Rev.  A,  vol.  38,  no.  6,  pp.  2883-2888.  Sept.  1988. 

[24]  E,  Jerby  and  G.  Bekefi,  “Cyclotron-maser  experiments  in  a  periodic 
waveguide,”  Phys.  Rev.  E,  vol.  48,  pp.  4637-4M1,  Dec.  1993. 


[25]  K.  P.  Maheshwari  et  ai.  "Cherenkov-cyclotron  resonance  interaction 
in  a  high-power  plasma-filled  backward-wave  osciUator."  Phys.  Rev.  A. 
vol.  45,  pp.  5866-5871,  Apr.  1992. 

[26]  W.  W.  Destler  and  S.  K.  Guharay,  Eds.,  High  Brightness  Beam  Jor 
Advanced  Accelerator  Applications.  New  York:  AIP  press,  1992. 


K.  Minami  was  bom  in  Japan  in  1938.  He  received 
the  B.S.  (EE)  degree  at  Nagoya  Institute  of  Tech¬ 
nology  in  1962.  the  M.S.  (EE)  degree  at  the  Tokyo 
Institute  of  Technology  m  1964.  and  the  Ph.  D 
degree  at  Nagoya  University  in  1969. 

Since  1986.  he  has  been  a  Professor  in  the  Electri¬ 
cal  Engineering  Depanment  of  Niigata  University. 
Japan.  His  research  interests  include  the  generation 
of  high-power  microwave  radiation,  and  interactions 
between  powerful  microwaves  and  plasmas. 


K.  Ogura  was  bom  in  Japan  in  1957.  He  receivec 
the  B.Sc.  degree  in  physics  from  Okayama  Univer¬ 
sity  in  1981.  and  the  M.Sc.  and  D.Sc.  degrees  in 
physics  from  Kyoto  University  in  1983  and  1989 
respectively. 

He  is  as  an  Associate  Professor  in  the  Elec 
trica!  and  Electronic  Engineering  Depanment  o: 
Niigata  University.  His  research  interests  include 
high-power  microwave  oscillators  and  their  appli 
cations  to  plasma  heating. 


Y.  Aiba,  photograph  and  biography  not  available  at  the  time  of  publication. 


M.  R.  Amin  was  bom  in  Rangpur.  Bangladesh  i 
1959.  He  received  the  B.Sc.  (EEE)  degree  from  ih 
University  of  Rajshahi  in  1984  and  the  M.Sc.  (EEE 
degree  from  Bangladesh  University  of  Engineenn 
and  Technology,  Dhaka,  Bangladesh  in  1987.  H 
is  pursuing  the  Ph.D.  degree  at  Niigata  Universitv 
Japan. 

He  is  now  on  sabbatical  from  Bangladesh  Instiiui 
of  Technology.  Rajshahi,  where  he  is  an  Assistar 
Professor  in  the  Electrical  and  Electronic  Engi 
^  neering  Department.  His  research  interests  includ 

Lheoretical  and  experimental  investigations  of  high-power  microwave  device; 
and  semiconductor  power  electronic  drives. 

Mr.  Amin  is  a  member  of  the  Institution  of  Engineers  Bangladesh  (lEB 
unA  Phvsical  Rncietv  of  Jaoan. 


X.  D.  Zheng  was  bom  in  Shandong,  China  in  1961 
He  received  the  B.S.  and  the  M.S.  degrees  in  optica 
engineenng  from  Zhejiang  University.  Hangzhou 
China  in  1982  and  1986,  respectively.  He  is  pursu 
ing  the  Ph.D.  degree  at  Niigata  University,  Japan. 

His  research  interests  include  theoretical  and  ex 
pcrimcntal  investigations  of  high-power  microwave 
devices,  and  optical  measurements. 


i3r 


IEEE  TRANSACTIONS  ON  PLASMA  SCIENCE.  VOL.  23.  NO.  Z  APRIL  1995 


T.  Vfmrntht  wm.  boar  ia  Ji^pm  1941.  He  re¬ 
ceived  the  B.S.  degrea  in  phy^  from  the  Tokyo 
Infdciite  of  Technology  in  19^,  and  the  M.S.  and 
Ph.D.  degrees  from  Nagoya  University  in  1966  and 
1969.  respectively. 

He  is  cunently  with  the  National  Institute  for 
Fusion  Science.  Nagoya.  His  research  interests  in¬ 
clude  nonlinear  plasma  dynamics,  and  a  generalized 
scheme  to  solve  numeric^ly  a  variety  of  systems  of 
nonlinear  equations. 


W.  W.  I>filkr(M*84-SM’90-F92).  photograph  and  biography  not  availabte 
at  the  time  of  publication. 


V.  L,  Granatstein  (S*59-M’64-Sm’86-F*92),  for  a  photograph  and  biogra¬ 
phy,  see  page  524  of  the  October  1994  issue  of  this  Transactions. 


Y.  Cannd  (S'6^M*69-^M'90).  for  a  photograph  and  biography,  see  page 
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Analysis  of  the  Electromagnetic  Waves  in  an 
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Abstract — The  electromagnetic  fields  of  the  higher  order  axial 
onant  modes  in  a  slow  wave  structure  are  analyzed  and  found 
I  have  considerably  different  characteristics  from  those  of  the 
aventional  fundamental  mode.  Here,  the  reflections  at  both 
produce  axial  resonant  modes  corresponding  to  axisym- 
ric  transverse  magnetic  (TM)  modes.  The  period  of  field 
dulation  of  some  of  the  higher  order  axial  modes  is  shorter 
that  of  the  usual  mode  in  a  cylindrical  waveguide,  which 
lid  be  of  practical  interest  for  higher  power,  higher  frequency 
Btion  of  backward  wave  oscillators.  A  perturbation  technique 
used  to  ascertain  the  field  distribution  inside  the  resonant 
Bvity,  and  the  numerical  results  thus  obtained  are  compared 
I  some  experimental  data. 


I.  Introduction 

HE  growing  need  for  coherent,  efficient,  and  high  power 
microwaves  has  led  to  the  development  of  a  number 
of  innovative  devices,  including  backward  wave  oscillators 
(BWO’s),  which  are  a  promising  class  of  devices  having 
a  number  of  useful  features,  namely:  high  spectral  purity 
microwave  power,  frequency  tunability,  high  efficiency,  etc. 
[I]-[5].  For  example,  for  an  overmoded  structure,  D/A  >  1, 
1%  frequency  tunability  in  the  frequency  range  5. 2-5.6  GHz 
has  been  reported  [6].  Here,  D  is  the  mean  diameter  of  the 
waveguide  and  A  is  the  wavelength.  Microwave  radiations  on 
the  order  of  I  GW  at  a  frequency  up  to  30  GHz  have  been 
obtained  [7],  [8].  Continuous  efforts  are  being  made  to  enhance 
the  power  and  frequency  level  of  the  devices. 

In  the  slow  wave  devices,  the  interaction  of  electromagnetic 
(EM)  quantities  takes  place  inside  the  slow  wave  structures 
(SWS),  In  order  to  understand  the  physics  of  the  mechanism 
involved,  it  is  necessary  to  analyze  the  SWS  in  a  realistic 
way.  Resonators  are  also  used  in  linear  accelerator  (linac), 
but  their  geometries  are  considerably  different  from  those 
used  in  BWO’s.  The  EM  behaviors  of  the  linac  cavities 
have  been  extensively  studied  by  many  researchers  using 
numerical  computational  techniques.  In  order  to  study  the 
accelerator  cavities  with  complex  geometries,  computational 
codes  such  as  SUPERFISH  [9],  URMEL-T  [10],  etc.,  have 
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been  developed.  These  codes  are  based  on  the  discretization 
techniques  of  Maxwell’s  and  Helmholtz's  equation,  in  general. 
The  computational  time  of  these  codes  depends  on  the  number 
of  mesh  points,  geometry  of  the  cavity,  boundary  conditions, 
and  accuracy  demanded  for  the  particular  purpose.  A  brief 
description  of  the  numerical  codes  generally  used  in  linac 
studies  is  presented  in  [  1 1  ]  by  Cooper  and  Jones.  In  this  work, 
an  attempt  has  been  made  to  explore  the  EM  quantities  of 
a  sinusoidally  corrugated  SWS  typically  used  in  the  BWO 
experiments  in  the  simplest  way,  by  using  analytic  equations 
derived  for  this  purpose.  In  the  BWO’s,  the  finite  axial 
dimension  of  the  SWS  causes  reflections  from  the  ends  and, 
thereby,  quantization  of  the  axial  wavenumber  results  in  axial 
resonant  modes.  The  finite  transverse  dimension  of  the  SWS 
limits  the  power  handling  capability  due  to  the  internal  RF 
breakdown.  In  order  to  overcome  such  a  problem,  one  can 
increase  the  mean  diameter  of  the  SWS,  thus  creating  an 
overmoded  system.  In  general,  many  higher  modes  can  be 
oscillated  in  such  an  overmoded  SWS.  It  is  still  possible  to 
operate  the  BWO  preferentially  at  a  particular  mode  with 
higher  frequency  by  carefully  selecting  the  beam  and  size 
parameters.  This  is  because  the  respective  modes  have  dif¬ 
ferent  starting  current  for  oscillation  with  each  other.  This 
difference  can  be  carefully  used  to  select  the  particular  mode 
for  oscillation.  In  fact,  recent  experiments  have  shown  that 
efficient  and  high  power  output  can  be  achieved  in  such 
devices  operating  in  a  single  mode  [6].  Despite  the  numerous 
studies  on  conventional,  weakly  relativistic  microwave  tubes, 
more  detailed  studies  of  the  resonant  modes  in  the  finite  length 
SWS’s  intended  for  operation  with  highly  relativistic  electron 
beams  are  required.  The  motivation  of  the  present  work  is 
the  recent  interest  in  the  generation  of  high  power  microwave 
radiation  employing  overmoded  slow  wave  systems.  Most  of 
the  analyses  of  such  systems  were  performed  assuming  infinite 
length  systems  or  perfectly  matched  finite  length  systems, 
which  are  far  from  the  actual  experiments.  Moreover,  their 
analyses  have  been  restricted  mainly  to  the  fundamental  mode 
[12]-[18].  In  this  paper,  we  model  a  finite  length  SWS 
consistent  with  real  experiments  and  include  higher  order 
modes  in  our  computations.  Specifically,  we  consider  the 
two  wave  interaction  process  with  100%  round-trip  reflection. 
Under  this  assumption,  our  SWS  becomes  an  SWS  cavity  with 
perfectly  shorted  ends.  Detailed  field  calculations  along  with 
higher  harmonic  analysis  of  the  axial  resonant  modes  in  the 
SWS  have  been  made.  For  the  higher  order  axisymmetric  TM 
modes  {TM^s^s  >1),  some  unusual  and  novel  phenomena 
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R(z)-Ro+hcos  (koz) 


FREQUENCY 


NORMALIZED  WAVENUMBER 
(b) 

Fig.  1.  Schematic  diagram  of  a  spatially  periodic  slow  wave  structure  (SWS) 
showing  the  directions  of  forwardly  propagating  wave  F  and  backwardly 
propagating  wave  B.  (a)  SWS  shorted  at  z  =  0  and  z  ■=  L  and  (b)  locations 
of  the  F  and  B  waves  in  the  dispersion  curve  in  the  first  Brillouin  zone. 


have  been  found  for  the  first  time  regarding  the  amplitude 
of  the  spatial  harmonic  components  of  the  resonant  modes  for 
particular  values  of  the  axial  wavenumber.  Cavity  perturbation 
technique  [19]  has  been  employed  to  calculate  the  resonant 
frequency  shift  of  the  cavity.  Some  of  the  present  numerical 
results  are  compared  to  exp)erimental  and  numerical  results 
obtained  by  using  SUPERFISH  in  [20],  and  are  found  to  be 
in  excellent  agreement. ' 

The  organization  of  the  paper  is  as  follows.  In  Section  II, 
we  present  the  mathematical  formulation  of  the  SWS.  Section 
III  describes  the  numerical  results  of  the  analysis.  The  cavity 
perturbation  technique  and  the  corresponding  numerical  and 
experimental  results  are  given  in  Section  IV.  In  Section  V, 
discussion  and  conclusion  of  the  present  works  are  presented. 

II.  Mathematical  Formulation 

First,  we  consider  an  infinite  length  SWS;  next,  after  im¬ 
posing  additional  axial  boundary  condition  to  the  system,  we 
will  obtain  the  EM  field  quantities  in  a  finite  length  SWS.  As 
depicted  in  Fig.  1(a),  the  SWS  is  assumed  to  be  sinusoidally 
corrugated  in  the  axial  direction  with  radius  R{z)  =  i2o  + 
hcoskoz,  where  fco  =  21: /zq.  Physical  quantities  associated 
with  an  EM  mode  are  represented  by  a  spatial  harmonic  series 
satisfying  Floquet’s  theorem.  For  axisymmetric  TM  modes, 
the  axial  electric  field  can  be  expressed  as  [14],  [16] 

f;  (1) 

n=— oo  '  ' 

where  Jo  is  the  Oth-order  Bessel’s  function  of  the  first  kind, 
x\  =  -  k^),  fc„  =  fc  +  nko,  k  is  the  axial  wavenum¬ 

ber,  and  n  is  an  integer.  For  slow  spatial  harmonic  waves 
with  <  0,  Bessel’s  function  Jq  becomes  the  modified 
Bessel’s  function  Iq.  Although  the  contributions  of  these 


-  n 1 2  0  7^/2  n 

NORMALIZED  WAVENUMBER  kzo 


Fig.  2.  Numerically  obtained  dispersion  relations  of  the  SWS 
for— TT  <  kzQ  <  TT  in  the  first  Brillouin  zone.  The  parameters  of  the 
SWS  are:  Ro  =  1.4499  cm,  zq  =  1.67  cm,  and  /i  =  0.406  cm.  Light  lines 
(u/k  =  c)  are  shown  by  dashed  lines. 


(b) 

Fig.  3.  Electric  field  patterns  of  the  TMo\  mode.  SWS  parameters  are: 
i?o  =  1.499  cm,  ZQ  =  1.67  cm,  and  h  =  0.406  cm.  Arrows  indicate 
the  direction  of  the  electric  field.  The  density  of  the  field  lines  indicates 
the  strength  of  electric  field  qualitatively,  (a)  TMoi(07r/6)  mode  and  (b) 
TA/oi(G7r/6)  mode. 


harmonics  are  substantial  inside  the  deep  corrugation  (r 
i?o  -h  /i),  this  region  is  so  small  that  we  can  still  determine 
the  EM  characteristics  inside  the  SWS  correedy.  The  other 
components,  Er  and  H&y  are  derived  from  E^.  The  dispersion 
relation  is  obtained  from  the  boundary  condition  at  the  wall 
of  the  structure,  i.e.,  the  tangential  component  of  electric  field 
should  be  zero  at  r  =  R{z),  The  mth  spatial  Floquet  harmonic 
components  of  this  boundary  condition  can  be  expressed  as 


^i[{n-m)kz] 


ikn  d  \ 

{uyc^-^kDdzJ 
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AXIAL  DISTANCE  2  (  cm  ) 
(a) 


Fig.  4.  Electric  field  patterns  of  the  T Mo^  mode.  SWS  parameters  are  the 
same  as  in  Fig.  3.  Arrows  indicate  the  direction  of  electric  field.  The  density 
of  the  field  lines  indicates  the  strength  of  the  electric  field  qualitatively,  (a) 
riV/o4(07r/G)  mode  and  (b)  TMo4{0v/6)  mode. 


In  order  to  evaluate  this  integral,  we  used  Taylor  series 
expansion  of  the  Bessel’s  function  around  R  =  Rq,  since 
direct  integration  takes  much  time  for  computation.  The  radial 
boundary  condition  is  imposed  into  the  following  matrix  form: 

[D]  •  [A]  =  0  (2) 

where  [A]  is  a  column  vector  with  elements  An,  and  [D] 
is  a  matrix  of  an  infinite  rank  with  each  element  given  by 
Dmn  =  [1  +  {tI^  TTl)Qn]^mn  WherC 

22q+\n-m\qUg^\n_m\y. 

<J=0 

Qn  =  -  kl)  and  a  =  h/R{).  The  dispersion 

relation  is  determined  from  the  condition  that  (2)  should  have 
nontrivial  solutions,  and  is  given  by 

a;)  =  det[D]  =  0.  (4) 

In  our  practical  calculation,  the  value  of  n  is  limited  to  —4 
<  n  <  4,  and  2q  +  |n  -  m\  <  10  is  chosen  in  (3).  By 
comparing  to  direct  integration,  we  have  confirmed  that  the 
Taylor  expansion  of  the  Bessel’s  function  converges  quite 
rapidly  and  that  the  numerical  errors  are  less  than  1%  for 
the  fundamental  mode  with  the  parameters  later  described. 

In  the  case  of  a  finite  length  SWS,  the  additional  boundary 
condition  at  both  ends  of  the  structure  must  be  included. 
Referring  to  Fig.  1(a),  a  forward  propagating  wave  F  in  the 
Z‘direction  is  reflected  at  z  =  L  and  becomes  a  backwardly 
propagating  wave  B.  The  locations  of  the  waves  on  the 
dispersion  curve  are  shown  in  Fig.  1(b).  Both  F  and  B  waves 
have  to  satisfy  (4).  The  F  wave  propagates  from  the  z  =  0 
boundary.  After  a  round-trip  with  reflection  zx  z  the 

resultant  wave  at  the  2  =  0  boundary  should  be  “single” 
valued.  This  is  a  two  wave  interaction  process,  and  the  axial 
boundary  condition  can  be  expressed  as  [16] 

ll^i(kp-kB)L  ^  ^  (5) 


(a) 


(b) 

Fig.  5.  Relative  magnitude  of  the  amplitude  of  the  Boquct  harmonics  An 
versus  normalized  wavenumber  for  0  <  kzo  <  tt.  (a)  TAfoi  mode  and  (b) 
rMo4  mode. 


where  R  is  the  total  reflection  coefficient  at  the  ends  of 
the  SWS,  and  kp  and  ke  are  the  wavenumbers  of  the  F 
and  B  waves,  respectively.  If  the  cavity  is  lossless  and  the 
ends  are  shorted  with  perfectly  conducting  metal  plates,  then 
R  =  1  and  kp  =  -ks  as  shown  in  Fig.  1(b).  Hence,  from 
(5),  kp  =  Ntt/L,  where  N  is  an  integer.  Except  for  the 
propagating  direction,  the  wave  B  is  the  same  as  the  wave 
F.  They  satisfy  identical  radial  boundary  conditions  and  have 
the  same  energy.  Therefore,  in  (1),  the  relationship  between 
the  coefficients  An  can  be  written  as  A^  =  By  summing 
the  F  and  B  waves,  the  expressions  of  the  EM  fields  in  the 
SWS  cavity  become 


E^{z,  r,  t)  =  e-‘"*  2AnJo(^^A  cos{knz)  (6) 

n=— oo  ^  ' 

Er{z,r,t)  =  IU>e-^*  f; 

n=-oo  n  \  / 

(7) 


—iujt 


He{zy  r,  t)  =  -ieoujRoe 
^  2A 

Xn 


COS  {knZ). 


(8) 


Once  the  dispersion  relation  (4)  is  solved  numerically,  the 
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Fig.  6.  Axial  profiles  of  Floqaci  harmonics  of  field  Ez  at  radial  position  r  =  0  for  TA/oi  and  TMqa  modes.  The  higher  order  harmonics  of  small  amplitudes 
are  not  shown  in  the  figures.  Thick,  solid  curves  represent  the  total  value  of  the  electric  field  Ez.  Thin,  dashed  and  chained  curves  represent  the  Floquet 
harmonics  as  indicated  by  n  in  the  figures,  (a)  TA/oi (Ott/G)  mode,  (b)  TAToi(67r/6)  mode,  (c)  TMoAiOir/G)  mode,  and  (d)  TAfo4(6^/6)  mode. 


relative  magnitudes  between  the  An  factors  are  determined 
from  (2).  With  a  known  set  of  values  of  A„,  the  normalized 
values  of  the  fields  and  He  can  be  calculated  from 

(6),  (7),  and  (8),  respectively. 

III.  Numerical  Results 

In  the  present  analysis,  a  six-period  SWS  having  the 
following  size  parameters  is  considered:  the  mean  ra¬ 
dius  Rq  =  1.499  cm;  the  period  of  corrugation  zq  = 
1.67  cm;  and  the  corrugation  amplitude  h  =  0.406  cm. 
These  dimensions  correspond  to  the  experimental  values 
used  by  our  research  group  at  the  University  of  Mary¬ 
land  [20].  Fig.  2  depicts  dispersion  relations  computed 
from  (4).  The  dashed  lines  in  the  figure  are  the  light 
lines  in  free  space.  The  end  reflections  cause  the  six- 
period  SWS  to  resonate  at  seven  distinct  frequencies 
corresponding  to  the  particular  TM  modes.  For  the  first 
Brillouin  zone  (-tt  <  kzQ  <  tt)  of  the  dispersion 
relation,  the  values  of  the  normalized  wavenumber  kzQ 
which  correspond  to  the  seven  resonant  axial  modes  are 
07r/6,7r/6,27r/6,37r/6,47r/6,57r/6,  and  Gtt/G,  respectively. 
We  will  designate  these  axial  modes  by  TMqs{Ntt/6) 
hereafter. 


A.  Electric  Field  Lines  of  the  Fundamental  Mode 

Using  the  derived  field  equations  in  the  SWS,  the  electric 
field  lines  are  calculated.  The  method  of  computation  is  identi¬ 
cal  to  that  developed  by  Ogura  et  al.  [18],  with  a  modification 
for  the  higher  order  modes.  Examples  of  the  electric  field 
patterns  of  the  TMoi  mode  are  depicted  in  Fig.  3(a)  and  (b), 
respectively,  for  TAfoi(0^/fi)  TMoi(67r/6)  modes.  The 
separations  Ar  between  the  field  lines  in  the  radial  direction 
are  so  chosen  that  Ez  Ar  =constant  at  the  axial  position 
Er  -  0.  The  density  of  the  field  lines  in  the  radial  direction 
represents  the  strength  of  the  electric  field  qualitatively.  For 
the  rMoi(07r/6)  mode  in  Fig.  3(a),  the  electric  field  is  fairly 
uniform  in  the  axial  direction,  and  the  radial  variation  of  the 
field  has  a  maximum  on  the  axis  of  the  structure.  The  axial 
pattern  of  the  fields  has  six  zero  points  (Ez  =  0)  for  the 
TMoi  (Gtt/G)  modes  as  shown  in  Fig.  3(b).  The  periodic  nature 
of  the  field  lines  depicted  in  Fig.  3(b)  can  be  explained  by 
the  contribution  of  the  Floquet  harmonics  for  n  =  0  to  the 
resonant  axial  modes. 

B.  Field  Patterns  of  Higher-Order  Modes 

The  numerical  calculation  of  the  field  patterns  for  the  higher 
order  modes  [TMo,(s  >  1)]  is  complicated.  We  followcd| 
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ft  improved  technique  from  the  case  of  the  fundamental 
ftode  [18],  and  have  tried  to  calculate  the  electric  field 
for  the  rMo2,TMo3,  and  TM04  modes.  Some  of  the 
Jesuits  for  the  TM04  mode  are  presented  in  Fig.  4(a)  and  (b), 
^respectively,  for  the  case  of  TMo4(07r/6)  and  TMo4(67r/6) 
.modes.  The  field  patterns  have  fine  and  peculiar  structures  and 
differ  considerably  from  those  of  the  TMoi  mode  in  Fig.  3. 
Ithe  radial  distribution  of  the  electric  field  lines  in  Fig.  4  is 
^localized  and  separated  into  two  parts:  the  peripheral  and 
Itbt  central  axis  regions.  The  possible  explanations  for  this 
difference  can  be  given  by  the  contribution  of  the  Floquet 
harmonics  involved  in  the  SWS.  This  will  be  discussed  in 
detail  in  the  following  subsection. 

C.  Floquet  Harmonics  of  the  Electromagnetic 
Fields  in  the  Structure 

The  number  of  Floquet  harmonics  to  represent  the  EM  fields 
are  practically  limited  by  the  computation  time  and  the  relative 
magnitude  of  the  amplitude  factor  i4„’s.  The  values  of  An 
differ  from  mode  to  mode  as  shown  in  Fig.  5(a)  and  (b), 
respectively,  for  TMqi  and  TM04  modes.  In  general.  An 
decreases  with  increasing  |n|,  namely,  Aq  within  the  first 
Brillouin  zone  of  the  dispersion  relation  is  the  largest.  In 
Fig.  5(a),  at  kzo  =  0,  the  relation  between  /l„’s  becomes 
\An\  =  I^4_„|  and  |^„|  «  |A_n-i|  at  kzo  =  tt.  At  r  = 
0,  the  nonzero  field  component  is  Ez,  and  it  is  proportional 
to  Ao  +  A^i  exp  [-i{2'K ! zq)z]  -f  Ai  exp  [i(27r/zo)z]  +  •  *  • .  As 
shown  in  Fig.  5(a),  the  contribution  of  the  n  =  0  Floquet 
harmonic  to  Ez  is  predominant,  and  the  amplitudes  of  the 
higher  order  harmonics  are  very  small  except  at  kzo  ^  tt  for 
the  rMoi(67r/6)  mode.  For  this  reason,  the  field  lines  of  the 
rMoi(07r/6)  mode,  as  shown  in  Fig.  3(a),  are  almost  straight 
lines  for  r  <  iio-  On  the  other  hand,  the  field  lines  of  the 
rMoi(67r/6)  mode  are  determined  by  the  n  =  0  and  n  =  -  I 
Floquet  harmonics.  These  two  harmonics  have  the  same  field 
variations  in  the  2-direction  with  a  period  of  2zo  and,  hence, 
we  obtain  the  field  pattern  with  the  periodicity  depicted  in 
Fig.  3(b). 

For  the  TM04  mode  in  Fig.  5(b),  the  values  of  A±i  become 
greater  than  Aq  for  small  values  of  k  within  the  first  Brillouin 
zone  of  the  dispersion  curve.  This  result  is  novel  and  was 
not  known  in  the  previous  works  [14H^b],  [18].  The  effects 
of  this  unusual  behavior  are  observed  in  the  electric  field 
patterns  of  the  structure  shown  in  Fig.  4.  In  Fig.  5(b),  as  k 
increases,  \  An  \  approaches  |A-n-i|  at  kzo  =  tt,  as  is  observed 
for  the  TMoi  mode  in  Fig.  5(a).  The  field  quantities  of  the 
TMo4(07r/G)  mode  in  Fig.  5(b)  are  nearly  proportional  to 
A^iexp[-i(27r/2o)2]-|-Aiexp  [i{2'K / zq)z],  which  is  periodic 
with  period  zq.  This  fact  can  be  seen  from  the  field  patterns 
of  the  !rMo4(07r/6)  mode  as  shown  in  Fig.  4(a).  As  kz^ 
approaches  tt,  the  field  patterns  of  the  TMo4(67r/6)  mode  near 
r  =  0  as  shown  in  Fig.  4(b),  have  a  periodicity  2zo,  which  is 
similar  to  that  for  the  TMoi(67r/6)  mode  in  Fig.  3(b). 

The  harmonic  components  of  Ez{t  =  0,  2)  are  shown  in 
Fig.  6(a)  and  (b)  for  the  TMoi  niode  and  in  Fig.  6(c)  and  (d) 
for  the  TMo4  mode.  The  EM  fields  are  normalized  by  Aq. 
The  thick,  solid  curves  in  the  figures  represent  the  total  value 


TMo4(0;r/6)  MODE 


(a) 

TMo4  {Sn/S)  MODE 


(b) 

TMo4  (6;r/6)  MODE 


(c) 

Fig.  7.  Resonance  frequency  shift  A/  of  the  axial  resonant  modes  for  the 
TMo4  mode  due  to  perturbation  by  a  spherical  metallic  bead,  (a)  Axial 
changes  in  A/  of  the  TMo4{07t f6)  mode  for  the  bead  at  r  =  0;  (b)  radial 
changes  in  A/  of  the  TA/o4(G7r/6)  mode  for  the  bead  at  r  =  A  20.  (c) 
axial  changes  in  A/  of  the  TA/o4(G7r/6)  mode  for  the  bead  at  r  =  0.  Here, 
AT  =1,  2,-  •,  6. 

of  Ez.  The  thin,  dashed  and  chained  curves  in  Fig.  6  express 
Ez  components  at  r  =  0  denoted  by  the  harmonic  number  n 
in  (1),  which  can  be  understood  if  one  compares  to  the  values 
of  An  at  A;2o  =  0  and  kzo  =  tt  in  Fig.  5(a)  and  (b).  The 
amplitudes  of  the  higher  order  harmonics  not  expressed  are 
very  small  compared  to  those  shown  in  Fig.  6,  The  2-direction 
periodicities  of  Ez — as  depicted  in  Fig.  6(a)-(d) — correspond, 
respectively,  to  those  of  Figs.  3(a)  and  (b),  and  4(a)  and  (b). 
In  Fig.  6(a),  the  contribution  of  the  n  =  ±  \  harmonics  to  Ez 
results  in  a  small  superimposed  perturbation  in  the  2-direction. 
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TMoi  (6rr/6)  MODE 


TMoi  (6;r/6)  MODE 


TMo3(3  7r/6)  MODE 


TMo3(3^/6)  MODE 


Fig.  8.  Comparison  of  numerical  and  experimental  results  of  frequency  shift  A/  by  the  displaced  bead  in  axial  and  radial  directions.  The  solid  curves 
indicate  the  numerical  results,  and  the  open  circles  represent  the  corresponding  experimental  results.  The  axial  data  were  measured  by  the  bead  at  a  radial 
position  r  =  0.543  cm.  (a)  Axial  changes  in  A/  for  the  TiV/oi (Ctt/G)  mode;  (b)  radial  changes  in  A/  for  the  rMoi(G7r/G)  mode;  (c)  axial  changes 
in  A/  for  the  TMo3{3tt/0)  mode;  (d)  radial  changes  in  A/  for  the  rA/o3{37r/G)  mode. 


The  effect  of  this  perturbation  can  easily  be  observed  in  the 
field  lines  for  the  TMoi(07r/6)  mode  in  Fig.  3(a). 

The  above  results  for  the  TMqa  mode  are  never  trivial  and 
have  not  been  clarified  in  the  past  in  the  literature  [14]-[I8]. 
In  usual  cases,  the  fundamental  Aq  term  is  dominant  in  (1)  as 
shown  in  Fig.  6(b)  and  (d);  the  shortest  axial  period  of  Ez  is 
2zo  for  kzQ  =  TT,  as  was  shown  by  Fig.  3(b)  and  Fig.  4(b). 
On  the  other  hand,  in  the  special  case  where  |i4±i|  is  larger 
than  Aq,  as  shown  in  Fig.  5(b),  the  axial  period  can  become 
as  small  as  zq  for  kzQ  ^  0,  which  is  much  shorter  than  the 
usual  minimum  period  of  2zo, 


IV.  PERTURBATION  TECHNIQUE 

The  perturbation  technique  is  a  powerful  method  to  measure 
EM  field  variations  inside  the  resonant  cavity  [19].  This 
technique  has  been  employed  in  the  following  analysis  of 
SWS  cavity,  and  the  numerical  results  have  been  compared 
to  experimental  ones  obtained  at  the  University  of  Maryland 
[20]. 

When  there  exists  a  small  metallic  sphere  bead  with  radius 
tq  in  the  cavity,  the  EM  fields  in  the  cavity  are  perturbed  and 
the  resonance  frequency  of  the  cavity  changes  by  an  amount  of 
A/.  If  the  perturbation  is  small  enough.  A/  can  be  calculated 
approximately  with  the  unperturbed  field  quantities  Ez,Er, 


and  He  given  by  (6)-(8).  The  resultant  expression  for  A/ 
becomes  [19],  [20] 


^  _  ^7rrg(i/iQ|//.|^-36o(|^,P  +  |£;.P)) 
/o  Jy{fioW  +  eo{\Ez\^-^\ErndV)  • 


From  (9),  the  field  quantities  are  closely  related  to  the  fre¬ 
quency  shift  of  the  SWS  cavity  due  to  the  bead.  In  the  case  of 
a  spherical  bead,  the  EM  field  quantities  cannot  be  determined 
separately,  because  an  unidirectional  perturbation  is  impossible 
to  achieve  with  the  spherical  bead  perturber.  The  perturbation 
in  E  is  always  accompanied  by  a  perturbation  in  H,  i 

The  calculated  values  of  the  frequency  shift  A/  for  the 
TMo4  mode,  with  a  bead  radius  of  0.1195  cm,  are  presented 
in  Fig.  7,  The  axial  changes  in  A/  in  Fig.  7  are  calculate4 
with  the  perturber  at  r  =  0.  The  axial  changes  in  A/  for 
the  TMo4(07r/6)  mode  are  shown  in  Fig.  7(a).  The  radial  and 
axial  changes  in  A/  for  the  TMo4(67r/6)  mode  are  shown 
in  Fig.  7(b)  and  (c),  respectively,  /d  r  =  0,  Er  =  He 
0  and  the  only  nonzero  field  component  in  (9)  is  Ez,  and* 
hence.  A/  in  the  axial  direction  as  shown  in  Fig.  7(a)  am 
(c),  respectively,  for  TMo4(07r/6)  and  TMo4(67r/6)  model 
are  directly  proportional  to  In  other  words,  the  results 

presented  in  Fig.  7(a)  and  (c)  are  proportional  to  the  sqm 
of  the  axial  profiles  of  total  Ez  in  Fig.  6(c)  and  (d).  Hem 
at  r  =  0,  we  can  obtain  the  quantitative  information  al 
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the  axial  electric  field  in  the  SWS  from  the  A/  data.  It 
should  be  noted  that  the  solid  curve  in  Fig.  6(c)  is  not 
symmetric  regarding  E,  (r  =  0)=  0.  The  curve  in  Fig.  7(a)  is 
therefore  not  simply  sinusoidal.  For  r  ^  0,  however,  using  the 
spherical  perturber,  it  is  impossible  to  measure  quantitatively 
the  individual  fields  involved  in  A/.  The  radial  variations  of 
A/  are  dependent  on  and  Ha.  In  the  case  of  r  / 

0,  at  the  axial  positions  where  Er  =  0,  both  contributions 
from  E.  and  to  A /  cannot  be  ignored  as  is  seen  from 
(9).  The  positive  region  in  the  radial  profile  ot  A/  for  the 
TMo4(67r/6)  mode,  in  Fig,  7(b),  is  ascribed  to  Hb>  The  axial 
variations  of  A/  for  the  TM(M(07r/G)  mode,  in  Fig.  7(a),  are 
very  rapid  compared  to  those  for  the  7'A7oi{G7r/G)  mode,  as 
shown  in  Fig.  8(a)  by  the  solid  curve.  Experimentally,  it  is 
difficult  to  measure  axial  changes  in  A/  for  the  TMo4(07r/6) 
mode,  in  Fig.  7(a),  because  of  the  very  rapid  variation  of  the 
fields  in  the  axial  direction.  However,  the  TMo4(G7r/G)  mode 
in  Fig.  7(c)  has  variations  of  A/  similar  to  those  observed 
for  the  TM()i(G7r/G)  mode  in  Fig.  8(a).  These  modes  have 
field  variations  which  are  not  very  rapid  in  the  axial  direction 
as  depicted  in  Fig.  6(b)  and  (d).  The  numerical  results  of 
A/  have  been  compared  to  experimental  measurements  [20], 
and  some  of  these  results  are  shown  in  Fig,  8.  The  axial  and 
radial  changes  in  A/  for  the  TMq i(G7r/G)  mode  are  shown  in 
Fig,  8(a)  and  (b),  respectively.  The  axial  and  radial  changes  in 
A/  for  the  rM()3(37r/6)  mode  are  shown  in  Fig.  8(c)  and  (d), 
respectively.  The  open  circles  represent  the  experimental  data, 
and  the  solid  curves  are  the  present  numerical  results.  The 
numerical  data  presented  in  Fig.  8(a)  and  8(c)  are  tor  the  bead 
at  a  radial  position  r  -  0.543  cm  from  the  axis.  By  comparing 
the  frequency  shift  to  the  field  profiles,  one  can  determine 
the  axial  resonant  modes  in  the  SWS  to  be  7’M()i  (Gtt/G)  and 
rMo3(37r/G)  modes.  The  agreement  between  the  numerical 
and  experimental  results  can  be  clearly  estimated  from  the 
figures.  The  discrepancy  between  them  lies  within  the  range 
of  10-20%. 

V.  Discussion  and  Conclusion 

We  have  numerically  analyzed  the  EM  quantities  of  fun¬ 
damental  and  higher  order  axial  TM  modes  in  a  finite  length 
SWS.  It  is  found  that,  for  the  higher  order  modes,  the  am¬ 
plitudes  of  the  Floquet  harmonics  show  unusual  behavior  as 
depicted  in  Fig.  5(b)  for  the  TM04  mode.  For  the  TMo4(B7r/6) 
mode  (/c2o  =  0)»  the  field  quantities  are  mainly  determined  by 
the  Floquet  harmonics  with  n  =  -1  and  I,  and  the  fields 
have  an  unusual  short  period  of  zq  as  shown  in  Figs.  4(a) 
and  6(c).  This  indicates  that  the  period  of  field  modulations 
can  decrease  to  zo,  which  is  small  compared  to  that  for  the 
usual  modes  in  cylindrical  waveguides  where  the  period  is 
larger  than  2zq.  Such  a  mode  with  short  field  modulations  may 
become  important  for  higher  frequency  operation  of  BWO’s. 

The  numerical  results  presented  in  this  paper  are  being 
verified  experimentally  at  the  University  of  Maryland  [20]. 
Generally  speaking,  it  is  impossible  to  reconstruct  the  field 
distributions  including  phase  change  by  the  frequency  shift 
measurements  only,  because  the  frequency  shift  is  related  to 
only  the  absolute  values  of  the  EM  fields,  as  is  seen  from 


(9).  In  some  cases,  however,  the  experimental  results  can  be 
compared  to  the  numerical  calculations  as  shown  in  Fig.  8.  It  is 
concluded  that  the  measured  resonant  modes  are  TMoi(G7r/6) 
in  Fig.  8(a)  and  (b),  and  TMo3(37r/G)  in  Fig.  8(c)  and  (d). 
The  agreement  between  the  numerical  and  experimental  data 
is  quite  satisfactory. 

To  determine  the  dispersion  characteristics  of  an  SWS 
cavity,  shorted  plates  are  placed  at  both  ends  and  the  resonant 
modes  are  excited  by  a  suitable  mode  launcher  at  one  of  the 
end  plates.  The  degree  of  coupling  between  the  mode  launcher 
and  the  SWS  determines  the  type  of  mode  launcher  to  be 
used  in  the  experiment.  If  the  coupling  between  the  mode 
launcher  and  the  cavity  is  strong,  complete  reflection  on  the 
input  antenna  side  cannot  be  expected.  For  a  disc-type  mode 
launcher,  the  reflection  coefficient  at  the  input  end  will  be 
small,  and  consequently  our  assumption  of  perfectly  shorted 
ends  of  the  SWS  becomes  inapplicable.  Hence,  we  prefer  a 
short  wire  antenna  at  the  center  of  the  plate  as  a  mode  launcher 
to  excite  the  cavity.  However,  it  is  difficult  to  excite  the  surface 
wave  modes  near  kzo  =  tt,  as  was  reported  in  [20]. 
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